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Abstract 

Hyperspectral  imaging  systems  passively  sense  radiant  electromagnetic  energy  from  a  remote 
scene  to  form  a  three  dimension  profile  of  the  remote  scene.  The  data  contained  in  this 
profile  describes  real  images  of  the  remote  scene  (two  spatial  dimensions)  for  a  certain  number  of 
spectral  wavelength  bands  across  a  finite  spectral  range  of  electromagnetic  radiation  (one  spectral 
dimension).  Typical  grating  type  hyperspectral  imaging  systems  collect  spectral  electromagnetic 
radiation  in  the  visible  and  near  infrared  spectral  range,  by  incrementally  scanning  across  the 
spatial  extent  of  the  remote  scene.  Due  to  the  very  low  optical  throughput  used  in  these  systems, 
a  legacy  of  the  optical  techniques  employed  to  scan  for  spectral  and  spatial  information,  adapting 
these  systems  for  spectral  ranges  extending  into  the  mid  wave  and  long  wave  infrared  is  difficult. 
Utilising  an  optical  arrangement  that  delivers  significantly  greater  optical  throughput  for  collecting 
radiant  electromagnetic  energy  from  remote  scenes,  the  chromotomographic  imaging  system  has 
been  reported  to  demonstrate  reasonable  quality  hyperspectral  images  for  finite  spectral  ranges  in 
the  visible,  mid  wave,  and  long  wave  infrared. 

This  thesis  investigates  the  design  for  a  software  model  that  emulates  the  attributes  of  the 
optics  system  loosely  detailed  in  reports  of  the  chromotomographic  imaging  system.  The  optical 
system  angularly  multiplexes  the  spatial  electromagnetic  radiation  collected  from  a  remote  scene 
onto  a  focal  plane  array  detector  using  a  rotating  prism  to  disperse  the  electromagnetic  radiation. 
Using  optics  principles  and  wave  optics  propagation  theory,  the  software  model  has  been  designed 
to  simulate  the  point  spread  function  characteristics  of  the  optics  system  in  the  chromotomographic 
real  image  plane.  The  software  model  can  be  used  as  an  aid  to  develop  and  test  the  characteris¬ 
tics  of  an  optics  system  through  simulation,  or  simulate  acquiring  actual  chromotomographic  real 
image  based  on  spectral  information  about  a  remote  scene,  or  explore  certain  hyperspectral  image 
reconstruction  techniques  suited  to  the  chromotomograhic  imaging  system. 
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SIMULATING  A  CHROMOTOMOGRAPHIC  SENSOR  FOR 
HYPERSPECTRAL  IMAGING  IN  THE  INFRARED 

I.  Introduction 

Interest  in  recording  visual  events  as  images  continues  to  evolve  as  advances  in  various  tech¬ 
nologies  and  the  application  of  these  technologies  continue  to  progress  according  to  the  pace 
of  changing  demands  and  desires  of  human  societies  around  the  world.  Photography  is  recognised 
as  the  first  advance  in  technology  by  mankind  for  directly  recording  real  world  scenery.  The  ap¬ 
plication  of  photography  over  the  years  has  influenced  innovation  and  inventions  with  commercial 
investment  to  explore  different  methods  for  essentially  recording  real  world  scenery.  In  the  engi¬ 
neering  discipline  of  remote  sensing,  photography  has  served  as  the  foundation  for  advancements 
that  have  been  made  in  remote  sensing  to  record  images  of  real  world  scenery. 

As  an  engineering  discipline,  remote  sensing  can  be  viewed  as  applying  scientific  techniques 
for  an  observer  to  gather  useful  information  about  a  remotely  located  subject.  In  this  context, 
photography  is  certainly  a  form  of  remote  sensing,  where  the  useful  information  acquired  from 
photography  is  typically  a  two-dimension  image  of  the  subject  remotely  located  within  a  line  of  sight 
from  the  observer.  One  particular  remote  sensing  application  that  continues  to  realise  potential  for 
further  advancement  is  the  field  of  hyperspectral  imagery  systems. 

1.1  Hyperspectral  Imaging 

Aerial  photography  is  probably  the  best  known  application  of  remote  sensing  where  hyper¬ 
spectral  imagery  systems  have  been  developed  for  the  purpose  of  capturing  and  recording  two 
dimension  images.  Individually  each  image  covers  different  finite  spectral  bands,  and  collectively 
form  a  data  set  of  real  images  across  a  finite  range  of  the  electromagnetic  spectrum.  This  data 
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Linear  dauctor  array 


Figure  1.1:  Hyperspectral  remote  sensing  techniques.  Aerial  scanning  techniques  that  employ 
dispersive  optical  systems  to  record  the  spectral  details  of  each  spatial  sample  of  the  remote  scene 
are  the  (a)  whiskbroom  (across  track)  [3:  311]  and  (b)  pushbroom  (along  track)  [3:  315].  These  tech¬ 
niques  incrementally  scan  across  the  entire  spatial  extent  of  the  remote  scene  over  time,  capturing 
the  hyperspectral  details  for  the  remote  scene. 


set  captures  the  spectral  and  spatial  fidelity  of  information  about  the  remotely  located  subject  or 
scene. 

During  the  process  of  sensing  and  recording  data  from  the  remote  scene,  the  collected  spatial 
and  spectral  details  must  be  in  a  format  at  the  time  of  recording  that  enables  spectral  images  of 
the  remote  scene  to  be  constructed  from  the  recorded  data.  The  images  of  the  same  remote  scene 
collected  over  a  spectral  range  of  the  electromagnetic  spectrum  collectively  form  a  hyperspectral 
data  set  for  the  remote  scene. 

During  the  past  development  of  hyperspectral  imaging  systems,  certain  techniques  have  en¬ 
dured  success  with  sensing  multispectral  remote  scenes  and  representing  this  information  as  hyper¬ 
spectral  data  sets.  The  whiskbroom  (across  track)  and  the  pushbroom  (along  track)  techniques  for 
collecting  multispectral  remote  scene  data  are  two  better  known  approaches  employed  for  hyper¬ 
spectral  image  remote  sensing.  Diagrams  representing  each  technique  appear  in  Figure  1.1. 
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These  techniques  scan  piecewise  subareas  from  a  large  area,  progressively  collecting  the  com¬ 
plete  spatial  extent  of  the  large  area  to  steadily  build  a  spectral  information  set  for  the  large  area. 
The  collection  of  hyperspectral  image  data  with  these  techniques  endures  a  time  penalty,  the  du¬ 
ration  of  time  required  to  scan  the  large  area  and  appropriately  record  the  spectral  information. 
This  time  requirement  has  confined  the  use  of  hyperspectral  imaging  systems  to  applications  where 
prominent  features  in  the  remote  scenery  do  not  change  over  extended  periods  of  time.  Platforms 
that  are  available  to  undertake  the  role  of  carrying  hyperspectral  imaging  systems  to  collect  this 
form  of  data  are  aircraft  and  satellites.  Consequently,  aerial  surveillance  of  geographic  regions  on 
Earth  has  pioneered  the  application  of  hyperspectral  imaging  systems,  where  community  or  rural 
centres  regions  may  be  surveyed  typically  with  aircraft  platforms,  and  large  regions  of  continents 
or  oceans  surveyed  with  satellite  platforms. 

1.2  Airborne  Visible  /  Infrared  Imaging  Spectrometer 

One  well  known  application  of  employing  a  hyperspectral  imaging  system  is  the  Airborne 
Visible  /  Infrared  Imaging  Spectrometer  (AVIRIS)  project.  The  AVIRIS  project  continues  to  be 
undertaken  by  NASA’s  Jet  Propulsion  Laboratory  (JPL)  to  identify,  measure  and  monitor  major 
constituents  of  the  Earth’s  surface  and  atmosphere  based  on  molecular  absorption  and  particle 
scattering  signatures. 

AVIRIS  uses  an  imaging  spectrometer  to  obtain  spectral  information  of  the  remote  scene,  by 
scanning  subregions  as  piecewise  components  of  the  entire  region  contained  in  the  remote  scene. 
Spectral  information  for  each  subregion  is  recorded  during  the  scan  period  allocated  for  each  sub- 
region,  so  progressively  the  hyperspectral  image  data  set  of  the  entire  remote  scene  is  recorded. 
Figure  1.2  is  a  representative  sketch  of  the  spectral  information  this  project  collects  from  remote 
scenes.  Figure  1.3  is  a  sketch  of  a  representative  spectroscopic  profile  for  a  subregion  or  pixel  from 
the  remote  scene.  Figure  1.4  is  a  hyperspectral  image  data  set,  represented  here  as  an  image  cube, 
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AVIRIS  CONCEPT 


EACH  SPATIAL  ELEMENT  HAS  A 
CONTINUOUS  SPECTRUM  THAT 
IS  USED  TO  ANALYZE  THE 
SURFACE  AND  ATMOSPHERE 


224  SPECTRAL  IMAGES 
TAKEN  SIMULTANEOUSLY 


Figure  1.2:  AVIRIS  concept  of  Hyperspectral  Imaging  [1] 


of  Moffett  Field  at  the  southern  end  of  San  Francisco  Bay,  California.  The  false  colour  image  on 
top  of  the  image  cube  has  been  established  to  highlight  the  structure  of  the  water  and  evaporation 
ponds  in  the  lower  right  region  of  the  remote  scene.  The  sides  of  the  cube  represent  spectral  slices  of 
the  respective  edges  of  the  remote  scene  represented  on  top,  across  the  spectrum  of  electromagnetic 
radiation  from  400  nm  (near  ultraviolet)  to  2,500  nm  (near  infrared). 


1.3  Conventional  Hyperspectral  Imaging  Limitation 

The  AVIRIS  approach  with  collecting  multispectral  remote  scene  image  data  highlights  the 
characteristic  limitation  of  information  throughput  with  conventional  hyperspectral  imaging  tech¬ 
niques  that  scan  across  one  dimension,  whether  spatially  or  spectrally,  to  obtain  the  complete  set 
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Figure  1.3:  AVIRIS  sample  profile  of  Intensity  for  a  single  pixel  on  the  detector  [1] 

of  spatial  and  spectral  information  from  the  remote  scene.  The  limitation  of  this  throughput  is 
further  compounded  by  the  time  interval  needed  to  completely  scan  the  remote  scene. 

Adapting  conventional  state-of-the-art  hyperspectral  imaging  systems  for  spectral  ranges  in 
the  mid  wave  and  long  wave  infrared  regions  of  the  electromagnetic  spectrum  is  compounded  by 
the  very  low  optical  throughput  techniques  used  these  systems.  At  infrared  spectral  wavelengths 
features  in  the  remote  scene  characteristically  have  low  contrast  between  energy  levels  of  electro¬ 
magnetic  radiation.  For  remote  sensing  in  the  infrared  and  overcome  the  low  contrast  between 
radiant  energy  levels,  a  hyperspectral  imaging  system  would  require  an  optical  system  arrangement 
that  offered  significantly  higher  optical  throughput. 

The  chromotomographic  imaging  system  was  developed  by  Jonathon  M.  Mooney  and  col¬ 
leagues  of  AFRL/SNHI  at  Hanscom  Air  Force  Base  [4-6],  utilising  an  optical  system  with  high  op¬ 
tical  throughput  that  angularly  multiplexes  the  spatial  electromagnetic  radiation  collected  from  a 


5 


Figure  1.4:  AVIRIS  Hyperspectral  Image  Cube  (Hyperspectral  Data  Set)  [1] 

remote  scene  onto  a  focal  plane  array  detector,  to  provide  significantly  greater  optical  throughput  for 
collecting  radiant  electromagnetic  energy  from  remote  scenes.  Comparisons  of  optical  throughput 
between  conventional  hyperspectral  imaging  systems  and  the  chromotomographic  imaging  system 
are  worth  noting  with  respect  to  the  hyperspectral  data  cube,  which  describes  the  three  dimension 
image  data  remotely  sensed  by  these  systems.  Figure  1.5  shows  this  comparison  for  the  hyperspec¬ 
tral  image  remote  sensing  techniques  of  spatial  scanning,  spectral  filtering,  and  chromotomography. 

1-4  Spectrotomography 

Adapting  the  application  of  hyperspectral  imaging  systems  for  other  uses  is  a  current  challenge 
in  the  remote  sensing  community,  particularly  over  the  past  decade.  Using  hyperspectral  imaging 
systems  for  surveying  features  in  remote  scenes  that  are  characteristically  transient  over  time  is 
the  challenge  being  tackled  as  the  next  advancement  in  the  application  of  hyperspectral  imaging 
systems. 
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Figure  1.5:  Optical  throughput  characteristics  with  conventional  hyperspectral  imaging  systems 
can  be  compared  from  the  hyperspectral  data  cube  where  (a)  spatial  scanning  (b)  spectral  filtering 
methods  deliver  low  optical  throughput.  The  spatial  and  spectral  slices  respectively  demonstrate 
this  throughput.  A  higher  optical  throughput  (c)  chromotomography  is  achieved  with  the  optical 
sensor  used  in  the  chromotomographic  imaging  system  [7]. 

Studies  by  Bulygin  et  al  [8]  and  Okamoto  et  al  [9]  pioneered  the  development  of  techniques 
that  demonstrate  how  multispectral  spatial  data  from  a  remote  scene  can  be  collected  and  recorded 
using  an  optical  system  in  a  manner  suited  to  hyperspectral  imaging.  A  computed  tomography 
approach  is  adopted  to  construct  the  spectral  real  images  of  the  remote  scene  across  the  spectral 
range  of  electromagnetic  radiation  obtained  with  the  optics  system.  Computed  tomography  was 
pioneered  in  the  medical  science  discipline  to  obtain  the  three  dimension  spatial  impression  of  a 
subject  from  the  collection  of  two  dimension  spatial  image  slices  of  the  subject.  Applying  the 
concept  of  computed  tomography  to  hyperspectral  imaging,  the  impression  of  the  subject  (remote 
scene)  is  obtained  as  a  collection  of  two  dimension  spatial  images  slicing  piecewise  through  a  spectral 
range  of  the  electromagnetic  spectrum  in  spectral  wavelength  bands.  This  computed  tomography 
collection  of  spatial  images  is  more  accurately  known  as  the  spectrotomography  collection  of  the 
remote  scene,  analogous  to  a  hyperspectral  data  set  of  images. 

Further  research  over  the  past  decade  in  spectral  imagery  that  employs  the  concept  of  spatially 
separating  spatial  images  of  a  subject  according  to  spectral  content  has  been  reported  in  various 
forums  by  Mooney  et  al  [4-6],  Brodzik  and  Mooney  [10],  Descour  and  Dereniak  et  al  [11,12],  Volin 
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et  al  [13],  Scholl  et  al  [14],  and  Ford  et  al  [15].  These  reported  works  contain  explanations  of 
the  top-level  detail  for  the  actual  systems  of  optics  respectively  used,  but  do  not  contain  sufficient 
detail  that  would  allow  these  efforts  in  spectrotomography  to  be  replicated  easily.  This  list  of 
works  is  by  no  means  an  exhaustive  representation  of  the  research  activity  experienced  in  the  field 
of  spectrotomography. 

1.5  Research  Motivation 

The  Center  for  Measurement  and  Signature  Intelligence  (MASINT)  Studies  and  Research 
(CMSR)  at  the  United  States  Air  Force  Institute  of  Technology  (AFIT)  is  interested  in  acquiring 
a  sensor  capable  of  hyperspectral  imaging  in  the  infrared  region  of  the  electromagnetic  spectrum 
to  capture  and  record  dynamic  real  world  remote  scenes  in  real  time.  A  sensor  that  employs  the 
techniques  of  spectrotomography  can  realistically  meet  this  capability.  The  sensor  could  also  be 
regarded  at  this  time  a  technologically  leading-edge  concept,  and  worthy  of  researching  with  a  goal 
of  designing  a  sensor  suitable  for  use  over  predetermined  spectral  ranges  of  the  electromagnetic 
spectrum. 

Improving  upon  the  MASINT  understanding  of  the  radiant  energy  phenomenology  in  the 
fireball  resulting  from  an  explosion  or  explosions  in  the  remote  scene  is  the  foundation  of  interest 
behind  this  research  of  the  chromotomographic  sensor.  To  capture  the  spectral  and  spatial  details 
from  transient  events  in  the  remote  scene,  such  as  explosions,  in  the  context  of  hyperspectral 
imagery,  a  sensor  with  the  optical  throughput  characteristics  of  the  chromotomographic  sensor  is 
needed. 

The  chromotomographic  sensor’s  high  optical  throughput  feature  offers  the  ability  to  collect 
the  spectral  electromagnetic  radiation  content  from  the  remote  scene  to  form  a  hyperspectral  real 
image  data  set  representing  the  remote  scene  in  a  significantly  lower  duration  of  time  (seconds) 
compared  with  other  grating  type  conventional  scanning  hyperspectral  imaging  systems  (hours). 


This  lower  time  duration  with  the  staring  chromotomographic  sensor  allows  for  speedier  sequential 
collection  over  time  of  radiant  energy  from  the  remote  scene  to  obtain  a  time  resolved  series  of 
hyperspectral  real  image  data  sets  representing  the  remote  scene.  Achieving  time-resolved  hyper- 
spectral  imagery  of  transient  events  elevates  the  dimensionality  of  for  this  type  of  imagery  with  the 
inclusion  of  time  as  a  dimension  to  a  four  dimension  real  image  data  set. 

The  work  reported  by  Mooney  et  al  [4-6]  explored  recording  remote  scenes  in  the  infrared 
range  of  the  electromagnetic  spectrum  as  hyperspectral  image  data  sets.  Recent  descriptions  of 
this  work  refers  to  the  complete  system  of  collecting  the  spatial  and  spectral  information  from  a 
multispectral  remote  scene  for  constructing  the  remote  scene  as  a  hyperspectral  image  data  set 
as  a  Chromotomographic  Imaging  System  (CTIS).  Instead  of  referring  to  spectrotomography,  this 
thesis  refers  to  chromotomography  and  the  Chromotomographic  Imaging  System  in  recognition  of 
the  preceding  work  by  Mooney  et  al. 

To  augment  AFIT’s  research  effort,  this  thesis  reports  the  research  undertaken  with  develop¬ 
ing  a  software  model  that  simulates  the  entire  optics  element  of  the  chromotomographic  imaging 
system.  This  software  model  employs  a  wave  optics  approach  for  representing  the  propagation  of 
electromagnetic  field  wavefronts  through  the  entire  optics  element  of  the  system.  Mooney  et  al 
presented  a  simplistic  description  of  the  optics  element  of  the  CTIS.  This  description  has  been  used 
loosely  as  an  aid  to  developing  the  software  model  reported  in  this  thesis. 

1 . 6  Research  Goals 

The  goal  of  the  research  effort  reported  in  this  thesis  is  to  develop  a  software  model  of  the 
entire  optics  element  of  the  chromotomographic  imaging  system.  The  software  model  is  intended  to 
be  a  software  tool  to  enable  understanding  of  the  optics  element  of  the  CTIS  through  exercising  the 
software  model  for  different  optics  configurations.  Gaining  a  deeper  understanding  and  appreciation 
of  constraints  and  freedoms  that  can  be  utilised  for  certain  design  parameters  in  developing  an 
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optics  element  using  the  software  model  is  intended  to  provide  the  user  a  realistic  expectation  of 
performance  from  a  physical  optics  element  with  equivalent  design  parameters. 

The  nature  of  constructing  hyperspectral  image  data  sets  with  the  post-processing  element 
of  the  chromotomographic  imaging  system  relies  on  obtaining  the  point  spread  functions  related 
to  the  collected  hyperspectral  data  set  of  the  optics  element  of  the  system.  With  an  established 
optic  element  baseline  of  design,  the  software  model  may  be  used  to  provide  this  information  for 
the  purpose  of  aiding  with  constructing  hyperspectral  image  data  sets. 

Furthermore  with  an  established  optic  element  baseline  of  design,  the  software  model  could 
be  utilised  to  obtain  simulated  a  priori  information  regarding  the  remote  scene  as  an  aid  for  recon¬ 
struction.  The  accuracy  of  this  information  would  depend  on  the  accuracy  to  which  the  software 
model  could  simulate  the  behaviour  of  the  real-world  sensor. 

1 . 7  Organisation 

Chapter  2  introduces  the  chromotomographic  imaging  system  chosen  for  the  basis  of  develop¬ 
ing  the  software  model  of  the  system  and  discusses  some  of  the  fundamentals  with  optics  considered 
in  the  synthesis  of  developing  the  software  model. 

Chapter  3  introduces  the  specifics  of  the  simulation  model,  separating  the  entire  sensor  op¬ 
tics  system  situated  in  an  environment  as  a  defined  system,  into  manageable  components  of  the 
simulation  model,  explaining  the  choices  and  decisions  made  with  developing  each  component  of 
the  model. 

Chapter  4  discusses  the  software  model,  examines  the  various  parameters  used  in  the  software 
model,  explains  the  approach  undertaken  to  evaluate  these  parameters  sufficiently  for  a  sensible 
configuration  of  parameters,  so  these  parameters  can  be  used  to  simulate  the  chromotomographic 
sensor  that  comprises  the  sensor  optics  system  of  the  chromotomographic  imaging  system.  The 
software  model  simulates  the  spectral  point  spread  functions  for  the  chromotomographic  sensor 
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for  a  scenario  configuration  that  comprises  the  sensor  optics  system  and  the  separation  distance 
between  the  sensor  and  the  remote  scene.  Evaluation  of  a  scenario  configuration  is  presented  to 
give  an  insight  about  the  type  of  information  the  software  model  provides  toward  understanding 
the  simulated  chromotomographic  sensor  system. 

Chapter  5  reviews  the  attributes  reported  in  this  thesis  that  contributed  during  the  research 
effort  for  designing  a  software  model  to  simulate  the  chromotomographic  sensor  system. 


11 


II.  Background 


This  chapter  introduces  the  chromotomographic  imaging  system  (CTIS)  as  presented  by 
Mooney  et  al  [4-6] .  The  research  presented  in  this  thesis  used  the  description  of  the  CTIS 
system  to  loosely  guide  the  development  of  the  software  model.  Geometric  optics  and  Fourier 
optics  are  introduced  as  alternate  approaches  to  modelling  the  optics  component  of  the  CTIS. 
Understanding  where  each  approach  has  relevance  in  the  design  of  the  software  model  is  important. 
Introducing  the  concept  of  propagating  light  as  a  spatial  and  spectral  waveform,  the  genesis  of  wave 
optics,  as  the  backbone  of  the  software  model. 

2.1  Chromotomographic  Imaging  System 

The  concept  of  the  CTIS  used  as  the  guide  for  developing  the  software  model  was  described 
originally  by  Mooney  et  al  [4]  as  an  Angularly  Multiplexed  Spectral  Imager  (AMSI) .  An  impression 
of  the  optical  system  arrangement  used  in  the  CTIS  is  shown  at  Figure  2.1. 


Figure  2.1:  Optical  system  arrangement  for  the  chromotomographic  imaging  system  [6] 

The  CTIS  as  a  complete  system  can  be  viewed  to  contain  two  sub-system  components,  the 
chromotomographic  sensor  system  and  the  digital  post-processing  and  image  construction  system, 
which  stand  alone  separately  as  significant  systems.  Previously  these  systems  were  discussed  as 
the  optics  component  and  the  digital  post-processing  component  of  the  CTIS  respectively.  The 
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Figure  2.2:  Arrangement  of  projection  planes  for  the  chromotomographic  imaging  system  [7, 10] 

separation  between  these  systems  is  described  in  Figure  2.2  from  the  viewpoint  of  projection  planes 
that  characterise  key  stages  of  the  information  process  for  the  chromotomographic  imaging  system. 


Although  the  process  of  achieving  hyperspectral  imaging  from  the  CTIS  has  been  described  [4- 
7,10],  the  science  used  for  designing  a  chromotomographic  sensor  system  to  achieve  image  spec¬ 
troscopy  by  spatially  separating  spectral  images  in  the  chromotomographic  real  image  plane  is  not 
clearly  outlined.  Research  reported  in  this  thesis  was  focussed  on  understanding  the  science  of  wave 
optics  propagation  and  establishing  a  software  model  that  accurately  represents  the  performance 
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of  a  physically  equivalent  chromotomographic  sensor  system.  In  the  software  model  design,  the 
sensor  optics  system  represents  the  chromotomographic  sensor  system. 

2.1.1  Description.  The  chromotomographic  sensor  system  is  comprised  of  a  double  lens 
arrangement  separated  with  a  field  stop  aperture,  followed  with  a  prism,  a  detector  focussing 
lens  and  focal  plane  array.  Figure  2.1  details  this  arrangement  where  the  double  lens  and  field 
stop  arrangement  collect  the  multispectral  electromagnetic  radiation  and  spatial  features  from  a 
remotely  located  subject.  The  electromagnetic  radiation  that  propagates  beyond  the  second  lens 
is  expected  to  be  collimated. 

These  collimated  wavefronts  enter  the  prism,  called  the  direct  vision  prism,  which  is  designed 
for  the  purpose  of  dispersing  the  spectral  wavefront  fields  contained  in  the  multispectral  electro¬ 
magnetic  radiation  in  angular  directions  with  respect  to  the  optic  axis.  The  design  of  the  prism 
affords  for  a  certain  wavelength  contained  in  the  finite  spectral  range  to  propagate  through  the 
prism  to  the  detector  focussing  lens  undeviated.  The  reference  to  undeviated  means  the  spectral 
field  for  the  specific  undeviated  wavelength  should  propagate  from  the  prism  to  the  focussing  lens 
parallel  to  the  optic  axis,  thus  exhibiting  no  net  deviation  in  angular  direction  of  propagation  with 
respect  to  the  optic  axis  after  propagating  through  the  prism.  All  other  spectral  wavefront  fields 
propagate  from  the  prism  adjacent  to  the  optic  axis  with  some  net  angular  change  in  direct  of 
propagation  with  respect  to  the  optic  axis. 

The  dispersed  spectral  wavefront  fields  propagate  from  the  prism  as  incident  radiation  on 
the  focussing  lens  of  the  detector.  Propagating  through  the  focussing  lens,  the  spectral  wavefront 
fields  are  focussed  at  the  detector  focal  plane  array  as  real  images  of  the  remote  scene,  separated 
spatially  according  to  the  spectral  wavelength  of  intensity  contained  in  each  real  image.  The  spatial 
separation  between  real  images  is  governed  by  the  angular  dispersion  encountered  from  propagation 
through  the  direct  vision  prism. 
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For  the  task  of  constructing  the  spectral  real  images  from  the  spatial  intensity  data  available 
at  the  focal  plane  array,  the  CTIS  uses  the  point  spread  function  (PSF)  characteristics  of  the 
chromotomographic  sensor  system.  The  PSF  for  each  spectral  wavelength  band  represented  with  a 
real  image  differ  in  location  on  the  focal  plane  array  for  each  angular  position  of  orientation  for  the 
prism.  At  the  chromotomographic  projection  plane,  the  point  spread  functions  for  each  spectral 
wavelength  band  that  uniquely  characterise  the  chromotomographic  sensor  system  are  dependent 
on:  (1)  the  spectral  wavelength  of  electromagnetic  radiation  for  spatial  extent  and  intensity  level; 
and  (2)  the  net  angular  direction  of  propagation  and  prism  angular  position  of  rotation  for  spatial 
location  in  the  chromotomographic  projection  plane. 

To  form  the  hyperspectral  image  data  cube,  the  CTIS  separately  convolves  each  spatial  point 
spread  function  with  the  spatial  characteristics  of  the  remote  scene  at  each  spectral  wavelength. 
For  each  angular  position  of  prism  rotation,  summing  the  outcomes  from  each  two  dimension 
convolution  operation  according  to  the  spatial  locations  of  the  point  spread  functions  in  the  chro¬ 
motomographic  projection  plane  results  in  real  images  of  the  remote  scene  that  are  spatially  offset 
according  to  the  spectral  wavelength  band  for  each  point  spread  function.  These  real  images  ar¬ 
ranged  spatially  represent  the  chromotomographic  real  image  (spectral  real  image  data)  for  a  single 
angular  position  of  prism  rotation  (frame  of  data  collection)  recorded  on  the  detector  focal  plane 
array.  Illustrating  this  explanation,  the  mathematical  expressions  used  to  describe  the  CTIS  [4] 
have  been  adapted  here  to  represent  four  spectral  wavelength  bands  and  four  angular  positions  of 
rotation  for  the  direct  vision  prism,  which  correspond  to  four  frames  of  chromotomographic  real 
images.  Equations  2.1  to  2.4  describe  this  relationship. 

ri(x,y)  =  Pi,i(x,  y)  (g)Ci(x,y)  +  pi,2(x,y)  <g>  c2(x,  y)  +  pij3(x,  y)  <g)  c3(x,  y)  +  pi,4(x,  y)  ®  c4(x,  y) 

(2.1) 

r2(x,y)  =  p2,i(x,  y)  ®Ci(x,y)  +  p2,2(x,y)  <g>  c2(x,  y)  +  p2>3(x,  y)  0  c3(x,  y)  +  p2j4(x,  y)  <g  c4(x,  y) 

(2.2) 
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r3(x,y)  =  p3,i(x,y)  0Ci(x,y)  +  p3,2(x,y)  0  c2(x,  y)  +  p3>3(x,  y)  0  c3(x,  y)  +  p3j4(x,  y)  0  c4(x,  y) 

(2.3) 

r4(x,y)  =  p4,i(x,y)  0Ci(x,y)  +  p4,2(x,y)  0  c2(x,  y)  +  p4>3(x,  y)  0  c3(x,  y)  +  p4j4(x,  y)  0  c4(x,  y) 

(2.4) 

Introducing  the  parameter,  pa>b  describes  the  point  spread  function  (denoted,  p)  for  the  spec¬ 
tral  wavelength  band  (denoted,  b)  resulting  from  the  specific  angular  position  of  prism  rotation 
(denoted,  a).  The  parameter,  ra(x,  y)  describes  the  chromotomographic  real  image  frame  recorded 
(denoted,  r),  and  this  real  image  frame  is  linked  to  the  specific  angular  position  of  prism  rotation. 
The  parameter  Cb(x,  y)  describes  the  spatial  details  of  the  remote  scene  (denoted,  c)  with  electro¬ 
magnetic  radiation  in  the  spectral  wavelength  band  (denoted,  b).  The  symbol  0  represents  the 
two-dimension  spatial  convolution  operation. 

The  hyperspectral  image  reconstruction  approach  described  for  the  CTIS  [4]  relies  on  ob¬ 
taining  the  system  transfer  function  matrix  from  the  complete  set  of  point  spread  functions.  The 
inverse  of  the  system  transfer  function  matrix  together  with  the  two-dimension  spatial  Fourier 
Transform  of  each  recorded  frame  of  chromotomographic  real  images  formed  on  the  detector  focal 
plane  array,  enables  the  spectral  real  images  of  the  remote  scene  to  be  computed.  The  outcome  is 
a  hyperspectral  image  data  set  that  spectrally  reconstructs  the  remote  scene. 

A  further  observation  to  note:  each  angular  position  of  prism  rotation  corresponds  to  a 
chromotomographic  real  image  frame  of  recorded  data.  There  must  be  at  least  as  many  frames  (one 
for  each  specific  angular  position  of  prism  rotation)  as  spectral  wavelength  bands  to  be  reconstructed 
for  the  digital  post-processing  approach  to  have  enough  parameters  for  successful  hyperspectral 
image  reconstruction. 

2.1.2  Linear  Shift  and  Time  Invariant  System.  The  CTIS  chromotomographic  sensor 
system  has  the  attributes  of  a  typical  optics  system,  with  lenses  and  a  prism  being  the  only  optical 
components  in  the  system.  These  components  do  not  introduce  electromagnetic  uncertainties  with 
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Figure  2.3:  Linear  System  Interpretation  of  Sensor  Optics  System 

field  propagation  that  would  prohibit  treating  the  system  as  a  linear  shift  and  time  invariant  (LSTI) 
system  [2:  519-556].  Figure  2.3  shows  the  linear  systems  viewpoint  of  the  chromotomographic  sensor 
system  with  the  remote  scene  included  to  completely  characterise  the  remote  sensing  scenario. 

From  a  linear  system  perspective,  each  point  spread  function  represents  a  separate  impulse 
response  characteristic  of  the  chromotomographic  sensor  system.  For  each  specific  angular  position 
of  prism  rotation,  the  corresponding  spectral  point  spread  functions  across  the  finite  spectral  range 
of  electromagnetic  radiation  summed  together  according  to  spatial  location,  forms  a  chromotomo¬ 
graphic  real  image  frame.  This  frame  may  also  be  interpreted  as  the  chromotomographic  point 
spread  function  for  the  angular  position  of  prism  rotation,  or  prism  prism  orientation. 

2.1.3  Geometric  optics  interpretation.  Geometric  optics  is  very  useful  for  interpreting 
the  stages  of  propagation  the  software  model  of  the  sensor  optics  system  is  expected  to  simulate. 
Electromagnetic  radiation  is  treated  as  rays  in  the  direction  of  propagation  with  geometric  optics, 
and  is  also  referred  to  as  ray  optics.  The  paths  of  propagation  for  the  electromagnetic  radiation 
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Figure  2.4:  Geometric  Optics  Interpretation  of  Sensor  Optics  System 

wavefront  are  captured  with  sufficient  accuracy  that  allow  for  the  direction  of  propagation  behaviour 
of  radiation  through  a  system  to  be  readily  identified. 

Rays  of  electromagnetic  radiation  propagating  through  each  stage  of  the  sensor  optics  sys¬ 
tem  can  be  expressed  with  a  matrices  that  characterise  the  propagation  in  terms  of  the  vertical 
displacement  and  angle  subtended  to  the  horizontal  plane  with  the  optic  axis.  More  precise  details 
about  the  propagation  of  electromagnetic  radiation,  like  phase  information  across  the  wavefront  of 
the  electromagnetic  field,  are  not  treated  with  a  geometric  optics  approach. 

Figure  2.4  shows  the  sensor  optics  system  from  the  software  modelling  viewpoint  and  the 
transformation  matrices  for  each  stage  of  propagation.  These  matrices  for  the  different  subsystems 
of  the  sensor  optics  system  are  described  in  Table  A.l  for  the  forward  optics  system,  Table  A. 2  for 
the  prism  optics  system,  and  Table  A. 3  for  the  detector  optics  system. 
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Figure  2.5:  Fourier  Optics  Interpretation  of  Sensor  Optics  System 

Phenomena  with  the  propagation  of  light  as  an  electromagnetic  wavefront,  like  interference 
patterns  in  the  wavefront  along  the  path  of  propagation  clue  to  the  diffraction  of  light  through  an 
aperture  and  the  incoherence  of  light  propagating  as  photon  radiation,  are  not  represented  with 
geometric  optics. 

As  a  technique  to  aid  with  modelling  the  sensor  optics  system,  geometric  optics  is  a  simple 
technique  for  implementing  the  appropriate  distances  between  lenses  and  imaging  planes  in  the 
system.  Establishing  the  transverse  plane  locations  along  the  optic  axis  where  real  images  of  the 
remote  scene  should  be  observed  is  readily  possible  with  use  of  the  Lensmaker’s  formula  for  lenses, 
and  Snell’s  Law  of  Refraction  for  the  prism.  Expressions  for  the  Lensmaker’s  formula  and  Snell’s 
Law  of  Refraction  are  presented  in  section  2.2. 


19 


2.1.4  Fourier  optics  interpretation.  Fourier  optics  is  another  alternate  approach  for 
describing  the  transitions  of  the  electromagnetic  field  for  the  electromagnetic  radiation  wavefront. 
Propagating  through  lenses  and  the  direct  vision  prism,  the  wavefront  experiences  changing  optical 
path  lengths  that  influence  whether  the  wavefront  beam  converges,  diverges,  or  is  collimated  for 
different  stages  of  propagation.  Figure  2.5  represents  the  sensor  optics  system  this  time  described 
from  the  Fourier  optics  viewpoint. 

Fourier  optics  employs  approximations  for  the  optical  path  lengths  for  points  along  the  wave- 
front  to  propagate  between  transverse  planes.  These  approximations  enable  the  stages  of  propaga¬ 
tion  characterised  by  Fourier  optics  to  be  evaluated  with  the  two-dimension  Fourier  Transform  for 
each  spectral  wavelength. 

From  a  software  modelling  viewpoint,  using  the  Fourier  optics  approach  to  model  these  stages 
of  propagation  offers  the  computational  advantage  of  speed  at  the  legacy  of  amount  of  data.  Simula¬ 
tion  packages  like  MatLab®  have  dedicated  fast  Fourier  Transform  algorithms  that  are  significantly 
quicker  than  separate  instructions,  however  using  these  algorithms  for  cascading  stages  of  propa¬ 
gation  can  attract  a  significant  increase  in  generated  electronic  data. 

Using  Fourier  Transforms  to  simulate  the  electric  field  propagation  of  electromagnetic  radia¬ 
tion  generally  requires  the  source  of  radiation  to  be  located  on  and  symmetric  with  respect  to  the 
optic  axis.  Using  Fourier  Transforms  as  the  basis  of  simulating  the  electric  field  propagation  in  a 
software  model  may  limit  the  versatility  of  the  software  model  with  computing  an  off-axis  point 
spread  function  for  the  simulated  sensor  optics  system. 

2.2  Propagation  with  wave  optics 

The  complex  electric  field  for  electromagnetic  radiation  is  used  in  wave  optics  models  to 
simulate  the  propagation  of  electromagnetic  radiation  through  a  medium.  During  propagation 
through  different  optical  components  the  complex  electric  field  across  the  wavefront  of  radiation 
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experiences  effects  optical  path  lengths  due  to  the  type  of  medium  and  length  of  propagation 
through  the  medium.  For  the  sensor  optics  system,  the  lenses  and  direct  vision  prism  have  mediums 
of  propagation  that  differ  from  air  and  will  impart  different  optical  path  lengths  of  propagation  for 
each  region  of  the  electric  field  wavefront. 

2.2.1  Propagation  through  media.  The  optical  path  lengths  are  influenced  by  the  speed  of 
propagation  and  the  duration  of  propagation  through  the  medium.  The  speed  of  propagation,  more 
commonly  described  as  the  speed  of  light,  through  a  medium  depends  on  the  index  of  refraction  for 
the  medium.  Index  of  refraction  for  a  certain  medium  describes  the  ratio  of  the  speed  of  light  in  a 
vacuum  c  =  -p===,  with  respect  to  the  speed  of  light  in  a  certain  medium  v  =  and  is  generally 
expressed  as  n  =  ^  =  \J  The  terms  in  these  expressions  are  index  of  refraction,  n,  dielectric 
constants  of  permeability  for  a  vacuum  e0i  permeability  for  a  certain  medium  e,  permittivity  for  a 
vacuum  poi  and  permittivity  for  a  certain  medium  /i.  The  duration  of  propagation  depends  on  the 
distance  of  propagation  through  the  medium. 

These  considerations  with  wave  optics  thus  captures  the  change  in  direction  of  the  propagation 
due  to  a  change  in  propagation  medium.  The  mathematical  relationship  that  describes  the  change 
in  wavefront  orientation  and  so  the  direction  of  propagation,  is  Snell’s  Law  of  Refraction  stated  in 
equation  2.5, 

ni  sin  6*i  =  n2  sin  O2  (2-5) 

where  the  angle  of  incidence  9\  in  the  first  media  at  the  surface  interface  between  two  different 
media,  and  the  index  of  refraction  of  the  first  medium  ni  where  the  angle  of  incidence  occurs,  is 
related  to  the  angle  of  refraction  @2  from  the  surface  interface  between  the  two  media,  and  the 
index  of  refraction  of  the  second  medium  n2  where  the  angle  of  refraction  occurs. 

The  Lensmaker’s  formula  uses  Snell’s  Law  of  Refraction  for  the  interface  between  propagation 
media  at  each  surface  of  the  lens  to  paraxially  approximate  the  optical  path  length  for  propagation 
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through  the  lens.  This  approximation  also  treats  the  lens  as  a  thin  lens,  and  so  the  Lensmaker’s 
formula  expressed  at  equation  2.6  is  also  referred  to  as  the  Gauss  Thin  Lens  formula  for  lenses, 


^observe 


3  image 


(2.6) 


where  f  denotes  the  focal  length  of  the  lens,  s0bserve  denotes  the  distance  the  subject  is  located  from 
the  lens,  and  Simage  the  distance  the  real  image  is  located  from  the  lens. 


2.2.2  Propagation  of  point  sources.  Interference  patterns  in  the  wavefront  along  the  path 
of  propagation  due  to  the  diffraction  of  light  passing  through  an  aperture,  is  a  phenomena  with  the 
propagation  of  light  that  can  be  emulated  with  simulation,  provided  an  appropriate  approach  is  used 
to  represent  the  propagation  of  light  as  a  complex  electromagnetic  field.  An  approach  that  captures 
the  essence  of  the  diffraction  experienced  in  the  electromagnetic  field  wavefront  propagation  of  light 
is  desirable  for  this  purpose. 

The  Rayleigh-Sommerfield  diffraction  theory  [16:  52]  describes  the  diffraction  effects  resulting 
from  the  propagation  of  a  spherical  electromagnetic  field  wavefront  through  an  aperture.  Correct 
use  of  the  Rayleigh-Sommerfield  diffraction  formula  representing  this  theory  and  described  at  equa¬ 
tion  2.7  in  the  software  model  should  achieve  simulating  the  diffraction  effects  inherent  with  this 
wavefront  propagation. 


Ui(xi,yi)  =  7^  jj  U°(x°’y°)  e  J  k  r  cosq  dxQ  dyo  (2.7) 

This  equation  describes  the  complex  electromagnetic  field  wavefront  Uo(x0,yo)  propagat¬ 
ing  from  the  transverse  plane  with  ordinates  xo  and  yo,  to  another  transverse  plane  with  ordi¬ 
nates  xi  and  yi,  resulting  in  the  complex  electromagnetic  field  wavefront  Ui(xi,yi).  The  dis¬ 
tance  of  propagation  for  a  point  source  on  the  wavefront  of  the  electromagnetic  radiation  is, 
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Figure  2.6:  Huygens’s  Principle  [2:  104] 


r(xo,yo,x1,y1)  =  y 1  (xq  —  x0)2  +  (yi  —  yo)2  +  z2  in  metres,  and  the  wave  number  for  each  wave¬ 
length  in  the  spectrum  of  electromagnetic  radiation  is,  k  =  ^  in  metres-1. 

To  appreciate  how  wave  optics  treats  a  propagating  electromagnetic  field  of  light,  a  sketch 
representing  Huygens’s  Principle  is  included  at  Figure  2.6.  This  figure  shows  the  wavefront  of  the 
electromagnetic  field  of  light  propagating  as  many  (uncountable)  point  sources  distributed  along 
the  wavefront,  each  emitting  spherical  wave  fields  of  electromagnetic  radiation. 

Although  the  propagation  of  a  field  of  electromagnetic  radiation  as  light  has  vectorial  prop¬ 
erties,  the  Rayleigh-Sommerfield  diffraction  theory  applies  scalar  properties  to  describe  the  prop¬ 
agation  of  this  field  in  an  optics  system  based  on  certain  key  simplifications  and  assumptions. 
Electromagnetic  interactions,  described  by  Maxwell’s  equations,  between  the  propagating  field  and 
other  electromagnetic  fields  at  surface  boundaries  and  apertures  of  optical  components  in  the  optics 
system  are  disregarded  in  the  Rayleigh-Sommerfield  diffraction  theory.  Coupling  between  electric 
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and  magnetic  fields  in  the  propagating  electromagnetic  radiation  are  also  disregarded  since  the 
propagation  media  at  all  stages  of  propagation  through  the  system  does  not  support  this  coupling, 
earlier  referred  to  as  electromagnetic  uncertainties  in  section  2.1.2. 

Treating  the  propagation  of  electromagnetic  radiation  as  a  scalar  phenomenon  is  acceptable 
provided  two  conditions  in  the  circumstances  of  propagation  are  satisfied.  These  conditions  are  (1) 
the  diffracting  aperture  must  be  significantly  large  compared  to  the  wavelength  of  electromagnetic 
radiation,  and  (2)  the  electromagnetic  field  resulting  after  the  diffracting  aperture  must  be  observed 
from  a  sufficiently  distant  location  from  the  diffracting  aperture  [16:  35]. 

The  criteria  for  the  terminology  significantly  large  and  a  sufficiently  distant  location  is  not 
specifically  defined.  For  a  typical  aperture  or  lens  diameter  from  2  cm  to  20  cm,  and  typical 
observation  distances  from  10  cm  to  100  cm  respectively,  in  the  infrared  spectrum  in  terms  of 
wavelength  from  2  ^m  to  5  /jin,  the  relationship  between  these  distances  and  wavelength  can  be 
described  in  terms  of  ratios  with  magnitudes  of  at  least  =  4,  000.  This  represents  at  least 

three  orders  of  magnitude  in  size  between  the  spectral  wavelength  and  the  distances  mentioned  in 
the  typical  conditions. 

Another  requirement  of  the  Rayleigh-Sommerfield  diffraction  theory  is  for  each  diffracting 
aperture  to  be  planar  [16:  51].  This  criteria  is  represented  in  the  expression  for  the  propagation 
distance  (denoted,  r)  and  the  term  cos 9  which  describes  the  orientation  of  the  direction  of  propa¬ 
gation  with  respect  to  the  surface  normal  of  the  aperture  plane.  This  condition  can  be  satisfied  at 
each  stage  of  propagation  for  the  software  model,  since  the  field  and  intensity  of  the  electromagnetic 
radiation  are  sampled  across  a  planar  cross-section  of  the  sensor  optics  system. 

2.3  Proposed  system  requirements  and  specifications 

With  developing  a  model  to  simulate  an  actual  system  of  interest  for  physical  design  and 
operation,  the  requirements  of  any  system  should  be  translated  to  the  design  of  the  model.  For  the 
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purpose  of  constructing  this  hyperspectral  imaging  system,  certain  attributes  of  the  system  may 
be  more  restrictive  on  choices  for  the  sensor  design  than  others,  so  it  is  important  to  capture  all 
the  attributes  that  may  influence  how  the  design  is  undertaken. 

The  principle  aim  of  this  software  model  is  to  retain  sufficient  flexibility  with  defining  system 
parameters,  so  that  design  choices  can  be  evaluated  effectively  using  the  model. 

2.3.1  Hardware.  Some  hardware  attributes  of  the  system  that  must  be  accommodated 
with  developing  the  model  of  the  sensor  optics  system  include  the  following  notes. 

2.3.1. 1  Detector  -  focal  plane  array.  The  infrared  detector  is  an  Indium  Anti- 
monide  (InSb)  focal  plane  array,  comprising  512  pixels  by  640  pixels  with  a  0.24  /jm  pitch  per  pixel 
for  height  and  width  equally  across  the  array. 

2. 3. 1.2  Detector  -  focussing  lens.  The  choice  focussing  lens  that  focuses  radiation 
incident  on  the  infrared  detector  to  the  InSb  focal  plane  array  can  be  selected  to  have  an  effective 
focal  length  of  either  50  mm,  100  mm,  150  mm,  200  mm,  250  mm,  and  500  mm. 

2. 3. 1.3  Stepper  Motor.  The  direct  vision  prism  will  be  mounted  in  a  high  pre¬ 
cision  axially  tubular  stepper  motor,  that  can  hold  and  axially  rotate  optical  components  that 
do  not  exceed  2  inches  in  diameter.  The  direct  vision  prism  diameter  must  not  exceed  2  inches 
(approximately  50.8  mm). 

2. 3. 1-4  Direct  Vision  Prism.  The  direct  vision  prism  diameter  is  one  inch.  The 
minimum  axial  distance  along  the  outer  axial  surface  for  each  bulk  media  stage  of  the  prism  is 
3  mm. 

2.3.2  Performance.  Some  hardware  attributes  of  the  system  that  must  be  accommodated 
with  developing  the  model  of  the  sensor  optics  system  include  the  following  notes. 
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2.3.2. 1  Spectral  coverage.  The  spectral  range  of  infrared  electromagnetic  radiation 
to  be  collected  and  imaged  by  the  sensor  optics  system  is  1.8  to  5.2  pm. 

2. 3. 2. 2  Undeviated  spectral  wavelength.  The  undeviated  spectral  wavelength  of 
infrared  electromagnetic  radiation  is  to  be  3.6  ±  0.1/rrn. 

2. 3. 2. 3  Usable  aperture.  The  diameter  of  the  collimated  beam  of  infrared  electro¬ 
magnetic  radiation  entering  the  direct  vision  prism  and  propagating  through  the  prism  is  to  be  at 
minimum  25  percent  of  the  prism  height  without  reflecting  off  the  outer  axial  surface  of  the  prism  . 

2. 3. 2. 4  Spatial  extent  of  images.  The  spatial  extent  of  the  imaging  area  for  each 
spectral  wavelength  band  is  not  to  exceed  128  pixels  by  128  pixels  on  the  focal  plane  array  of  the 
infrared  detector. 


26 


III.  Methodology 


This  chapter  covers  the  design  of  the  software  model  for  simulating  the  sensor  optics  sys¬ 
tem.  The  functional  attributes  of  the  system  are  covered  in  general  at  section  3.1.  A 
generic  treatment  of  the  spatial  sampling  requirements  in  order  to  accurately  sample  the  simulated 
electromagnetic  radiation  wavefront  is  outlined  at  section  3.2.  The  system  coordinates  and  sign 
conventions  defined  for  the  software  model  are  presented  at  section  3.3.  The  impetus  for  pursu¬ 
ing  the  software  model  is  restated  at  section  3.4  before  the  functional  sub-systems  of  the  software 
model  are  discussed  at  section  3.5  (forward  optics  system),  section  3.6  (prism  optics  system)  and 
section  3.7  (detector  optics  system). 

3.1  General  Description 

Designing  the  software  model  of  the  sensor  optics  system  involved  dividing  the  whole  system 
into  subsystems  that  were  logically  separable.  Fortunately  most  optics  systems  can  be  viewed  in 
simple  terms  as  cascading  stages  of  electromagnetic  field  wavefront  propagation,  and  the  sensor 
optics  system  is  an  excellent  example  of  this. 

The  principal  influences  on  the  remote  scene  radiation  performed  by  the  system  are  collec¬ 
tion,  dispersion,  and  recording,  so  a  logical  approach  was  to  divide  the  LSTI  system  described 
at  section  2.1.2  into  three  sub-systems.  Figure  3.1  shows  the  sensor  optics  system  divided  into 
three  separable  subsystems  designated,  forward  optics  system,  prism  optics  system,  and  detector 
optics  system,  according  to  the  influences  of  collecting,  dispersing  and  detecting  the  remote  scene 
radiation  . 

3.1.1  Forward  Optics  System.  From  collecting  all  the  incident  fields  of  electromagnetic 
radiation  from  the  remote  scene  and  surrounding  environment  propagating  through  the  forward 
lens  is  the  incoming  wavefront  of  electromagnetic  radiation  to  the  forward  optics  system.  This 
electromagnetic  field  wavefront  is  focussed  in  the  field  stop  transverse  plane.  The  aperture  size  set 
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Figure  3.1:  Description  of  the  CTIS  Sensor  Optics  System  as  separate  subsystems 


at  the  field  stop  determines  the  angular  extent  of  the  collected  electromagnetic  radiation  propagates 
from  the  field  stop  through  to  the  aft  lens.  With  the  distances  between  each  lens  and  the  field  stop 
set  appropriately  so  a  minified  real  image  of  the  remote  scene  is  formed  in  the  field  stop,  the  field 
of  electromagnetic  radiation  propagating  from  the  aft  lens  should  be  a  collimated  wavefront.  As 
described  here,  the  task  of  collecting  electromagnetic  radiation  from  the  remote  scene  is  achieved 
with  the  forward  optics  system. 


3.1.2  Prism  Optics  System.  The  field  of  electromagnetic  radiation  leaving  the  aft  lens  is 
the  incoming  wavefront  of  electromagnetic  radiation  to  the  prism  optics  system.  The  description  of 
the  wavefront  as  a  field  here  implies  there  are  many  spectral  wavelengths  contained  in  the  wavefront, 
effectively  many  (uncountable)  fields  of  discrete  electromagnetic  radiation.  The  incoming  wavefront 
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propagates  through  the  different  stages  of  media  as  material  for  the  direct  vision  prism.  The 
separation  of  air  between  the  prism  and  the  aft  lens  and  detector  focussing  lens  are  also  considered 
stages  of  propagation  in  the  prism  optics  system.  The  materials  that  comprise  the  prism  are  selected 
as  a  media  of  propagation  to  offer  more  favourable  transmission  properties  and  some  dispersion  for 
the  range  of  spectral  wavelengths  of  interest  for  detection.  The  collimated  wavefront  entering  the 
prism  media  experiences  dispersion  according  to  the  spectral  wavelengths  in  one  direction  along 
the  transverse  plane  of  the  prism  optics  system  in  the  course  of  propagating  through  the  direct 
vision  prism.  The  wavefront  incident  at  the  detector  focussing  lens  is  the  result  of  dispersion  on 
the  wavefront  of  collected  electromagnetic  radiation  from  the  remote  scene,  due  to  propagation 
through  the  prism  optics  system.  This  system  achieves  the  task  of  dispersing  the  collected  remote 
scene  electromagnetic  radiation. 

The  prism  used  in  this  system  is  permitted  to  rotate  around  the  optic  axis  and  disperses  the 
collected  radiation  across  the  transverse  plane  along  an  axis  of  the  prism  cross  section  aligned  with 
the  transverse  plane  to  the  optic  axis.  The  orientation  of  the  direction  of  dispersion  depends  on 
the  orientation  of  the  prism  with  respect  to  the  sensor  optics  system  coordinates.  As  the  prism  is 
rotated  angularly  around  the  optic  axis,  the  real  image  of  the  remote  scene  formed  on  the  detector 
focal  plane  array  as  a  dispersed  image,  also  rotates  around  the  optic  axis.  The  spectral  content  of 
the  collected  radiation  is  dispersed  in  the  transverse  plane  in  both  directions  from  the  optic  axis 
according  to  wavelength,  with  the  exception  of  one  wavelength  band.  This  exclusive  wavelength 
band  of  the  remote  scene  forms  a  real  image  of  the  remote  scene  centred  on  the  optic  axis,  and  is 
referred  to  as  the  undeviated  wavelength.  For  this  wavelength  band,  the  real  image  of  the  remote 
scene  does  not  change  location  on  the  focal  plane  array  as  the  angular  orientation  of  the  prism 
changes. 

3.1.3  Detector  Optics  System.  The  field  of  electromagnetic  radiation  entering  the  detector 
focussing  lens  is  the  incoming  wavefront  of  electromagnetic  radiation  to  the  detector  optics  system. 
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The  dispersion  contained  in  this  wavefront  establishes  different  angular  directions  of  propagation 
with  respect  to  the  optic  axis  for  the  spectral  wavefront  fields  of  electromagnetic  radiation  contained 
in  the  incident  wavefront  of  electromagnetic  radiation  at  the  detector  focussing  lens.  The  lens 
focusses  the  incident  spectral  wavefront  fields  on  the  detector  focal  plane  array. 

The  detector  focussing  lens  interrogates  each  spectral  wavefront  field  of  electromagnetic  ra¬ 
diation  in  a  transverse  plane  to  the  horizontal  plane  of  the  sensor  optics  system  that  contains  the 
optic  axis.  The  view  of  the  wavefront  obtained  by  the  detector  focussing  lens  differs  according  to 
the  angular  direction  of  propagation  for  each  spectral  wavelength  of  electromagnetic  fields  in  the 
wavefront.  In  essence  this  view  conveys  a  certain  gradient  or  tilt  in  phase  (phase  tilt)  across  the 
wavefront  in  the  direction  of  dispersion  from  the  direct  vision  prism.  The  amount  of  phase  tilt  for 
the  electromagnetic  field  of  each  spectral  wavelength  can  be  gauged  by  comparing  these  wavefronts 
incident  on  the  detector  focussing  lens  with  equivalent  wavefronts  propagating  with  no  angular 
direction  of  propagation  with  respect  to  the  optic  axis  (ie.  along  the  optic  axis).  The  location  in 
the  transverse  plane  of  the  detector  focal  plane  array  where  the  real  image  of  the  remote  scene  for 
a  spectral  wavelength  of  electromagnetic  radiation  is  influenced  by  the  phase  tilt  on  the  spectral 
wavefront  field. 

The  detector  optics  system  achieves  the  task  of  recording  the  collected  remote  scene  electro¬ 
magnetic  radiation  in  a  format  that  separates  the  spectral  wavelength  content  of  the  remote  scene 
to  enable  retrieval  of  real  images  of  the  remote  scene  for  each  spectral  wavelength  of  electromag¬ 
netic  radiation  recorded.  For  an  upright  oriented  direct  vision  prism,  the  angular  directions  of 
propagation  for  every  spectral  wavefront  field  of  electromagnetic  radiation  in  the  wavefront  change 
along  the  vertical  axis  of  the  system,  applying  vertical  phase  tilts  to  the  spectral  wavefront  fields. 
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3.2  Sampling  Requirements 


The  purpose  of  the  software  model  is  to  represent  a  discrete  version  of  the  collected  fields 
of  electromagnetic  radiation  from  the  remote  scene  that  propagate  through  the  sensor  optics  sys¬ 
tem.  To  ensure  these  fields  are  accurately  represented  as  discretely  sampled  spectral  fields,  the 
requirements  of  spatial  sampling  across  the  apertures  contained  in  the  sensor  optics  system  must 
be  assessed. 

3.2.1  General  Description.  Imaging  systems  typically  consist  of  lenses  and  apertures,  and 
the  optical  performance  of  these  systems  is  limited  by  the  propagation  of  electromagnetic  fields 
through  the  imaging  system.  These  fields  carrying  the  remote  scene  spatial  feature  information, 
and  the  size  of  apertures  and  lenses  influences  the  optical  performance  of  imaging  systems  with 
reproducing  the  real  image  of  the  remote  scene. 

Using  Fourier  optics  analysis  between  two  transverse  planes,  the  spatial  frequency  response 
of  an  imaging  system  can  be  evaluated  with  some  important  fundamentals  for  diffraction-limited 
imaging  [16:  134-144].  An  Amplitude  Transfer  Function  (ATF)  describes  the  spatial  frequency 
response  for  a  coherent  imaging  system,  whereas  an  Optical  Transfer  Function  (OTF)  describes 
the  spatial  frequency  response  for  an  incoherent  imaging  systems.  So  for  one  physical  configuration 
of  optics,  the  optical  performance  related  to  the  diffraction-limited  spatial  frequency  cut-off,  can 
be  assessed  from  the  ATF,  whether  the  optics  system  is  used  in  a  coherent  or  incoherent  imaging 
system. 

The  diffraction-limited  spatial  frequency  cut-off  for  a  field  of  electromagnetic  radiation  with 
spectral  wavelength,  A,  propagating  through  an  aperture  with  diameter,  dia,  in  a  direction  adjacent 
to  the  optic  axis  to  some  location  in  a  transverse  plane  to  the  optic  axis,  with  the  distance,  z, 
separating  the  aperture  and  the  transverse  plane  measured  along  the  optic  axis,  is  described  at 
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equation  3.1  for  coherent  imaging,  and  at  equation  3.2  for  incoherent  imaging, 


dia 

fco(f)  ~~  2Az 


(3.1) 


dia 

fco(i)  “  ~Xz 


(3.2) 


where  fco(f)  describes  the  cut-off  spatial  frequency  for  the  field  amplitude  related  to  the  ATF  of  the 
optics  system,  and  fco(n  describes  the  cut-off  spatial  frequency  for  the  field  intensity  related  to  the 
OTF  of  the  optics  system.  Comparing  these  cut-off  spatial  frequencies,  fco(j)  =  2  fco(f). 


3.2.2  Software  Model.  The  sensor  optics  system  of  the  CTIS  is  used  for  incoherent 
imaging,  so  the  sampling  requirements  for  the  software  model  to  represent  the  incident  fields  of 
electromagnetic  radiation  as  discretely  sampled  fields  were  treated  for  incoherent  imaging.  The 
principles  of  geometric  optics  will  be  used  to  define  the  relationships  for  obtaining  the  discretely 
sampled  fields. 

Viewing  the  arrangement  of  the  software  model  for  the  sensor  optics  system  between  the 
remote  scene  and  the  field  stop  with  the  forward  lens  situated  between  these  transverse  planes  as 
shown  in  Figure  3.2,  the  geometric  optics  principle  that  rays  propagating  through  the  centre  of  the 
lens  pass  through  undeviated,  establishes  the  following  relationship  at  equation  3.3  between  the 
sampling  intervals  at  the  field  stop  transverse  plane  ds3,  and  the  remote  scene  transverse  plane  dsi. 


dsi  ds3 
zi  z2 


(3.3) 


For  the  aperture  of  the  forward  lens  with  a  diameter  dia2,  the  cut-off  spatial  frequency  at  the 
field  stop  is  expressed  at  equation  3.4  for  a  diffraction-limited  system  based  on  coherent  imaging. 


_  dia2  1 

W)  - 


(3.4) 
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Figure  3.2:  Interpretation  of  spatial  sampling  intervals  for  the  forward  lens  with  geometric  optics 

Modifying  this  cut-off  spatial  frequency  for  an  optics  system  used  for  incoherent  imaging,  the 
cut-off  spatial  frequency  for  forming  a  real  image  at  the  field  stop  described  at  equation  3.5. 


f  _of  _  dia2 
l3co(i)  —  ^3co(f)  —  \ 

AZ2 


(3.5) 


The  spatial  rate  of  sampling,  termed  the  spatial  sampling  frequency,  must  record  this  cut¬ 
off  spatial  frequency,  so  obeying  the  Nyquist  criteria  for  spatial  sampling  so  the  complete  extent 
of  phase  variations  in  the  fields  of  electromagnetic  radiation  are  accurately  recorded,  the  spatial 
sampling  frequency  as  a  minimum  must  be  twice  the  cut-off  spatial  frequency  as  described  in 
equation  3.6. 

f3s  =  2f2co(i)  =  ^  (3.6) 

AZ2 
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The  spatial  sampling  interval  in  the  transverse  plane  of  the  field  stop  is  the  inverse  of  this 


spatial  sampling  frequency  and  is  described  in  equation  3.7. 

_  _  Az2 

S3  f3s  2dia2 


(3.7) 


Revisiting  the  relationship  between  the  spatial  sampling  intervals  dsi  and  ds3  at  equation  3.3, 
the  spatial  sampling  interval  in  the  transverse  plane  of  the  remote  scene  can  be  described  in  equa¬ 


tion  3.8. 


dsi 


=$■  dsi 


Azi 

2dia2 


(3.8) 


An  aperture  located  in  the  remote  scene  can  be  notionally  formed  around  the  extent  of  the 
remote  scene  observed  through  the  sensor  optics  system,  for  the  purpose  of  evaluating  the  number 
of  samples  required  to  discretely  represent  the  remote  scene,  and  described  to  have  a  diameter, 
dia3  =  Ndsi.  The  cut-off  spatial  frequency  at  the  forward  lens  f2co(f)  for  the  notional  aperture  in 
the  remote  scene  is  described  in  equation  3.9,  for  a  diffraction-limited  system  based  on  coherent 
imaging.  The  number  of  samples  N,  represents  the  total  number  of  samples  across  the  extent  of 
the  remote  scene. 


^2co(f)  — 


=>  f2co(f)  = 


(3.9) 


This  expression  represents  the  field  amplitude  cut-off  spatial  frequency  at  the  forward  lens. 
For  sampling  the  fields  of  electromagnetic  radiation  at  the  forward  lens,  this  cut-off  spatial  frequency 
can  be  retained  because  there  is  no  real  image  formed  at  the  forward  lens.  Obeying  the  Nyquist 
sampling  criteria,  the  spatial  sampling  frequency  as  a  minimum  is  set  to  twice  this  cut-off  spatial 
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frequency,  as  described  in  equation  3.10. 


f2s  =  2f. 


2co(f) 


Ndsi 

Azi 


(3.10) 


The  sampling  spatial  interval  at  the  forward  lens  is  the  inverse  of  this  spatial  sampling  fre¬ 
quency,  as  described  in  equation  3.11. 


ds2  =  —  = 
l2s 


Azi 

Ndsi 


(3.11) 


Noting  the  spatial  sampling  interval  at  the  remote  scene  dsi  from  equation  3.3  and  the  number 
of  spatial  samples  across  the  spatial  extent  of  the  remote  scene  N,  the  size  of  the  remote  scene  can 
be  described  as  Ndsi,  and  used  to  evaluate  the  spatial  sampling  interval  at  the  forward  lens  ds2. 
The  expression  used  to  evaluate  ds2  appears  at  equation  3.12. 


ds2 


Azi 

Ndsi 


=>  ds2 


/  Azi  A  (  2dia2  A  2dia2 

(it)  (tst)  dS2  =  — 


(3.12) 


The  diameter  of  the  forward  lens  reveals  the  total  number  of  spatial  samples  M  across  the 
diameter  of  the  forward  lens  required  to  sample  the  incident  fields  of  electromagnetic  radiation 
described  at  equation  3.13  represents  half  the  number  of  spatial  samples  required  to  discretely 
sample  the  remote  scene. 


N 

2 


(3.13) 


With  the  real  image  of  the  remote  scene  formed  in  the  field  stop,  the  ray  propagation  principles 
of  geometric  optics  were  used  to  establish  the  spatial  sampling  intervals  for  the  aft  lens  transverse 
plane.  Propagating  rays  of  electromagnetic  radiation  incident  on  the  field  stop  converge  at  the 
field  stop,  then  diverge  past  the  field  stop  to  propagate  toward  the  aft  lens  as  shown  in  Figure  3.3. 
The  geometry  of  this  propagation  for  the  respective  spatial  sampling  intervals  at  the  forward  lens 
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aperture  ds2  and  aft  lens  aperture  ds4  is  expressed  in  equation  3.14. 


dS2  dS4 
Z2  z3 


(3.14) 
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Figure  3.3:  Interpretation  of  spatial  sampling  intervals  for  the  field  stop  with  geometric  optics 

Using  the  forward  lens  expression  from  equation  3.12  and  sampling  the  spectral  wavefront 
fields  of  electromagnetic  radiation  at  the  aft  lens  for  the  same  reasons  described  previously  for  the 
forward  lens,  the  spatial  sampling  interval  in  the  transverse  plane  of  the  aft  lens  ds4  was  evaluated 
using  equation  3.15. 

ds<  =  d82(S)  ^  dS4  =  (^r)  (!)  (3d5) 

The  spatial  sampling  intervals  have  a  dependence  on  the  spectral  wavelength  of  every  spectral 
field  of  electromagnetic  radiation  being  discretely  sampled.  To  sample  these  fields  sufficiently  for  all 
spectral  wavelengths  the  spatial  sampling  interval  required  for  the  minimum  spectral  wavelength 
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of  electromagnetic  radiation  is  used  in  the  software  model  for  all  spectral  wavelengths,  ensuring  all 


spectral  wavefront  fields  are  sufficiently  sampled. 


3.3  System  coordinates  and  sign  conventions 

A  system  defining  the  orientations  of  position  in  an  arbitrary  transverse  vertical  plane,  and 
angles  of  trajectory  with  respect  to  the  optic  axis  lying  in  an  axial  horizontal  plane,  was  required 
for  developing  the  software  model.  The  reference  system  used  throughout  the  software  model  and 
shown  at  Figure  3.4  defines  the  positions  and  angles  of  propagation  for  the  simulate  electromagnetic 
field,  by  orientation  and  sign  convention. 


y 
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Figure  3.4:  Sensor  optics  system  sign  conventions  and  coordinate  system  defined  for  the  designing 

the  software  model 


The  system  coordinates  chosen  for  the  software  model  align  the  x-axis  horizontal,  the  y-axis 
vertical,  so  both  axes  lie  in  a  transverse  plane  to  the  horizontal  optic  axis  oriented  as  the  z-axis. 
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The  sign  convention  for  locations  of  position  follows  the  right-hand  rule  used  in  electromagnetic 
theory,  where  the  direction  of  propagation  is  defined  as  the  positive  direction  along  the  z-axis. 

The  sign  convention  for  angles  of  propagation  follows  typical  geometric  optics  sign  convention. 
A  plane  of  propagation  oriented  vertically  (in  yz-plane),  where  the  angular  direction  of  dispersion 
varies  vertically  according  to  the  relationships  between  angles  of  incidence  and  angles  of  refraction, 
was  used  to  establish  the  sign  convention  use  in  the  software  model.  Positive  angles  represent  angles 
subtended  to  the  propagation  ray  of  electromagnetic  radiation  from  the  optic  axis  in  the  upward 
direction,  and  negative  angles  represent  angles  subtended  to  the  propagation  ray  from  the  optic 
axis  in  the  downward  direction. 

3-4  Software  Model  Development 

Development  of  the  sensor  optics  system  software  model  closely  followed  the  subsystems 
identified  previously  in  section  3.1.  The  mathematics  and  graphics  laboratory  tool  MatLab®  was 
chosen  as  the  facilitating  software  to  perform  the  instructions  scripted  for  simulating  the  system. 
Figure  3.5  describes  the  simplified  block  diagram  of  the  software  model  and  includes  the  forward 
optics  system,  prism  optics  system,  and  detector  optics  system  as  subsystems  of  the  entire  system. 

The  significance  of  defining  the  point  spread  functions  for  the  sensor  optics  system  to  aid  with 
characterising  the  system  and  aid  with  reconstruction  of  the  remote  scene  as  a  hyperspectral  real 
image  set  were  discussed  at  section  2.1.  The  software  model  simulates  the  point  spread  functions 
for  the  system.  The  source  of  electromagnetic  radiation  used  for  the  software  model  simulates  fields 
of  electromagnetic  radiation  across  the  range  of  spectral  wavelengths  from  2  /mi  to  5  n m  emitted 
from  the  remote  scene  as  spherical  wavefronts  of  electromagnetic  radiation  from  a  non-discernible 
point  source. 
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Figure  3.5:  CTIS  Sensor  Optics  System  description  for  software  model 
3.5  Forward  Optics  System 

The  design  of  the  forward  optics  system  as  a  software  model  is  represented  as  a  simple  block 
diagram  at  Figure  3.6. 

The  point  source  was  simulated  as  a  non-discernible  source  emitting  spherical  wavefronts  of 
electromagnetic  radiation  in  the  discrete  spectral  wavelength  region  between  2  /jin  and  5  /iin.  Prop¬ 
agation  of  this  radiation  as  spherical  electromagnetic  wavefronts  is  simulated  using  the  Rayleigh- 
Sommerfield  diffraction  formula  through  the  various  stages  of  propagation  in  the  forward  optics 
system.  The  discrete  version  of  the  Rayleigh-Sommerheld  diffraction  formula  used  to  compute  dis¬ 
crete  samples  of  the  electromagnetic  field  for  each  spectral  wavelength  at  the  destination  aperture 
of  propagation,  is  represented  at  equation  3.16  where  Nx  and  Ny  represent  the  number  of  discrete 
samples  dxo  and  dyo  across  the  forward  lens  aperture  of  the  electromagnetic  field  wavefront  in 
directions  along  the  x-axis  and  y-axis  respectively.  Nx  and  Ny  were  chosen  to  be  odd  integers  for 
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Figure  3.6:  Forward  Optics  System  description  for  software  model 


symmetry  around  the  optic  axis. 


Ui(xi,yi,  A) 


jA 


i(Nx-l)  |(Ny-l) 

E  E 

x=-I(Nx-l)ny=-i(Ny  — 1) 


U0(nxdx0,  nvdy0,  A) 


J  w 


cos 6  nxdxp  nydyo  (3.16) 


The  distance  of  propagation  r,  to  the  destination  transverse  plane  from  the  emitting  transverse 


plane  is  described  at  equation  3.17. 


r  =  r(x0,y0,x1,y1)  =  r(nxx0,  nyy0,  xl5  y^  =  sj (x:  -  x0)2  +  (yi  -  y0)2  +  z2  (3.17) 


Simulating  the  electromagnetic  fields  discretely  with  the  Rayleigh-Sommerfield  diffraction  for¬ 
mula  in  this  way  provides  a  wave  optics  treatment  of  the  propagation  of  electromagnetic  radiation 
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through  the  forward  optics  system.  Once  the  sampling  requirements  are  established  for  each  trans¬ 
verse  plane,  and  the  image  distance  for  the  spectral  wavelengths  evaluated  for  locating  the  field  stop 
along  the  optic  axis,  the  wave  optics  propagation  of  the  electromagnetic  radiation  was  simulated 
for  the  stages  of  propagation:  (1)  between  the  point  source  and  the  forward  lens  transverse  planes; 
(2)  between  the  forward  lens  and  the  field  stop  transverse  planes;  and  (3)  between  the  field  stop 
and  the  aft  lens  transverse  planes,  using  the  Rayleigh-Sommerfield  diffraction  formula. 

Simulating  the  wave  optics  propagation  through  the  forward  and  aft  lenses  respectively  in¬ 
volved  computing  additional  phase  to  each  sample  of  the  electromagnetic  field.  Using  the  lens 
thickness  formula  for  each  lens,  the  optical  path  length  for  each  sample  of  the  electromagnetic  field 
was  evaluate  and  used  to  compute  the  additional  phase.  The  lens  thickness  formula,  including 
the  paraxial  approximation  used  for  representing  a  spherical  lens  surface  as  a  parabolic  surface,  is 
described  at  equation  3.18. 


Ai(x0,y0) 


(3.18) 


Using  this  lens  thickness  expression,  the  additional  phase  for  all  spatial  samples  across  ev¬ 
ery  spectral  wavefront  field  of  electromagnetic  radiation  incident  on  the  lens  Uo(xo,yo,A)  after 
propagation  through  the  lens  Ui(xo,yo,A)  was  evaluated  using  equation  3.19. 

Ui(x0,y0,A)  =  U0(x0,y0,  A)e  J  ^  Al(x°’yo)  (niens_1)  (3.19) 


3.6  Prism  Optics  System 

The  general  overview  for  designing  the  prism  optics  system  as  a  software  model  is  represented 
by  the  simple  block  diagram  at  Figure  3.7. 

The  fields  of  electromagnetic  radiation  entering  the  prism  optics  system  are  simulated  as  a 
combined  collimated  wavefront  of  electromagnetic  radiation.  The  wavefront  passes  through  each 
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Figure  3.7:  Prism  Optics  System  description  for  software  model 

stage  of  prism  bulk  media,  experiencing  the  effects  of  refraction  in  the  course  of  propagating  through 
the  prism  optics  system.  The  direction  of  propagation  for  each  spectral  field  of  electromagnetic 
radiation  is  altered  at  each  interface  between  different  prism  bulk  media  according  to  Snell’s  Law 
of  Refraction  described  at  equation  2.5.  Altering  the  direction  of  propagation  for  the  collimated 
wavefront  is  simulated  by  varying  the  angle  of  propagation  with  respect  to  the  optic  axis  along  the 
direction  of  dispersion.  Each  spectral  field  of  electromagnetic  radiation  in  the  collimated  wavefront 
experiences  a  different  change  in  direction  of  propagation  at  each  interface  between  prism  bulk 
media. 
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The  prism  optics  system  in  the  software  model  consists  of  five  separate  bulk  media  propagation 
elements,  meaning  there  are  four  separate  interfaces  in  the  software  model.  The  air  separation 
between  the  aft  lens  and  direct  vision  prism,  also  the  detector  focussing  lens  and  direct  vision  prism, 
are  also  treated  as  propagation  bulk  media  in  the  simulated  prism  optics  system.  The  optical  path 
lengths  for  each  sample  across  each  spectral  held  of  electromagnetic  radiation  are  computed  so  the 
properties  of  wave  optics  are  represented  in  the  simulated  spectral  fields  computed  in  the  software 
model  for  the  prism  optics  system. 

Dispersion  of  the  wavelengths  (spreading  in  location  of  separate  wavelengths)  contained  in 
the  broadband  electric  held  is  achieved  from  two  effects  resulting  from  the  propagation  through  the 
prism.  The  refraction  that  occurs  at  the  surface  interface  between  two  different  propagation  media, 
resulting  from  the  contrast  in  index  of  refraction  between  the  two  propagation  media,  is  one  effect. 
The  variation  of  index  of  refraction  for  the  propagation  media  with  respect  to  wavelength,  for  the 
broadband  electric  held  is  the  other  effect. 

The  angles  measured  with  respect  to  the  surface  normal  for  each  surface  interface  between 
separate  propagation  bulk  media  in  the  prism  are  designated  (3.  Figure  3.8  and  Figure  3.9  show 
the  arrangement  of  these  angles,  with  the  electromagnetic  held  propagation  represented  by  the  ray. 
The  set  of  equations  at  equation  3.20  describe  these  angles. 
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Figure  3.8: 


Aperture  Aperture 


(aft  lens)  Direct  Vision  Prism  (focussing  lens) 


Angles  of  refraction  through  the  forward  section  of  the  direct  vision  prism 


Pi  a  —  OL  l 

Pib  =  ~  arcsin  (^)  sin  /3ia 
/?2a  =  ~Plb 

fob  =  -  arcsin  (j^)  sin  /32a 
/33a  =  ~  (ctl  +  /32b) 

/33b  =  -  arcsin  (^)  sin  /33a 
/34a  =  a3  —  /33b 
/34b  =  -  arcsin  (^)  sin  /34a 

(3.20) 
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Aperture 
(aft  lens) 


Direct  Vision  Prism 


Aperture 
(focussing  lens) 


Figure  3.9:  Angles  of  refraction  through  the  aft  section  of  the  direct  vision  prism 

The  angles  measured  with  respect  to  the  optics  axis  are  designated  as  e.  The  set  of  equations 
at  equation  3.21  describe  these  angles. 
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Cla  —  0 


eib  =  ai  +  (3  ib 

C2a  =  — eib 
C2b  =  ai+  02b 
C3a  =  ^2b 

e3b  =  /?3b 
C4a  =  —  C3b 
^4b  =  0:3+  /?4b 


(3.21) 


The  speed  of  propagation  through  bulk  media,  v,  for  an  electromagnetic  field  is  related  to 
the  index  of  refraction  for  the  media,  n,  and  the  speed  of  propagation  in  free  space,  c,  expressed  in 
equation  3.22.  Knowing  the  distance  of  propagation,  d,  and  the  speed  of  propagation,  the  duration 
of  propagation,  t,  or  time  of  flight  for  the  electromagnetic  radiation  can  be  determined. 


v  = 


c 

n 


d 

t 


(3.22) 


The  time  of  flight  for  each  sample  of  electromagnetic  field  across  the  collimated  wavefront 
is  used  to  evaluate  the  optical  path  lengths  for  each  sample  propagating  through  the  prism  optics 
system.  The  time  of  flight  for  the  entire  prism  optics  system  can  be  represented  by  taking  the  sum 
of  all  times  of  flight  through  the  separate  bulk  media  propagation  elements  for  each  sample  of  the 
collimated  wavefront,  as  described  in  equation  3.23. 

to  =  tl  +  t2  +  t3  +  t4  +  t5  (3.23) 
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Substituting  the  time  of  flight  expression  from  equation  3.22,  the  time  of  flight  for  samples 


across  each  spectral  field  of  electromagnetic  radiation  can  be  alternately  expressed  by  equation  3.24. 


nidi  n2d2  n3d3  n4d4  n5d5 

to  = - 1 - 1 - 1 - 1 - 

c  c  c  c  c 


(3.24) 


The  optical  path  length  (OPL)  for  each  sample  across  the  spectral  field  propagating  through 
the  separate  bulk  media  propagation  elements  is  related  by  the  product  between  speed  of  light  in 
free  space  and  the  time  of  flight  to  propagate  the  length  of  media  in  free  space,  as  described  in 
equation  3.25, 

ct0  =  n4di  +  n2d2  +  n3d3  +  n4d4  +  n5d.5  (3.25) 

where  the  expression,  nidi,  used  separately  for  each  spectral  field,  represents  the  OPLs  of  each 
sample  of  the  field  for  the  first  bulk  media  propagation  element  in  the  prism  optics  system. 

The  distances  of  propagation  for  each  spectral  wavelength  field  of  electromagnetic  radiation 
simulated  through  the  prism  optics  system  were  evaluated  according  to  the  system  coordinates  and 
sign  conventions  stated  in  section  3.3.  The  distances  of  propagation  for  each  stage  of  propagation 
through  the  prism  optics  system  is  described  from  section  3.6.1  to  section  3.6.5.  Figure  3.10  shows 
the  propagation  for  the  first,  second  and  third  bulk  media  stages  of  propagation,  and  Figure  3.11 
shows  the  propagation  for  the  fourth  and  fifth  stages  of  propagation  through  the  prism  optics 
system. 

3.6.1  First  bulk  media  propagation.  Propagation  from  the  aft  aperture  of  the  forward 
optics  system  to  the  first  surface  interface  of  the  prism.  Treat  the  propagation  of  the  wavefront  for 
each  wavelength  as  a  collimated  beam  propagating  along  the  optic  axis  of  the  system. 
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Aperture  Aperture 

(aft  lens)  Direct  Vision  Prism  (focussing  lens) 


Figure  3.10: 


Paths  of  propagation  through  the  forward  section  of  the  prism  optics  system 


'h (offset.)  yi(totai)  tan  Ofl 
zl(total)  “  zl(ref)  zl(offset) 
l'l  Z1  (total) 

(3.26) 


The  OPLs  for  this  propagation  element  of  the  prism  optics  system  are  described  by  the 
expression,  niri. 
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3.6.2  Second  bulk  media  propagation.  Propagation  through  the  first  media  of  the  prism, 
which  is  a  thin  uniform  coating  of  sapphire  of  thickness  dsurface  on  the  surface  of  the  bevelled  forward 
face  of  the  prism.  Treat  the  propagation  of  the  wavefront  for  each  wavelength  as  a  collimated  beam 
propagating  at  an  angle  to  the  optic  axis  of  the  system. 


f2 


^surface 
COS  /3ib 


y2(total)  —  r2  sin(a!  +  /3ib)  +  yi(total) 


(3.27) 


The  OPLs  for  this  propagation  element  of  the  prism  optics  system  are  described  by  the 
expression,  112^. 

3.6.3  Third  bulk  media  propagation.  Propagation  through  the  second  media  of  the  prism, 
which  is  the  first  actual  bulk  media  stage  of  the  prism.  The  forward  surface  for  this  stage  of 
propagation  is  set  at  an  angle  oq  with  respect  to  the  optic  axis  of  the  system,  therefore  giving  the 
appearance  of  a  bevelled  prism  face.  Treat  the  propagation  of  the  wavefront  for  each  wavelength 
as  a  collimated  beam  propagating  at  an  angle  to  the  optic  axis  of  the  system. 

The  equation  list  described  at  equation  3.28  reveals  the  detail  of  computation  used  to  obtain 
the  length  of  propagation  r3  and  the  cumulative  vertical  displacement  of  the  simulated  electromag¬ 
netic  radiation  wavefront  y3(totai)  after  completing  this  stage  of  propagation. 
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^2(ref)  =  2  dictprism  tan  O' i  Z2(min) 

^2(ofFset)  =  Y2(total)  tail  O' i 
^2(total)  ^2(ref)  “1“  ^(offset) 


r  _  z2(total) 

3  —  cos(c*i+/32b) 

y 3(total)  =  r3  sin(ai  +  /?2b)  +  Y2(total) 


(3.28) 


The  OPLs  for  this  propagation  element  of  the  prism  optics  system  are  described  by  the 
expression,  113^. 

3.6.4  Fourth  bulk  media  propagation.  Propagation  through  the  third  media  of  the  prism, 
which  is  the  second  actual  bulk  media  stage  of  the  prism.  The  aft  surface  for  this  stage  of  propaga¬ 
tion  is  set  at  an  angle  <23  with  respect  to  the  optic  axis  of  the  system,  also  giving  the  appearance  of 
a  bevelled  prism  face.  Treat  the  propagation  of  the  wavefront  for  each  wavelength  as  a  collimated 
beam  propagating  at  an  angle  to  the  optic  axis  of  the  system. 


_  _ £3 _  _  Z3(ref)  +  z3(offset)  /g  29) 

4  cos(a3  +  /?3b)  cos(a3  +  /33b) 

An  expression  for  computing  z3(0ffset)  directly  could  not  be  determined,  so  an  alternate  ap¬ 
proach  using  the  geometry  of  the  direct  vision  prism  was  developed  to  compute  the  length  of 
propagation  14.  The  equation  list  described  at  equation  3.30  reveals  the  detail  of  computation 
required  to  obtain  and  the  cumulative  vertical  displacement  of  the  simulated  electromagnetic 
radiation  wavefront  y4(t0tai)  after  completing  this  stage  of  propagation. 
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^3(ref)  —  2  ^iclprism  tail  CM 3  +  Z3(min) 

_  cos  p3  b 
r4(ref)  -  Z3(ref) 

Y4(ref)  ^4(ref)  sin  /^3b  H-  y3(total) 

_  Y4(ref)  sin  <23 

r4(offSet)  -  sin(f_a3+/33b) 
r4  l'l(ref)  r4(offset) 


Y4(total)  —  r4  Sin  /?3b  “1“  Y3(total) 


(3.30) 

The  OPLs  for  this  propagation  element  of  the  prism  optics  system  are  described  by  the 
expression,  114^. 

3.6.5  Fifth  bulk  media  propagation.  Propagation  from  the  aft  surface  of  the  prism  to  the 
focussing  lens  of  the  Infrared  detector.  Treat  the  propagation  of  the  wavefront  for  each  wavelength 
as  a  collimated  beam  propagating  at  an  angle  to  the  optic  axis  of  the  system. 

The  equation  list  described  at  equation  3.31  reveals  the  detail  of  computation  used  to  obtain 
the  length  of  propagation  r5  and  the  cumulative  vertical  displacement  of  the  simulated  electromag¬ 
netic  radiation  wavefront  y.5(totai)  after  completing  this  stage  of  propagation. 
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^5(ofFset)  y4(total)  tail  Q. 3 


^5(total)  —  ^5(ref)  “t"  ^(offset) 

„  _  z5(total) 

5  —  cos(a3+/34b) 

y5(total)  =  r5  sin(a3  +  /34b)  +  y3(total) 


(3.31) 


The  OPLs  for  this  propagation  element  of  the  prism  optics  system  are  described  by  the 
expression,  nsrs. 

Need  to  evaluate  the  phase  across  the  transverse  plane  to  the  optic  axis,  and  modify  the 
exit  angle  of  the  propagation  so  that  this  phase  is  constant  across  the  wavefront  lying  in  the 
transverse  plane.  To  modify  the  exit  angle  of  the  propagation,  the  angles  between  optic  axis  and 
the  respective  entry  and  exit  faces  of  the  prism  along  the  direction  of  dispersion  can  be  modified 
to  tune  the  propagation  through  the  prism  for  the  centre  wavelength. 

The  physical  appearance  of  wedge  shaped  prism  optical  components  was  used  for  the  design  of 
the  software  model  to  represent  the  construct  of  the  direct  vision  prism  bulk  media.  To  simulate  the 
electromagnetic  field  that  propagates  through  each  different  bulk  media,  the  software  model  utilises 
the  design  expectations  that  the  wedge  shaped  bulk  media  has  a  certain  diameter,  and  is  located 
in  the  sensor  optics  system  so  that  the  centre  axis  of  the  bulk  media  is  located  with  the  optics 
axis.  The  wedge  shaped  optical  components  are  specified  typically  by  the  angle  that  the  face  of  the 
wedge  subtends  with  the  transverse  plane  to  the  centre  axis,  and  the  minimum  length  dimension  of 
bulk  media  parallel  to  the  optic  axis,  described  as  zmjn  in  the  software  model.  The  parameters  are 
physical  constraints  of  geometry  with  the  bulk  media  that  must  be  considered  when  evaluating  the 
propagation  of  the  electromagnetic  radiation.  These  parameters  are  used  in  the  software  model  to 
evaluate  the  actual  path  lengths  of  propagation  for  all  samples  of  each  spectral  field.  The  actual 
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lengths  of  propagation  are  used  to  compute  the  optical  path  length  for  the  respective  samples  of 
each  spectral  field. 

3. 6. 6  Prism  rotation.  Rotating  the  direct  vision  prism  should  result  in  circular  translation 
of  the  spectral  real  image  of  the  remote  scene,  based  on  the  spectral  wavelength  electromagnetic 
radiation  specific  to  the  spectral  real  image. 

Sampling  of  the  spectral  field  wavefronts  of  electromagnetic  radiation  from  the  aft  lens  is 
used  as  the  input  spectral  field  wavefronts  for  the  prism  optics  system.  As  the  angular  position 
of  the  direct  vision  prism  changes  due  to  rotation,  these  wavefronts  incident  on  the  prism  entry 
surface  aperture  remain  unchanged.  The  characteristics  of  these  wavefronts  are  modified  by  the 
different  orientations  of  the  prism  entry  surface  aperture  due  to  Snell’s  Law  of  Refraction,  for  each 
angular  position  of  rotation.  Material  properties  of  the  index  of  refraction  and  transmission  for 
bulk  media  used  in  the  direct  vision  prism,  the  incident  angles  of  the  spectral  field  wavefronts  with 
each  respective  prism  bulk  media  interface  are  treated  in  the  wave  optics  software  model  to  model 
the  rotation  of  the  direct  vision  prism. 

This  is  done  by  using  a  separate  coordinate  system,  referred  to  as  prism  coordinates,  for 
the  propagation  of  the  spectral  field  wavefronts  of  electromagnetic  radiation.  With  the  direct 
vision  prism  oriented  in  the  upright  angular  position,  the  prism  coordinates  align  with  the  system 
coordinates. 

3.7  Detector  Optics  System 

The  design  of  the  detector  optics  system  as  a  software  model  is  represented  as  a  simple  block 
diagram  at  Figure  3.12. 

The  detector  optics  system  includes  a  focussing  lens  and  a  focal  plane  array.  The  simulated 
field  of  electromagnetic  radiation  developed  by  the  software  model  only  represents  the  field  incident 
on  the  transverse  plane  where  the  focal  plane  array  is  situated.  The  focal  plane  array  used  for  the 
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Figure  3.12:  Detector  Optics  System  description  for  software  model 


physical  system  must  have  responsivity  properties  suited  to  detecting  photons  at  infrared  spectral 
wavelengths  of  electromagnetic  radiation. 

The  software  model  is  designed  to  simulate  the  point  spread  function  for  each  spectral  wave¬ 
length  band.  Thus  modelling  the  detector  optics  system  must  include  a  process  for  locating  where 
the  point  spread  functions  must  be  in  the  focal  plane  array  transverse  plane.  Without  an  imple¬ 
mentation  in  the  software  model  that  performs  this  task,  assembling  all  spectral  wavelength  point 
spread  functions  as  a  chromotomographic  real  image  would  be  very  difficult  to  achieve. 

The  task  of  locating  the  point  spread  function  is  further  challenged  by  the  sampling  phe¬ 
nomenon  identified  from  the  Whittaker-Shannon  sampling  theorem  for  representing  continuous 
aperiodic  signals  as  discrete  signals  [16:  22-26].  Representing  the  continuous  fields  of  the  point 
spread  functions  as  discretely  sampled  simulated  fields,  the  software  model  duplicates  the  point 
spread  functions  of  the  discrete  fields  in  the  transverse  plane  of  the  focal  plane  array. 

The  software  model  of  the  detector  optics  system  is  required  to  locate  the  discrete  point 
spread  function  for  each  spectral  wavelength  that  represents  the  continuous  field  point  spread 
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function  from  these  duplicate  discrete  point  spread  functions.  Using  the  principles  of  Fourier  optics 
that  relate  the  fields  in  the  transverse  planes  of  the  detector  focussing  lens  and  focal  plane  array 
respectively,  the  centre  of  the  point  spread  function  for  each  spectral  wavelength  can  be  determined 
which  closely  match  the  realistic  physical  locations  expected  from  analysing  the  continuous  fields 
of  the  point  spread  functions. 

3. 7. 1  Real  image  locations.  The  time  delay  Fourier  Transform  relationship  is  analogous 
to  the  relationship  between  these  transverse  planes  using  Fourier  optics.  The  time  delay  Fourier 
Transform  identity  is  stated  in  equation  3.32, 


g(t  -  a)  =>  e-jwQG(w)  (3.32) 

where  G(o>)  =  F{g(t)}  =  g(t)e_jwtdt,  is  the  expression  for  the  Fourier  Transform  used  in  the 
equation.  This  equation  states  that  a  delay  in  time,  a ,  of  the  signal,  g(t),  relates  to  a  constant 
gradient  of  phase  change  over  the  frequency  content  of  the  signal. 

Using  Fourier  optics  to  describe  the  relationship  of  the  electromagnetic  fields  between  both 
transverse  planes,  the  time  delay  Fourier  Transform  relationship  is  modified  to  describe  the  rela¬ 
tionship  between  these  fields,  as  described  at  equation  3.33. 


g(x-  au,y  -  av)  =>  e-i(“"“"+“'a'>GK,Wv)  =  e-^(™u+vQv)Gk;Wv) 


(3.33) 


For  the  detector  optics  system,  this  expression  indicates  that  a  constant  gradient  of  phase 
introduced  to  the  field  that  characterising  the  spatial  frequency  content  for  the  remote  scene  (right 
side  of  equation  3.33),  is  related  to  the  position  offset  location  from  the  central  location  (left  side 
of  equation  3.33)  for  the  image  of  the  remote  scene.  The  expression  g(x  —  au,y  —  av)  represents 
the  field  present  in  the  transverse  plane  for  the  the  detector  focal  plane  array.  The  expression 
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e  j(wuau+wvav)Q^Wui  a;vj  represents  the  field  present  in  the  transverse  plane  at  the  detector  focussing 
lens,  prior  to  propagating  through  the  lens. 

The  prism  optics  system  has  been  described  to  introduce  a  phase  tilt  to  the  electromagnetic 
fields,  based  on  the  angular  direction  change  in  propagation  from  on-axis  propagation  along  the 
optic  axis.  The  prism  optics  system  can  be  approximated  as  a  collection  of  phase  tilts  on  the 
collimated  wavefront  entering  the  prism  optics  system,  where  these  phase  tilts  are  established 
according  to  the  angle  of  dispersion  calculated  using  Snell’s  Law  of  refraction  at  equation  2.5. 
Relating  each  phase  tilt  to  a  position  in  the  detector  focal  plane  array  where  the  real  image  of  the 
remote  scene  is  formed,  serves  the  purpose  for  approximating  the  prism  optics  system. 

3.7.2  Direct  vision  prism  approximation.  The  angular  change  in  direction  of  propagation 
for  the  electromagnetic  field  propagating  through  the  direct  vision  prism  due  to  dispersion  of  each 
spectral  field,  can  be  represented  as  a  constant  gradient  of  phase  (phase  tilt)  across  the  collimated 
wavefront  in  the  direction  of  dispersion. 

Commencing  with  the  electromagnetic  field  at  the  detector  focussing  lens  U5.  This  field  can 
be  represented  as  the  Fourier  Transform  of  the  electromagnetic  field  at  the  detector  focal  plane 
array  that  forms  an  upright  and  minified  real  image  of  the  remote  scene.  The  field  at  the  detector 
focussing  lens,  leaving  the  prism  optics  system,  U5  can  be  described  by  equation  3.34. 


U5(u,v)  =  G(u,  v)  =  -F{g(x,y)} 


(3.34) 


The  two-dimensional  Fourier  Transform  of  an  arbitrary  image  as  a  two-dimensional  signal 
g(x,  y)  ,  representing  the  real  image  at  the  focal  plane  array  for  a  specific  wavelength  of  electro¬ 
magnetic  radiation,  is  described  at  equation  3.35, 


F{  g(x,y)} 


r  /*+oo 

//  g(x,y)e-j(“uX+w''y)dxdy 

J  J  —CO 


(3.35) 
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where  the  angular  spatial  frequency  expressions  uju  and  u>v  represent  the  spectrum  of  electromag¬ 
netic  radiation  in  terms  of  spatial  frequencies,  fu  and  fv,  as,  wu  =  27rfu  and  wv  =  27rfv  .  The 
two-dimensional  Fourier  Transform  expression,  described  at  equation  3.36  with  the  use  of  these 
expressions  becomes, 

r  /*+oo 

F{g(x,y)}=  //  g(x,y)e-j2^f“x+f^Mxdy  (3.36) 

J  J  —  OO 

and  the  expressions  for  these  spatial  frequencies  respectively  are,  fu  =  ^  and  fv  =  .  Sub¬ 

stituting  for  these  expressions,  the  two-dimensional  Fourier  Transform  expression  described  at 
equation  3.37  becomes, 

f{g(x>y)}=//  g(x,y)e_J^(ux+vy)dxdy  (3.37) 

J  J  —  OO 


Examine  the  electromagnetic  field  of  the  same  real  image,  now  offset  a  displacement  x  =  a  and 
y  =  b  from  the  previous  position  in  the  focal  plane  array  transverse  plane.  The  two-dimensional 
Fourier  Transform  of  this  offset  real  image  is  described  at  equation  3.38. 


Gpt(u,v)  =  F{g(x  —  a,  y  —  b)} 


+oo 


g(x  -  a,  y  -  b)e"j  ^(ux+vy)dxdy 


(3.38) 


Grouping  the  system  ordinates  variables  and  the  offset  displacements  as  new  variables,  where 
x  —  a  =  xa  and  y  —  b  =  yb,  and  the  derivative  with  respect  to  the  system  ordinates  variables  as, 
dx  =  dxa  and  dy  =  dyb,  substituting  the  new  variables  for  the  system  ordinate  variables  in  the 
expression  is  described  at  equation  3.39. 


F{g(x-a,y-b)}  =  F{g(xa,yb)} 


r  r  +  CO 

/  /  g(xa,yb)e-jSMxa+a)  +  v(yb+b)}dXadyb 

J  J  —  OO 


(3.39) 
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Rationalising  terms  from  the  integration  expression  leads  to  the  two-dimensional  Fourier 


Transform  expression  as  described  at  equation  3.40. 


F{g(x^a,y-b)} 


p-jS(au+bv) 


+  OO 


,  WjS(uxa+vyb)  rK' 


(3.40) 


Recognising  the  expression  contained  within  the  integral  to  be  similar  to  the  expression  at 
equation  3.35,  changing  the  notation  of  variables  within  this  integral  where  x  =  xa,  y  =  yb  and 
dx  =  dxa,  dy  =  dyb,  the  two-dimensional  Fourier  Transform  expression  of  the  real  image  offset  in 
position,  as  described  at  equation  3.41  becomes, 


F{g(x-a,y-b)} 


=-j  If  (au+bv) 


r  /*+ oo 

//  g(x,y)e-^(ux+vy)dxdy 
•  '  J  —  OO 


(3.41) 


Expressing  the  Fourier  Transform  in  terms  of  the  original  real  image,  as  described  by  equa¬ 
tion  3.42, 


Gpt(u,  v)  =  F{g(x  —  a,  y  —  b)}  =  e-iS(-u+bv)jr{g(x>y)}  =  e-i£(.»+bv)G(U}  v)  (3.42) 

the  Fourier  Transform  for  the  offset  real  image  can  be  expressed  in  terms  the  Fourier  Transform 
of  the  image  at  the  original  location  with  a  constant  phase  gradient  included,  as  described  at 
equation  3.43,  where  the  phase  gradient  is  related  to  the  displacement  offset  between  the  images. 

U6(x  —  a,  y  —  b)  =  g(x-a,y-b)  =►  e-j£<“+bv>  Ub(u,v)  =  g-iCa^+b^)  u5(u,v)  (3.43) 

3.7.3  Phase  tilts  for  prism  rotation.  The  expressions  derived  earlier  use  the  principles 
of  Fourier  Transforms  to  relate  the  two-dimension  electromagnetic  field  between  the  focussing  lens 
and  the  focal  plane  array  in  the  detector  optics  system.  The  software  model  was  designed  to 
use  composite  phase  tilts  of  the  electromagnetic  wavefront  as  the  reference  for  predetermining  the 
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location  of  the  centre  for  each  point  spread  function.  The  composite  phase  tilts  exhibited  by  the 
simulated  electromagnetic  field  wavefronts  along  the  x-axis  and  y-axis  directions  with  respect  to  the 
system  coordinates  are  required  for  every  angular  position  of  rotation  for  the  direct  vision  prism. 
Predetermining  these  locations  with  the  software  model  allows  the  electromagnetic  field  wavefronts 
to  be  realistically  simulated  to  form  real  images  of  the  remote  scene  in  the  transverse  plane  of  the 
detector  focal  plane  array. 

The  angular  direction  of  propagation  with  respect  to  the  optic  axis  64b,  for  each  wavefront  of 
electromagnetic  field  can  be  expressed  in  physical  dimensions  as  described  from  equation  3.44  to 
equation  3.46. 


sin  e4b  = 


a 

c 


(3.44) 


COS  64b  = 


b 

c 


(3.45) 


tan  £44,  =  ^  (3.46) 

b 

Representing  these  angular  directions  or  trajectory  angles  for  each  wavefront  as  composite 
trajectory  angles  £4bx  and  C4bY  and  relate  to  the  phase  tilt  for  each  wavefront  with  respect  to  x-axis 
and  y-axis  of  system  coordinates. 

As  the  prism  rotates,  the  composite  trajectory  angles  C4bx  and  C4bY  are  related  to  the  dis¬ 
persion  trajectory  angle  64b  as  described  in  equations  3.47  to  equation  3.54  as  follows. 


tan  9  = 

b 


(3.47) 
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(3.48) 


tan(|  -  0)  = 


d2 

b 


Using  the  Sine  rule. 


sin  9  sin  f  1 

di  ei  ei 


(3.49) 


sin(|  —  9)  sin  |  1 

9l2  e2  e2 


(3.50) 


Noting  the  composite  angles  C4bx  and  64bY  respectively. 


tan  e4bx  =  —  = 


asin(|  —  9)  asin(|  —  9) 


btan(|  —  9)  b 


=  rcos(x  ~e) 


(3.51) 


a  a  sin  9 
tan  e4bY  =  —  =  — ; - 

ei  di 


a  sin  9 
btan  9 


a  a 
—  cos  9 
b 


(3.52) 


Observing  the  relationship,  cos(|  —  9)  =  —  sin#,  and  revisiting  the  expression  for  the  disper¬ 
sion  trajectory  angle,  tan  64b  =  ?,  the  composite  trajectory  angles  e4bx  and  64b  y  can  be  expressed 
as  follows, 


£4bx  =  arctan{— sin0tane4b}  (3.53) 


£4bY  =  arctan{cos6)tane4b}  (3.54) 
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Apply  these  expressions  to  specifically  determine  the  location  of  the  real  image  in  the  focal 
plane  array  transverse  plane  of  the  object  according  to  the  Dispersion  for  each  wavelength  of  optical 
EM  radiation. 

The  optical  EM  radiation  is  always  examined  across  a  transverse  plane  to  the  optic  axis. 
Examining  the  dispersion  of  optical  EM  radiation  resulting  from  propagation  through  the  prism 
across  a  transverse  plane  leaves  a  constant  change  in  phase  (phase  tilt)  in  the  direction  of  dispersion, 
encoded  across  the  entire  wavefront  of  optical  EM  radiation. 

Examining  the  phase  tilt  in  terms  of  the  composite  trajectory  angles,  the  phase  tilt  across 
the  transverse  plane  along  the  x-axis  direction  in  system  coordinates  is  described  at  equation  3.55, 
and  the  phase  tilt  across  the  transverse  plane  along  the  y-axis  direction  in  system  coordinates  is 
described  at  equation  3.56. 


ej2frx  =ej2fuSine4bx  (3.55) 

ej^ry  =  ej2fvsine4bY  (3.56) 

Using  these  expressions  to  represent  the  composite  phase  tilts  evident  on  the  electromagnetic 
field  wavefront  at  the  detector  focussing  lens,  after  propagating  through  the  prism  optics  system, 
the  expression  for  the  fields  at  the  focussing  lens  U5  and  the  focal  plane  array  U6,  are  described 
from  equation  3.57  to  equation  3.64. 


U5(u,  v)  =  ej^usine4bxej^vsin£4bYF{U6(x,y)} 


(3.57) 


U5(u,  v)  =  e-j(-^usine4bx)e-j(-^vsine4bY)F{U6(x,y)}  (3.58) 
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U5(U,V)  =  e-j(-^usin£^x)e-j(-2fvSin£4bY)  /  /  U6(x,y)e_j(w“X+u'’'y)dxdy 


+  OO 


(3.59) 


/»  /»  -(“OO 

U5(u, v)  =  e-j(-Hfusine4bx)e-j(-2fvsine4bY)  /  /  U6(x,y)e-j2,r(f“x+f''y)dxdy 

J  J  —  OO 


(3.60) 


✓*  /*+oo 

U5(u,v)  =  e-j(-^usine4bx)e-j(-^vsin<r4bY)  /  /  U6(x,y)e-j27r(^x+^y)dxdy 

J  J  —  OO 


(3.61) 


+oo 


U5(u,v)=  //  U6(x,y)e-j2S!(x-zsin£4bx)e-j2S:(y-zsine4bY)dxdy 


(3.62) 


Using  new  variables,  s  =  x  — zsine^x  and  t  =  y  —  z sinewy,  together  with  the  respective 
derivatives,  els  =  dx  and  dt  =  dy,  the  expression  becomes, 


+oo 


U5 (u, v)  =  //  Ue(s  +  z sin t  +  z sin e4by)e  “zs+ Az^dsdt  (3.63) 


U5(u,  v)  =  U{U6(s  +  zsine4bx,t  +  zsine4bY)}  (3.64) 

The  expression  at  equation  3.64  offers  the  ability  to  locate  the  real  image  of  the  point  spread 
functions  in  the  focal  plane  array,  due  to  the  angle  of  propagation  due  to  dispersion  through  the 
prism.  This  expression  together  with  equation  3.53  and  equation  3.54  provide  a  mathematical 
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outcome  that  was  used  for  the  software  model  to  commence  locating  the  centres  of  the  real  images 
for  the  point  spread  functions. 
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IV.  Analysis 


This  chapter  describes  the  main  considerations  encountered  with  designing  the  software  model 
to  simulate  the  sensor  optics  system.  Considerations  with  the  area  of  each  aperture,  at 
section  4.1,  the  optical  component  parameters,  at  section  4.2,  and  numerically  simulating  wave 
optics  propagation,  at  section  4.3  are  discussed.  Manipulation  of  the  MatLab®  functions  and 
program  scripts  comprising  the  suite  of  software  developed  during  this  research  is  presented  at 
section  4.5  to  demonstrate  the  design  of  a  software  model  that  simulates  a  particular  configuration 
of  the  sensor  optics  system. 

1 . 1  Apertures 

The  spatial  extent  of  the  apertures  throughout  the  sensor  optics  system  should  be  validated 
before  using  these  dimensions  for  each  respective  aperture  in  the  construct  of  the  simulated  sensor 
optics  system  software  model.  An  explanation  of  the  sequence  of  apertures  selected  for  validating 
the  software  model,  each  numbered  to  enable  the  reader  to  keep  track  of  the  dependencies  between 
certain  apertures,  appears  in  the  following  paragraphs. 

Based  on  the  area  of  the  direct  vision  prism  entry  aperture  (1),  the  area  of  the  aft  lens  aperture 
(2)  was  established  for  some  spatial  freedom  to  accommodate  the  spectral  fields  of  electromagnetic 
radiation  dispersing  through  the  prism.  Validating  the  area  of  the  direct  vision  prism  exit  aperture 
(3)  was  based  on  the  range  of  spectral  wavelengths  of  electromagnetic  radiation  dispersed  through 
the  direct  vision  prism.  Each  spectral  held  in  the  software  model  represents  the  respective  wave¬ 
length  band  of  electromagnetic  radiation  collected  by  the  physical  sensor  optics  system.  Validating 
the  area  of  the  detector  focussing  lens  aperture  (4)  was  based  on  the  distance  of  separation  along 
the  optic  axis  between  the  direct  vision  prism  exit  surface  and  the  detector  focussing  lens.  The 
angular  directions  propagation  for  every  spectral  wavefront  held  of  electromagnetic  radiation  must 
be  collected  within  the  spatial  extent  of  the  focussing  lens  aperture. 
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To  simulate  the  sensor  optics  system  these  dependencies  between  system  parameters  required 
validation  before  a  software  model  configuration  could  be  used.  Validation  of  the  system  parameters 
to  be  used  in  the  software  model  was  based  on  the  collective  influence  of  the  respective  areas  of 
each  aperture,  the  separation  between  the  direct  vision  prism  and  detector  focussing  lens,  the 
choice  of  undeviated  spectral  wavelength  of  electromagnetic  radiation,  the  finite  spectral  range  of 
electromagnetic  radiation  of  interest  for  the  design  of  the  sensor  optics  system.  Geometric  optics 
was  used  in  the  software  model  to  discern  the  spatial  extent  requirements  for  each  aperture  and 
validate  the  system  parameters  used  for  the  simulated  sensor  optics  system. 

4-1.1  Forward  Optics  System.  The  spatial  extent  of  the  forward  aperture  and  lens  should 
be  selected  to  allow  the  collection  of  sufficient  electromagnetic  radiation  from  the  remote  scene  for 
hyperspectral  imaging.  The  forward  aperture  collects  (sampling)  the  radiant  energy  incident  on  the 
sensor  optic  system  from  the  remote  scene,  represented  as  spectral  electromagnetic  field  wavefronts 
in  the  software  model,  and  the  forward  lens  focusses  focusses  this  radiant  energy  to  form  a  real 
image  of  the  remote  scene  in  the  field  stop.  Spatial  sampling  criteria  is  important  for  the  software 
model  to  ensure  continuous  spectral  fields  are  represented  accurately  as  discrete  spectral  fields  with 
the  software  model.  Spatial  sampling  at  the  forward  lens  establishes  the  representation  of  the 
spectral  fields  of  electromagnetic  radiation  used  through  the  remaining  simulated  propagation  of 
these  fields. 

The  spatial  extent  of  the  aft  aperture  and  lens  is  required  to  be  smaller  than  the  spatial 
extent  of  the  direct  vision  prism  entry  aperture.  The  spatial  extent  of  the  aperture,  and  number 
characteristics  of  the  aft  lens  are  established  in  the  software  model  so  the  radiant  energy  collected 
by  the  forward  aperture  also  propagates  through  the  aft  aperture. 

To  reduce  the  wavefront  error  resulting  from  the  propagation  of  the  electromagnetic  radi¬ 
ation  through  both  lenses,  the  software  model  represents  both  lenses  with  the  same  F #  number 
characteristics. 
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4-1-2  Prism  Optics  System.  The  spatial  extent  of  the  direct  vision  prism  exit  aperture 
must  be  sufficient  to  accommodate  for  dispersing  all  spectral  wavefront  fields  of  electromagnetic 
radiation  from  the  remote  scene  to  be  used  in  the  hyperspectral  real  image  data  set.  The  spatial 
offset  in  the  transverse  plane  and  angular  direction  of  propagation,  both  with  respect  to  the  optic 
axis  for  each  spectral  wavefront  differs  according  to  Snell’s  Law  of  Refraction  and  the  bulk  media 
configurations  used  in  the  direct  vision  prism  for  each  stage  of  propagation  in  the  prism  optics 
system.  Each  spectral  wavelength  no  longer  propagates  symmetrically  along  the  optic  axis,  instead 
each  spectral  wavelength  has  some  distinct  angular  direction  of  propagation  with  respect  to  the 
optic  axis.  In  the  software  model,  the  presence  of  electromagnetic  radiation  at  the  outermost 
locations  with  respect  to  the  optic  axis,  in  the  transverse  plane  of  the  exit  surface  of  the  direct 
vision  prism  defines  the  minimum  aperture  area  required  for  the  direct  vision  prism. 

The  spatial  extent  of  the  detector  focussing  lens  must  be  sufficient  to  collect  the  incident 
fields  of  electromagnetic  radiation  dispersed  from  propagation  through  the  direct  vision  prism. 
The  separation  along  the  optic  axis  between  the  exit  surface  of  the  direct  vision  prism  and  the 
detector  focussing  lens  influences  the  required  spatial  extent  of  the  detector  focussing  lens.  In  the 
software  model,  the  presence  of  electromagnetic  radiation  at  the  outermost  locations  with  respect 
to  the  optic  axis,  in  the  transverse  plane  of  the  focussing  lens  defines  the  minimum  aperture  area 
required  for  the  focussing  lens. 

4-1.3  Detector  Optics  System.  The  detector  focussing  lens  used  in  the  software  model  to 
focus  the  electromagnetic  radiation  collected  through  the  focussing  lens  aperture  on  the  detector 
focal  plane  array,  is  configured  with  the  same  F #  number  characteristics  as  the  previous  lenses  in 
the  forward  optics  system. 
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4-2  Optical  Components 


The  parameters  used  in  the  software  model  to  represent  the  optical  components  of  the  sensor 
optics  system  were  selected  to  enable  use  of  equations  based  on  the  theory  of  optics  stemming  from 
geometric  optics,  Fourier  optics  and  the  Rayleigh-Sommerfield  diffraction  formula.  The  system 
components  represented  in  the  software  model  were  the  forward  lens,  the  aft  lens,  the  direct  vision 
prism,  and  the  detector  focussing  lens.  All  intrinsic  characteristics  of  the  detector  focal  plane  array 
such  as  material  sensivities,  non-uniformities,  responsivity,  are  not  represented  in  the  software 
model.  The  spectral  fields  of  electromagnetic  radiation  incident  on  the  detector  focal  plane  array 
are  computed  in  the  software  model. 

4-2.1  Lenses.  All  lenses  in  the  sensor  optics  system  were  represented  with  the  Lens 
Maker’s  equation  in  the  software  model.  These  lenses  are  approximated  as  thin  lenses  with  this 
equation  also  known  as  the  Gauss  thin  lens  formula,  and  using  this  formula  implies  the  constraint 
with  choice  of  numbers  to  characterise  the  lenses.  The  maximum  thickness  and  the  focal  length 
of  each  lens  are  two  parameters  used  in  the  Lens  Maker’s  equation,  where  the  maximum  thickness 
was  arbitrarily  chosen  for  each  lens,  and  the  focal  length  defined  so  the  same  F#  number  was 
used.  The  choice  of  F^  number  used  for  all  lenses  in  the  software  model  obeying  the  criteria 
of  F#  number  >  10. 

The  software  model  simulates  the  point  spread  function  of  the  sensor  optics  system  for  a  point 
source  located  at  some  distance  from  the  forward  lens  of  the  system  emitting  spectral  wavefront 
fields  as  spherical  wavefronts  of  electromagnetic  radiation.  Each  incident  spectral  wavefront  has  a 
different  phase  curvature  across  the  aperture  of  the  forward  lens,  resulting  with  the  real  image  of 
the  point  spread  function  for  each  spectral  wavefront  located  at  different  locations  along  the  optic 
axis.  If  these  spectral  wavefronts  appeared  planar  (hence  no  phase  curvature)  across  the  aperture 
of  the  forward  lens,  the  real  image  of  the  point  spread  functions  of  each  spectral  wavefront  would 
be  located  at  the  focal  point  of  the  lens. 
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Figure  4.1:  Spherical  wavefront  propagation  through  a  lens  [2:  164].  The  point  source  loca¬ 
tion  along  the  optic  axis  from  the  lens  influences  where  on  the  optic  axis  the  spherical  wavefront 
converges  to  form  a  real  image  of  the  point  source. 

The  forward  lens  aperture  would  be  arbitrarily  chosen  to  collect  sufficient  radiant  electromag¬ 
netic  energy  and  enable  recording  of  spectral  real  image  data  from  the  remote  scene.  The  aft  lens 
aperture  would  be  chosen  to  ensure  this  energy  is  permitted  to  propagate  and  disperse  through  the 
prism  optics  system.  The  distance  between  the  remote  scene  and  the  forward  aperture  influences 
the  angular  extent  of  the  emitted  spherical  wavefront  incident  on  the  forward  aperture,  and  further 
influences  the  focal  distance  of  the  real  image  for  the  point  source  at  the  field  stop.  So  the  location 
of  the  remote  scene  cannot  be  within  the  focal  length  of  the  forward  lens.  Figure  4.1  illustrates 
the  location  of  the  point  source  at  various  distances  from  the  forward  aperture  and  lens  that  result 
in  real  images  being  formed.  Observe  from  this  figure  the  location  of  the  remote  scene  cannot  be 
within  the  focal  length  of  the  forward  lens  in  order  for  a  real  image  of  the  remote  scene  to  be  formed 
in  the  field  stop. 

Choosing  the  same  F^  number  for  all  lenses  offers  an  inherent  attribute  of  minimising  the 
wavefront  error  encountered  on  the  spectral  wavefront  fields  of  electromagnetic  radiation  with 
propagating  through  lenses.  Near  complementary  magnitude  and  phase  errors  due  to  propagation 
through  the  forward  optics  system  are  introduced  to  the  spectral  wavefront  fields  when  the  same 
F#  number  is  used  for  the  forward  and  aft  lenses,  thereby  minimising  wavefront  error.  Figure  4.2 
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Figure  4.2:  Point  source  location  along  the  optic  axis  of  the  sensor  optics  system.  The  location  of 
the  remote  scene  is  defined  by  the  location  of  the  point  source  in  the  software  model.  The  location 
of  the  forward  lens  with  respect  to  the  held  stop  is  related  to  the  location  of  the  point  source  along 
the  optic  axis. 


illustrates  the  arrangement  of  the  forward  and  aft  lenses  used  in  the  software  model  to  obtain  a 
real  image  of  the  point  source  in  the  held  stop,  and  propagate  the  spectral  wavefront  helds  to  the 
aft  lens. 


4-2.2  Direct  Vision  Prism.  Optical  materials  chosen  to  represent  the  bulk  media  and 
thin  coating  stages  of  propagation  through  the  direct  vision  prism  were  selected  because  they 
exhibited  high  transmissive  properties  over  the  medium  wave  infrared  (MWIR)  spectral  range 
of  electromagnetic  radiation.  The  indexes  of  refraction  for  these  materials  at  specihc  spectral 
wavelengths  were  obtained  and  linear  interpolation  between  each  specihed  spectral  wavelength  was 
used  in  the  software  model  to  develop  index  of  refraction  data  for  these  materials  in  one  nanometre 
increments.  Indexes  of  refraction  characteristics  of  these  optical  materials  suited  to  the  MWIR 
spectral  range  of  electromagnetic  radiation  and  available  in  the  software  model  are  detailed  at 
Appendix  B.  Any  birefringent  properties  with  these  materials  were  treated  in  the  software  model 
only  for  the  fast  axis  indexes  of  refraction. 
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4-3  Wave  Optics  Propagation 


The  propagation  of  each  spectral  field  wavefront  of  electromagnetic  radiation  through  the 
lenses  and  direct  vision  prism  was  represented  as  changes  in  phase  to  each  piecewise  sample  of 
these  wavefronts  in  the  software  model.  These  optical  components  were  treated  with  transmissive 
properties  of  100%  transmission  for  the  MWIR  spectrum  of  electromagnetic  radiation. 

All  spectral  fields  of  electromagnetic  radiation  from  the  multispectral  point  source  were  rep¬ 
resented  as  spherical  wavefronts  with  the  Rayleigh-Sommerfield  Diffraction  theory.  Equation  2.7 
describe  this  theory  required  some  manipulation  to  represent  the  discrete  samples  of  the  spectral 
fields  of  electromagnetic  radiation  through  each  stage  of  spherical  wavefront  propagation  in  the 
software  model  of  the  sensor  optics  system.  The  stages  of  spherical  wavefront  propagation  evident 
in  the  software  model  were:  (1)  from  the  multispectral  point  source  location  to  the  forward  lens 
as  a  spatial  diverging  wavefront;  (2)  from  the  forward  lens  to  the  field  stop  as  a  spatial  converging 
wavefront;  (3)  from  the  field  stop  to  the  aft  lens  as  a  spatial  diverging  wavefront;  and  (4)  from  the 
focussing  lens  to  the  focal  plane  array  within  the  detector  as  a  spatial  converging  wavefront. 

The  numeric  form  of  the  Rayleigh-Sommerfield  Diffraction  formula  used  to  represent  the 
effects  of  diffraction  across  the  wavefront  each  spectral  field  wavefront  due  to  spherical  wavefront 
propagation,  developed  further  from  equation  3.16,  is  described  at  equation  4.1, 

i(Nx-l)  |(Ny  — 1) 

Ui(xi,yi,A)  =  Ai(xi,yi,  A)  ^  ^  °  *n^Jny ’  eJ ^ r  tiens (x„x ,  yny ,  A)  (4.1) 

„x=-i(Nx-l)ny=-i(Ny-l) 

where  the  discrete  amplitude  function,  A1(x1,y1,A),  is  described  in  equation  4.2,  the  discrete 
path  length  of  propagation  function,  r  =  r(xi,yi,xnx,yny),  is  described  in  equation  4.3,  the  lens 
transformation  function,  tiens(xnx, yny,  A),  and  lens  thickness  function,  A(x„x,yny),  are  described 
in  equation  4.4. 

Ai(xi,yi,  A)  =  j^sinc  (-^ydx0)  sine  (-^ydy0)  (4-2) 
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xi)2  +  (yny  -  yi)2 


(4.3) 


r 


r(xi,yi,xnx,yny)  =  yz2  +  (s 


^lens 


,yn  ,A)  =  ej  Tc(nlens-l)A(x„x,y„y)  _  gj  ^(niens-l)  A0  -  (xQ2 +y02)  (  3^- -  jjA-)  ] 


(4.4) 


4-4  Design  the  Sensor  Optics  System 

The  software  suite  used  for  developing  the  software  model  of  the  simulated  sensor  optics 
system  is  listed  at  Appendix  C.  Knowing  certain  parameters  about  the  optical  components  to 
be  used  for  the  sensor  optics  system  is  only  the  start  for  developing  the  software  model.  Other 
parameters  such  as  spatial  sampling  intervals,  separation  distances  along  the  optic  axis  between 
components,  validation  of  aperture  sizes  are  derived  from  the  optical  components  characteristics 
to  form  the  software  model  and  simulate  the  design  of  the  sensor  optics  system.  One  program 
script  was  used  as  the  software  model  to  simulate  the  wave  optics  propagation  of  a  point  source  of 
electromagnetic  radiation  through  the  sensor  optics  system  appears  at  Appendix  C. 

4-4-1  System  parameters.  The  parameters  of  the  software  model  are  defined  and  evalu¬ 
ated  using  the  MatLab®  function  parameters,  m.  The  spectrum  of  electromagnetic  radiation  used 
for  representing  the  remote  scene  is  represented  an  array  of  discrete  wavelengths,  where  each  wave¬ 
length  is  the  centre  wavelength  for  a  specific  spectral  wavelength  band  of  electromagnetic  radiation. 
Parameters  that  are  specified  for  the  system  design  are  represented  in  the  MatLab®  function  pa¬ 
rameters  -02-define.m,  and  based  on  these  parameters,  other  system  parameters  computed  for  the 
system  design  are  evaluated  with  the  MatLab®  function  parameterS-03-Compute.m.  The  distances 
along  the  optic  axis  between  the  forward  lens  and  the  field  stop,  and  between  the  field  stop  and 
the  aft  lens,  plus  the  spatial  sampling  intervals  for  each  principal  transverse  plane  in  the  simulated 
sensor  optics  system  are  computed  based  on  the  defined  system  parameters. 

4-4 -2  Undeviated  spectral  wavelength.  Two  pivotal  parameters  in  the  system  design  are 
the  undeviated  spectral  wavelength  and  the  spectral  range  of  electromagnetic  radiation  around  this 
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wavelength  required  to  represent  the  spectral  detail  of  radiant  energy  collected  from  the  remote 
scene.  Both  parameters  are  specified  by  the  system  designer,  and  both  contribute  to  the  dispersion 
characteristics  for  the  system.  From  the  specified  diameter  of  the  direct  vision  prism,  the  dispersion 
characteristics  influence  the  aperture  sizes  required  at  the  aft  lens  (forward  optics  system)  and  the 
focussing  lens  (detector  optics  system). 

The  design  of  the  direct  vision  prism  must  be  tailored  for  the  undeviated  spectral  wavelength, 
since  the  spectral  real  images  of  the  remote  scene  are  spatially  distributed  along  a  radial  line 
(representing  the  direction  of  dispersion)  in  the  transverse  plane  at  the  detector  focal  plane  array. 
The  undeviated  spectral  wavelength  spectral  real  image  is  formed  with  the  centre  of  the  image  on 
the  optic  axis.  Spectral  real  images  for  spectral  wavelength  bands  below  the  undeviated  spectral 
wavelength  are  formed  with  centres  along  the  radial  line  located  away  from  the  optic  axis  in  one 
direction,  and  the  images  for  spectral  wavelength  bands  above  the  undeviated  spectral  wavelength 
are  formed  in  the  opposite  direction  along  the  radial  line.  The  spatial  separation  between  each 
spectral  real  image  is  derived  from  the  dispersion  attained  with  the  direct  vision  prism  for  each 
corresponding  spectral  wavelength  band. 

Using  Snell’s  Law  of  Refraction,  the  undeviated  spectral  wavelength  of  electromagnetic  ra¬ 
diation  is  evaluated  using  geometric  optics  through  each  stage  of  propagation  through  the  prism 
optics  system,  to  arrive  at  the  detector  focussing  lens  with  an  angular  direction  of  propagation 
with  respect  to  the  optic  axis  of  zero.  The  parameter  64b  is  used  in  the  software  model  to  describe 
this  angular  direction.  Based  on  the  index  of  refraction  characteristics  of  the  bulk  media  materials 
used  for  the  direct  vision  prism,  the  parameters  of  the  direct  vision  prism  available  for  tailoring  the 
system  design  are  the  separate  prism  surface  angles.  By  retaining  one  prism  surface  angle  constant, 
and  thus  a  specified  parameter,  the  other  prism  surface  angle  can  be  evaluated  iteratively  by  testing 
against  the  angular  direction  of  propagation  criterion,  so  the  prism  surface  angle  solution  achieves 
an  angular  direction  of  propagation  with  respect  to  the  optic  axis  of  zero. 
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The  entry  prism  surface  angle,  defined  as  ai,  can  be  evaluated  by  using  the  MatLab®  function 
evaluate-alpha-1  .m,  with  the  exit  prism  surface  angle  predefined.  The  exit  prism  surface  angle, 
defined  as  03,  can  be  evaluated  by  using  the  MatLab®  function  evaluate-alphaS.m,  with  the  entry 
prism  surface  angle  predefined. 

4-4-3  Spatial  extent  of  apertures.  With  the  surface  angles  of  the  direct  vision  prism 
established  so  the  prism  optics  system  is  tailored  for  the  specified  undeviated  spectral  wavelength, 
the  next  task  in  designing  the  software  model  is  to  validate  the  spatial  extent  of  all  apertures  in  the 
prism  optics  system.  These  apertures  are:  (1)  the  prism  entry  surface;  (2)  the  prism  exit  surface; 
(3)  the  aft  lens;  and  (4)  the  detector  focussing  lens. 

The  direction  of  dispersion  encountered  by  the  spectral  wavefront  fields  of  electromagnetic 
radiation  occurs  along  the  vertical  orientation  in  prism  coordinates.  Using  an  upright  orientation 
of  the  direct  vision  prism,  the  direction  of  dispersion  for  these  spectral  wavefront  fields  also  occurs 
in  the  vertical  orientation  in  system  coordinates. 

Geometric  optics  and  Snell’s  Law  of  Refraction  are  used  to  evaluate  the  propagation  path 
distances  through  each  stage  of  propagation  in  the  prism  optics  system.  The  minimum  lengths 
parallel  to  the  optic  axis  of  the  direct  vision  prism,  defined  Z2min  for  the  entry  surface  and  Z3min 
for  the  exit  surface,  plus  the  diameter  of  the  direct  vision  prism,  are  system  parameters  defined 
by  the  system  designer.  Working  with  these  parameters,  sensible  propagation  path  distances  can 
be  determined.  By  testing:  (1)  the  propagation  distances  with  the  minimum  lengths  parallel  to 
the  optic  axis;  and  (2)  the  vertical  position  of  each  spatial  sample  of  the  spectral  wavefront  fields 
as  these  positions  change  due  to  dispersion  with  the  radius  of  the  direct  vision  prism,  sensible 
spatial  extent  sizes  for  the  entry  surface  and  exit  surface  apertures  of  the  direct  vision  prism  can 
be  determined  by  the  system  designer.  Once  sensible  minimum  spatial  extent  sizes  for  the  direct 
vision  prism  apertures  are  determined,  the  sensible  minimum  spatial  extent  size  for  the  detector 
focussing  lens  aperture  can  be  determined. 
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Should  the  tests  for  propagation  distances  or  vertical  positions  fail  during  these  evaluations, 
the  system  designer  has  the  option  of  changing  the  defined  system  parameters  of:  (1)  direct  vision 
prism  diameter,  to  effect  a  change  in  entry  surface  and  exit  surface  apertures;  (2)  minimum  lengths 
of  each  bulk  media  component  of  the  direct  vision  prism  diameter,  to  effect  changes  in  propagation 
distance  and  vertical  position  to  accommodate  for  the  available  spatial  extent  of  the  entry  surface 
and  exit  surface  apertures;  (3)  separation  between  the  direct  vision  prism  exit  aperture  and  the 
detector  focussing  lens  along  the  optic  axis,  defined  z5ref,  to  effect  a  change  in  the  spatial  extent  for 
all  dispersed  spectral  wavefront  fields  collected  at  the  focussing  lens  aperture;  (4)  diameters  of  the 
aft  lens  and  focussing  lens,  to  effect  adjustments  in  the  spatial  extent  of  the  spectral  wavefront  fields 
commencing  propagation  through  the  prism  optics  system  (with  the  aft  lens)  or  effect  adjustments 
in  the  spatial  extent  available  at  the  detector  for  collecting  these  fields  (focussing  lens). 

Deciding  how  these  defined  system  parameters  should  be  varied  is  left  to  the  system  designer, 
since  each  parameter  is  not  entirely  independent  from  influencing  the  design  requirements  with 
another  of  these  parameters.  Each  design  change  made  needs  to  be  tested  against  the  propagation 
distances  and  vertical  position  criteria  so  that  a  sensible  arrangement  of  apertures  in  the  prism 
optics  system  is  found  for  the  software  model. 

Propagation  of  the  spectral  wavefront  fields  of  electromagnetic  radiation  through  the  prism 
optics  system,  is  simulated  with  the  MatLab®  function  propagate-04-prism.m.  The  prism  group 
of  metrics  obtained  from  this  function  deliver  values  for  each  stage  of  propagation  as:  propagation 
distances,  defined  prism. r_l,  prism. r_2,  prism. r_3,  prism. r_4,  prism. r_5  respectively;  and  vertical 
position,  defined  prism.yl .total;  prism. y2_total;  prism. y3_total;  prism. y4_total;  prism.y5_total  re¬ 
spectively.  These  metrics  are  available  for  the  system  designer  to  determine  the  appropriate  choices 
in  system  parameters  for  achieving  a  sensible  arrangement  of  apertures  in  the  system  design.  The 
MatLab®  function  propagate-04-prism.m  can  be  used  with  manipulation  from  the  system  designer 
to  evaluate  different  choices  of  defined  system  parameters. 
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Table  4.1:  Software  Model  Design  Parameters  used  to  simulate  the  CTIS  Sensor  Optics  System 


Description  of  Design  Parameter 

Values 

Number  of  spectral  wavelength  bands 

31 

Centre  wavelengths  of  spectral  wavelength  bands 

2.0  /jin  to  5.0  /iin  in 
100  nm  intervals 

Undeviated  spectral  wavelength 

3.6  /jin 

Number  of  spatial  samples  across  the  remote  scene 

66 

Power,  emitted  for  each  spectral  wavelength  band 

1  Watt 

Range,  distance  between  remote  scene  and  forward  lens 

1000  m 

Forward  Lens  Aperture 

0.1  m 

Forward  Lens  Diameter 

0.1  m 

Forward  Lens  Focal  Length 

1.0  m 

Aft  Lens  Aperture 

0.022  m 

Aft  Lens  Diameter 

0.022  m 

Aft  Lens  Focal  Length 

0.220  m 

Focussing  Lens  Aperture 

0.050  m 

Focussing  Lens  Diameter 

0.050  m 

Focussing  Lens  Focal  Length 

0.500  m 

Prism  Diameter 

0.030  m 

Prism  Entry  Surface  Angle,  op 

30°  00'  00" 

Prism  Exit  Surface  Angle,  a 3 

23°  56'  00" 

Prism  Entry  Bulk  Media,  n3 

Lithium  Fluoride  (LiF) 

Prism  Exit  Bulk  Media,  n.4 

Barium  Fluoride  (BaF2) 

4-5  Simulate  the  Sensor  Optics  System 

Once  a  sensible  arrangement  for  propagation  distances  and  vertical  offsets  were  obtained  for  a 
set  of  system  parameters,  these  parameters  were  retained  to  define  the  sensor  optics  system  design 
configuration  to  be  represented  through  simulation  with  the  software  model.  The  software  model 
was  used  to  propagate  a  multispectral  point  source  through  the  system  and  obtain  the  point  spread 
functions  for  the  specific  system  design  configuration  using  wave  optics  propagation  theory.  From 
the  linear  systems  viewpoint  these  point  spread  functions  represent  the  impulse  response  of  the 
simulated  sensor  optics  system.  The  system  design  configuration  parameters  used  for  the  software 
model  presented  in  the  following  discussion  are  listed  at  Table  4.1.  A  limited  extent  of  the  output 
parameters  available  from  the  software  model  are  presented  in  the  following  discussion,  as  a  means 
of  providing  some  insight  on  the  extent  of  metrics  computed  in  simulating  the  sensor  optics  system. 
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Figure  4.3:  Point  Spread  Functions  for  2.0  /an,  3.0  /im,  4.0  ^m,  5.0  /jm  at  the  Field  Stop  Aper¬ 
ture.  An  increase  in  spectral  wavelength  is  coupled  with  an  increase  in  spatial  area  and  a  decrease 
in  peak  intensity  for  the  Point  Spread  Functions 


4-5.1  Field  Stop  Aperture  (Forward  Optics  System).  The  intensity  of  electromagnetic 
radiation  at  each  spatial  sample  in  the  field  stop  transverse  plane  of  the  forward  optics  system  is 
represented  at  Figure  4.3.  The  point  spread  functions  at  the  field  stop  increase  in  area  across  the 
transverse  plane  for  increased  spectral  wavelength.  Since  the  power  for  each  spectral  wavelength 
band  emitted  from  the  multispectral  point  source  is  the  same,  the  trend  of  a  decrease  in  peak 
intensity  for  increased  spectral  wavelength.  These  complementary  relationships  reveal  the  total 
energy  for  each  spectral  wavelength  band  from  the  point  source  is  conserved,  thus  demonstrating 
the  software  model  also  obeys  the  laws  of  conservation  of  energy. 
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The  cross-section  cut-away  through  the  centre  of  the  point  spread  functions  for  the  2.0  pm, 
3.0  pm,  4.0  pm,  and  5.0  pm  spectral  wavefront  fields  of  electromagnetic  radiation  are  superimposed 
on  the  same  diagram  for  comparison  at  Figure  4.4.  The  ability  to  view  the  spatial  profile  of 
individual  spectral  point  spread  functions  at  the  field  stop  is  available  with  the  software  model, 
however  this  observation  would  not  be  possible  with  the  actual  system.  Observing  the  spatial 
profile  of  the  sum  of  all  spectral  point  spread  functions  at  the  field  stop  would  only  be  possible 
with  the  actual  system  since  the  centres  of  these  spectral  point  spread  functions  are  collocated  in 
the  held  stop. 

4-5.2  Aft.  Lens  Aperture  (Forward  Optics  System).  The  phase  curvature  of  the  spectral 
wavefront  fields  of  electromagnetic  radiation  after  propagating  through  the  aft  lens  are  displayed 
at  Figure  4.5.  Each  spectral  wavefront  held  has  less  than  half  a  wavelength  of  phase  deviation 
across  the  aft  lens  aperture.  These  spectral  wavefront  helds  form  the  inputs  to  the  prism  optics 
system.  The  wavefront  propagation  from  the  aft  lens  aperture  is  simulated  as  a  collimate  wavefront. 
The  low  variation  in  phase  across  the  wavefront  simulated  through  wave  optics  propagation  theory 
supports  this  next  stage  of  propagation  for  the  spectral  wavefront  helds. 

4-5.3  Focussing  Lens  Aperture  (Prism  Optics  System).  Results  from  the  software  model 
for  simulating  the  wave  optics  propagation  through  the  direct  vision  prism  for  3.6  /tin  and  3.8  /tin 
spectral  wavefront  helds  of  electromagnetic  radiation  are  shown  at  Figure  4.6.  The  phase  tilt 
evident  on  each  spectral  wavefront,  for  the  prism  angular  position  of  rotation  oriented  in  the 
upright  position.  The  3.8  gm  spectral  held  shows  noticeable  phase  tilt  across  the  wavefront  along 
the  direction  of  dispersion,  which  is  vertical  due  to  the  angular  orientation  of  the  prism.  The 
3.6  /jm  spectral  held  shows  no  significant  phase  tilt,  demonstrating  this  spectral  wavelength  is  the 
undeviated  wavelength  of  the  simulated  sensor  optics  system. 

Results  from  the  software  model  for  simulating  rotation  of  the  3.8  /mi  wavelength  band  of 
spectral  wavefront  held  at  the  detector  focussing  lens  aperture  are  shown  in  Figure  4.7.  These 
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Field  Intensities  at  the  Field  Stop 


Figure  4.4:  Point  Spread  Function  cross  sections  for  2.0  pm,  3.0  pm,  4.0  pm,  5.0  pm  at  the  Field 
Stop.  Shows  the  relationship  between  peak  intensity  and  the  spatial  cross  section  coverage  (spatial 
area)  for  different  spectral  wavelengths,  where  2.0  pm  has  maximum  peak  intensity  and  5.0  pm  has 
most  spatial  cross  section  coverage. 


results  show  the  sampled  phases  obtained  through  simulated  wave  optics  propagation  using  the 
software  model.  The  trend  of  phase  variation  reveals  the  phase  tilt  across  the  wavefront  aligns  with 
the  angular  orientation  of  the  direct  vision  prism  with  respect  to  the  system  coordinates.  Note  the 
change  location  at  the  detector  focussing  lens  aperture  of  each  sampled  element  of  the  wavefront 
with  different  angular  orientations  of  the  direct  vision  prism.  The  phase  tilts  across  the  wavefront, 
the  angular  orientation  of  the  wavefront  with  respect  to  the  system  coordinates,  and  the  change  in 
spatial  location  of  the  wavefront  with  different  angular  orientations  of  the  prism,  are  results  that 
demonstrate  the  use  of  wave  optics  propagation  theory  in  the  software  model. 
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Figure  4.5:  Phase  of  the  spectral  wavefront  fields  for  2.0  /an,  3.0  pm,  4.0  pm,  5.0  pm  at  the  Aft 
Lens.  Each  spectral  wavefront  field  has  less  than  half  a  wavelength  of  phase  deviation  across  the 
aft  lens  aperture. 

4-5-4  Focal  Plane  Array  (Detector  Optics  System).  Obtaining  the  point  spread  functions 
for  each  spectral  wavelength  of  electromagnetic  radiation  at  the  focal  plane  array  was  achieved 
at  first  for  a  stationary  upright  oriented  direct  vision  prism  configuration  of  the  simulated  sensor 
optics  system.  Propagation  theory  based  on  Fourier  optics  was  used  to  identify  locations  in  the 
transverse  plane  at  the  focal  plane  array  where  the  centres  of  every  spectral  point  spread  function 
were  anticipated.  A  spatial  area  of  arbitrary  size  around  the  location  of  these  spectral  point  spread 
function  centres  was  chosen  to  simulate  the  wave  optics  propagation  to  the  focal  plane  array  for  each 
respective  spectral  field.  The  outcome  from  performing  this  simulation  is  shown  in  Figure  4.8  (a) 
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Spectral  Wavelength  =  3.8  pm  Phase  at  Detector  Focussing  Lens  [  radians  ] 
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Figure  4.6:  Phase  of  the  spectral  wavefront  field  for  wavelengths  (a)  3.6  /jin  and  (b)  3.8  /jin,  with 
the  prism  oriented  in  the  upright  position  with  respect  to  the  system  coordinates  of  the  software 
model.  These  results  from  the  software  model  show  noticeable  phase  tilt  across  the  wavefront  for 
3.8  /jm  spectral  field,  and  no  noticeable  phase  tilt  across  the  wavefront  for  3.6  /(in  spectral  field. 
The  3.6  /(in  spectral  wavefront  field  corresponds  to  the  undeviated  spectral  wavelength  for  this 
configuration  of  the  simulated  sensor  optics  system. 


with  a  spectral  point  spread  function  formed  according  to  the  propagation  of  the  corresponding 
spectral  field.  The  spectral  fields  of  infrared  electromagnetic  radiation  were  uniformly  distributed 
over  the  finite  spectral  range  between  2.0  /im  and  5.0  /im  in  100  nm  increments.  The  spatial 
separation  between  the  respective  spectral  point  spread  functions  in  the  focal  plane  array  was 
not  uniform  along  the  radial  line  aligned  in  the  direction  of  dispersion  and  indicating  the  angular 
orientation  of  the  direct  vision  prism.  This  non-uniform  spatial  separation  between  spectral  point 
spread  functions  is  a  testament  to  the  trigonometric  expressions  at  equation  3.53  and  equation  3.54. 

Incorporating  the  rotating  direct  vision  prism  feature  in  simulating  the  sensor  optics  system 
was  the  next  step  in  designing  the  software  model.  Fourier  optics  propagation  theory  was  used 
to  identify  the  spectral  point  spread  function  centre  locations  for  an  upright  orientation  of  the 
prism  commensurate  with  each  angular  orientation  of  the  direct  vision  prism.  With  a  rotational 
transform  relationship  from  linear  algebra  for  every  discretely  sampled  angular  position  of  prism 
rotation,  these  spectral  point  spread  function  centre  locations  were  used  to  compute  the  anticipated 


81 


Figure  4.7:  Phase  of  the  3.8  /an  spectral  wavefront  field  at  the  detector  focussing  lens  aperture. 
Four  angular  orientations  of  the  3.8  /jm  spectral  wavefront  field  are  shown  corresponding  to  direct 
vision  prism  angular  positions  of  rotation.  This  example  is  representative  of  the  angular  orientation 
with  respect  to  the  optic  axis  simulated  for  all  spectral  wavefront  fields  in  the  software  model. 


spectral  point  spread  function  centre  locations  in  the  focal  plane  array  transverse  plane  resulting 
from  direct  vision  prism  rotation. 

Performing  the  wave  optics  propagation  using  solely  this  process  revealed  favourable  outcomes 
for  prism  angular  orientations  that  situated  the  directions  of  dispersion  tending  along  a  vertical 
orientation  and  unfavourable  outcomes  for  prism  angular  orientations  that  situated  the  directions 
of  dispersion  tending  along  a  horizontal  orientation.  Favourable  outcomes  correspond  to  prism 
angular  orientations  near  0  degrees  and  180  degrees,  while  unfavourable  outcomes  correspond  to 
prism  angular  orientations  near  90  degrees  and  270  degrees. 
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Figure  4.8:  Point  Spread  Functions  for  all  spectral  wavelengths  as  a  Chromotomographic  real 
image  for  Prism  Angular  Position  at  (a)  0  degrees,  and  (b)  92.9032  degrees  from  upright  orientation. 
Locations  for  the  Point  Spread  Functions  in  the  Detector  Focal  Plane  Array  were  determined  using 
Fourier  optics  principles  only. 

The  outcome  obtained  for  the  90  degrees  case  of  angular  prism  orientation  is  shown  in  Fig¬ 
ure  4.8  (b).  The  spectral  point  spread  functions  for  electromagnetic  fields  at  spectral  wavelengths 
to  either  side  of  and  near  the  undeviated  spectral  wavelength  appear  favourable,  while  the  spectral 
point  spread  functions  for  fields  at  spectral  wavelengths  further  from  the  undeviated  spectral  wave¬ 
length  appear  unfavourable.  An  attempt  to  display  the  unfavourable  outcomes  for  both  extremes  of 
the  finite  spectral  range  of  wavelengths  appear  at  Figure  4.9  (a)  for  electromagnetic  fields  at  spec¬ 
tral  wavelengths  between  2.0  pan  and  3.2  pm,  and  Figure  4.9  (b)  for  fields  at  spectral  wavelengths 
between  3.8  pm  and  5.0  pm. 

In  simulating  the  dispersed  spectral  wavefront  fields  of  electromagnetic  radiation  at  the  de¬ 
tector  focussing  lens,  the  vertical  and  horizontal  intervals  between  discrete  samples  of  the  spectral 
wavefront  fields  no  longer  remains  constant.  These  intervals  vary  according  to  (1)  spectral  wave¬ 
length  as  a  result  of  non-uniform  dispersion  in  the  prism,  indicated  along  the  radial  line  of  spectral 
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Figure  4.9:  Magnified  view  of  Point  Spread  Functions  for  spectral  wavelengths  as  a  Chromo- 
tomographic  real  image  for  Prism  Angular  Position  at  92.9032  degrees  from  upright  orientation. 
Locations  for  the  Point  Spread  Functions  in  the  Detector  Focal  Plane  Array  were  determined  us¬ 
ing  Fourier  optics  principles  only  where,  (a)  represents  spectral  wavelengths  between  2.0  /jin  and 
3.2  /jin,  and  (b)  represents  spectral  wavelengths  between  3.8  /im  and  5.0  /jin. 

point  spread  functions  at  the  detector  focal  plane  array,  and  (2)  angular  orientation  of  the  di¬ 
rect  vision  prism,  where  prism  angular  orientations  deviating  from  an  upright  or  vertical  position 
introduce  non-uniform  dispersion  along  the  x-axis  in  system  coordinates.  For  the  upright  prism 
orientation  the  x-axis  was  simulated  with  spatial  sampling  intervals  that  were  uniform  across  spec¬ 
tral  wavefront  fields.  Introducing  non-uniform  variation  in  spatial  sampling  intervals  along  the 
x-axis,  due  to  (1)  simulating  the  angular  rotation  of  the  direct  vision  prism  and  (2)  different  an¬ 
gular  directions  of  propagation  with  respect  to  the  optic  axis  for  spectral  wavefront  fields  at  the 
detector  focussing  lens,  alters  the  accuracy  of  determining  the  spectral  point  spread  function  centre 
locations  using  the  rotational  transform  relationship.  Remember  the  centre  locations  for  spectral 
point  spread  functions  were  obtained  for  an  upright  angular  orientation  of  the  direct  vision  prism 
using  Fourier  optics  propagation  theory  only. 
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Some  approach  was  required  to  locate  the  spectral  point  spread  function  centres  successfully 
for  all  discretely  sampled  prism  angular  orientations.  A  radial  search  function  was  developed  and 
incorporated  in  the  software  model  to  perform  the  task  of  accurately  pin-pointing  the  respective 
spectral  point  spread  function  centre  locations  for  every  prism  angular  orientation.  The  radial 
search  function  employs  Fourier  optics  propagation  theory  with  the  rotational  transform  linear 
algebra  relationship  to  identify  the  vicinity  of  the  respective  spectral  point  spread  function  centre 
locations  in  the  focal  plane  array  transverse  plane,  then  assembles  a  finite  length  radial  line  of 
spatial  sample  intensity  points  aligned  with  the  prism  angular  orientation.  Wave  optics  propagation 
is  then  used  to  obtain  the  intensity  profile  along  this  radial  line.  Searching  along  the  radial  line  of 
spatial  sample  points  for  the  peak  intensity,  and  identifying  the  location  of  the  spatial  sample  point 
with  the  peak  intensity  completes  the  task  of  defining  the  favourable  location  of  all  spectral  point 
spread  function  centres  for  all  prism  angular  orientations.  These  spectral  point  spread  function 
centre  locations  were  used  as  the  centres  for  the  spatial  areas  defined  to  simulate  the  wave  optics 
propagation  of  the  spectral  wavefront  fields  from  the  detector  focussing  lens  to  the  focal  plane  array. 

Incorporating  the  radial  search  function  in  the  software  model  resulted  in  a  dramatic  improve¬ 
ment  with  accurately  simulating  the  spectral  point  spread  functions  across  the  finite  spectral  range 
of  electromagnetic  radiation  at  the  focal  plane  array.  Figure  4.10  displays  the  improvement  in  the 
spectral  point  spread  functions  for  the  90  degree  angular  orientation  of  prism  rotation.  Figure  4.11 
offers  a  comparison  to  assess  the  improvement  made  with  respect  to  the  unfavourable  results  shown 
in  Figure  4.9  for  both  spectral  extremes  of  spectral  point  spread  functions. 

To  appreciate  the  impact  this  anomaly  had  on  the  accuracy  of  simulated  power  for  samples 
in  the  spectral  point  spread  function,  Figure  4.12  shows  the  power  for  each  sample  computed 
over  the  spatial  inaccurately  define  areas  used  for  simulating  the  2.0  /im,  3.0  /jm,  4.0  /.an,  5.0  /an 
spectral  point  spread  functions.  After  incorporating  the  radial  search  method  to  locate  the  centres 
of  the  spectral  point  spread  functions,  the  wave  optics  propagation  resulted  in  significantly  superior 
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Figure  4.10:  Point  Spread  Functions  for  all  spectral  wavelengths  as  a  Chromotomographic  real 
image  for  Prism  Angular  Position  at  (a)  0  degrees,  and  (b)  92.9032  degrees  from  upright  orientation. 
Locations  for  the  Point  Spread  Functions  in  the  Detector  Focal  Plane  Array  were  determined  by 
combining  Fourier  optics  principles  with  the  radial  search  technique. 

outcomes  for  the  simulated  spectral  point  spread  functions.  These  outcomes  for  the  the  2.0  /jin, 
3.0  /jin,  4.0  /jm,  5.0  /jm  spectral  point  spread  functions  are  shown  in  Figure  4.13  with  the  peak 
intensity  of  the  spectral  point  spread  functions  located  in  the  centre  of  each  spatial  area,  and 
coinciding  with  the  off-centre  locations  where  peak  intensity  would  be  expected  in  Figure  4.12. 

The  spectral  point  spread  function  for  the  undeviated  wavelength  remained  unaffected  from 
these  anomalies,  as  evidenced  in  Figure  4.14  showing  the  spectral  point  spread  function  for  a 
90  degrees  angular  orientation  of  the  direct  vision  prism. 

4-5.5  Energy  Distribution  under  the  Point  Spread  Functions.  Look  at  the  shape  of  the 
point  spread  function  for  each  spectral  wavelength  of  electromagnetic  radiation.  The  intensity 
profile  between  spectral  wavelengths  of  electromagnetic  radiation  at  the  focal  plane  array. 
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Figure  4.11:  Magnified  view  of  Point  Spread  Functions  for  spectral  wavelengths  as  a  Chromo- 
tomographic  real  image  for  Prism  Angular  Position  at  92.9032  degrees  from  upright  orientation. 
Locations  for  the  Point  Spread  Functions  in  the  Detector  Focal  Plane  Array  were  determined  by 
combining  Fourier  optics  principles  with  the  radial  search  technique.  Note  the  centres  of  each  Point 
Spread  Function  are  accurately  located,  enabling  the  imaging  area  to  capture  the  entire  spread  of 
spectral  wavelength  intensity. 

The  point  spread  functions  for  the  spectral  fields  of  electromagnetic  radiation  located  at  the 
detector  focal  plane  array  can  be  compared  with  the  point  spread  functions  for  these  spectral  fields 
located  at  the  field  stop.  Viewing  the  spectral  point  spread  functions  at  the  field  stop  in  Figure  4.3, 
and  comparing  these  with  the  point  spread  functions  at  the  focal  plane  array  in  Figure  4.13  for 
spectral  wavelengths  2.0  /iin,  3.0  /im,  4.0  /iin  and  5.0  /iin,  the  peak  intensities  at  the  focal  plane 
array  are  approximately  one  order  of  magnitude  below  the  peak  intensities  at  the  focal  plane  array. 
The  decrease  in  peak  intensity  between  the  field  stop  and  the  focal  plane  array  is  also  coupled 
with  an  increase  in  spatial  coverage  of  the  distribution  of  intensity  for  each  spectral  point  spread 
function,  demonstrating  the  conservation  of  energy  for  each  spectral  point  spread  function.  The 
effects  of  diffraction  on  each  spectral  wavefront  field  of  electromagnetic  radiation  due  to  propagation 
through  the  aft  lens,  direct  vision  prism,  and  focussing  lens  optical  components  are  also  evident  by 
viewing  the  differences  between  these  spectral  point  spread  functions. 
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Figure  4.12:  Point  Spread  Functions  for  2.0  pm,  3.0  pm,  4.0  pm,  5.0  ptm,  with  Prism  Angular 
Position  at  92.9032  degrees  from  upright  orientation.  Locations  for  each  Point  Spread  Function  in 
the  Detector  Focal  Plane  Array  were  determined  using  Fourier  optics  principles  only.  The  centre  of 
each  imaging  area  does  not  accurately  locate  the  centres  for  every  spectral  point  spread  function. 
Further  evaluation  of  the  centre  locations  for  every  spectral  point  spread  function  is  required  from 
the  software  model. 

f.6  Summary  of  Assumptions 


The  development  of  the  software  model  to  simulate  the  sensor  optics  system  of  the  chro- 
motomographic  imaging  system  has  involved  a  number  of  assumptions.  Some  assumptions  have 
been  disclosed  during  discussion  of  the  software  model  while  others  have  not  been  mentioned.  This 
section  is  an  attempt  to  identify  the  assumptions  that  made  a  contribution  to  the  development  of 


the  software  model. 
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Figure  4.13:  Point  Spread  Functions  for  2.0  pm,  3.0  pm,  4.0  pm,  5.0  /iiri,  with  Prism  Angular 
Position  at  92.9032  degrees  from  upright  orientation.  Locations  for  each  Point  Spread  Function  in 
the  Detector  Focal  Plane  Array  were  determined  by  combining  Fourier  optics  principles  with  the 
radial  search  technique.  The  centre  of  each  imaging  area  accurately  locates  the  centres  for  every 
spectral  point  spread  function. 


Lenses  are  modelled  as  thin  lenses  only,  using  the  paraxial  approximation  for  the  surface 
contour  of  each  lens  as  parabolic  in  curvature  instead  of  a  spherical  curvature.  Another  de¬ 
sign  constraint  was  imposed  using  the  thin  lens  representation  regarding  the  F^number  crite¬ 
ria  for  each  lens.  Demonstration  of  the  software  model  presented  in  this  thesis  obeyed  the  con¬ 
straint  F#  number  >  10. 

All  intrinsic  characteristics  of  the  detector  focal  plane  array  such  as  material  sensivities,  non¬ 
uniformities,  responsivity,  are  not  represented  in  the  software  model.  The  simulated  intensity  at 
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Prism  Angular  Position  =  92.9032  degrees,  Spectral  Wavelength  =  3.6pm  x 
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Figure  4.14:  Point  Spread  Function  for  3.6  /jin  with  Prism  Angular  Position  at  92.9032  degrees 
from  upright  orientation.  Location  for  the  imaging  area  centre  Point  Spread  Function  in  the 
Detector  Focal  Plane  Array  was  determined  using  Fourier  optics  principles  only. 


the  focal  plane  array  available  for  detection  with  a  focal  plane  array  has  been  represented  in  the 
software  model. 

Transmission  influences  of  optical  materials  used  for  the  lenses  and  bulk  media  for  the  direct 
vision  prism  were  only  considered  from  an  influencing  phase  perspective.  Attenuation  effects  were 
not  considered. 

Atmospheric  transmission  characteristics  related  to  extinction  of  electromagnetic  radiation 
due  to  atmospheric  absorption  and  scattering  were  not  considered. 
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Field  Intensities  at  the  Detector  Focal  Plane  Array 
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Figure  4.15:  Intensity  profiles  of  the  Point  Spread  Functions  for  all  spectral  wavelength  bands  of 
electromagnetic  radiation  between  2.0  /iin  and  5.0  /.an. 

The  effects  of  a  random  atmosphere  influencing  the  wavefront  of  electromagnetic  radiation 
during  the  initial  stage  of  propagation  was  not  considered. 

The  remote  scene  scenario  was  simulated  without  the  presence  of  background  electromagnetic 
radiation. 

The  wave  optics  propagation  was  simulated  for  a  stationary  direct  vision  prism  for  each  an¬ 
gular  orientation  of  the  prism.  In  collecting  the  spectral  point  spread  functions  through  simulating 
wave  optics  propagation  for  one  entire  rotation  cycle  of  the  prism,  a  stationary  direct  vision  prism 
was  simulated  at  each  different  angular  orientation  of  the  prism.  The  collection  of  wave  optics  prop- 
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agations  through  the  system  for  different  angular  orientations  simulated  the  wave  optics  simulation 
of  the  system  for  a  rotating  prism. 
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V.  Conclusions 


Hyperspectral  imaging  systems  passively  sense  radiant  electromagnetic  energy  from  a  remote 
scene  to  form  a  three  dimension  profile  of  the  remote  scene.  The  data  contained  in  this 
profile  describes  real  images  of  the  remote  scene  (two  spatial  dimensions)  for  a  certain  number  of 
spectral  wavelength  bands  across  a  finite  spectral  range  of  electromagnetic  radiation  (one  spectral 
dimension).  Typical  grating  type  hyperspectral  imaging  systems  collect  spectral  electromagnetic 
radiation  in  the  visible  and  near  infrared  spectral  range,  by  incrementally  scanning  across  the  spatial 
extent  of  the  remote  scene.  Due  to  the  very  low  optical  throughput  used  in  these  systems,  a  legacy 
of  the  optical  techniques  employed  to  scan  for  spectral  and  spatial  information,  adapting  these 
systems  for  spectral  ranges  extending  into  the  mid  wave  and  long  wave  infrared  is  difficult.  Since 
features  in  the  remote  scene  at  infrared  spectral  wavelengths  have  characteristically  low  contrast 
between  radiant  energy  levels  emitted  from  remote  features,  hyperspectral  imaging  in  the  infrared 
demands  higher  optical  throughput  techniques  to  overcome  the  low  contrast  challenges  with  existing 
hyperspectral  imaging  systems. 

The  chromotomographic  imaging  system  was  developed  by  Jonathon  M.  Mooney  and  col¬ 
leagues  of  AFRL/SNHI  at  Hanscom  Air  Force  Base  [4-6],  utilising  an  optical  system  that  angularly 
multiplexes  the  spatial  electromagnetic  radiation  collected  from  a  remote  scene  onto  a  focal  plane 
array  detector,  to  provide  significantly  greater  optical  throughput  for  collecting  radiant  electromag¬ 
netic  energy  from  remote  scenes.  A  rotating  prism  disperses  the  collected  electromagnetic  radiation 
according  to  the  finite  spectral  range  of  wavelength  bands  to  be  used  with  obtaining  spectral  real 
images  as  a  hyperspectral  real  image  data  set  of  the  remote  scene.  The  CTIS  has  been  reported 
to  demonstrate  reasonable  quality  hyperspectral  images  for  finite  spectral  ranges  in  the  mid  wave 
and  long  wave  infrared. 

This  thesis  investigates  the  design  for  a  software  model  that  emulates  the  attributes  of  the 
optics  system  loosely  detailed  in  reports  of  the  CTIS.  Using  optics  principles  and  wave  optics 
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propagation  theory,  the  software  model  has  been  designed  to  simulate  the  point  spread  function 
characteristics  of  the  optics  system  in  the  chromotomographic  real  image  plane.  The  software  model 
can  be  used  as  an  aid  to  develop  and  test  the  characteristics  of  an  optics  system  through  simulation, 
or  simulate  acquiring  actual  chromotomographic  real  image  based  on  spectral  information  about 
a  remote  scene,  or  explore  certain  hyperspectral  image  reconstruction  techniques  suited  to  the 
chromotomograhic  imaging  system. 

5.1  Software  Model 

The  software  model  developed  in  this  research  simulates  the  spectral  point  spread  functions 
and  chromotomographic  point  spread  functions  characterising  the  chromotomographic  sensor  for 
the  chromotomographic  imaging  system.  The  design  of  the  software  model  was  founded  on  the 
optics  principles  related  by  geometric  optics,  Fourier  optics,  and  wave  optics  propagation. 

Geometric  optics  principles  were  used  to  examine  the  propagation  of  the  electromagnetic 
radiation  using  ray  tracing  principles.  These  principles  enabled  an  appreciation  of  the  aperture 
sizes,  focussing  lengths  for  real  imaging,  and  dispersion  angles  through  the  sensor  optics  system  to 
be  evaluated.  Use  of  Snell’s  Law  of  Refraction  and  the  Lensmaker’s  formula  also  aided  with  gaining 
an  appreciation  of  the  system. 

Fourier  optics  principles  were  used  to  guide  the  expectation  of  wavefront  behaviour  at  certain 
apertures  throughout  the  software  model  design,  plus  serve  as  the  foundation  to  establish  the  spatial 
sampling  requirements  for  each  significant  aperture  in  the  sensor  optics  system.  In  establishing  the 
spatial  sampling  requirements,  general  optics  principles  describing:  where  real  images  at  the  remote 
scene  would  be  expected,  at  the  field  stop  and  the  focal  plane  array;  the  differences  between  coherent 
and  incoherent  imaging  systems,  where  the  cut-off  spatial  frequency  for  an  incoherent  imaging 
system  is  double  the  cut-off  frequency  for  a  coherent  imaging  system,  results  in  the  number  of 
spatial  samples  required  to  sample  the  electromagnetic  field  across  the  forward  lens  evaluated  as 
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half  the  number  of  spatial  samples  evident  in  the  remote  scene.  This  outcome  achieved  a  noticeable 
reduction  in  computation  workload  for  wave  optic  propagation  using  the  Rayleigh  Sommerheld 
diffraction  formula. 

Fourier  optics  principles  proved  useful  only  as  a  start  for  obtaining  the  true  relationship 
between  the  detector  focussing  lens  and  the  focal  plane  array.  The  software  model  delivered  a 
wave  optics  simulated  wavefront  with  characteristics  that  deviated  sufficiently  from  Fourier  optics 
principles  to  required  further  processing  of  spatial  data  with  the  radial  search  technique  to  locate 
each  spectral  point  spread  function  accurately. 

Wave  optics  principles  were  used  to  represent  the  simulated  wavefront  propagation  through  all 
stages  of  the  simulated  chromotomographic  sensor,  characterised  by  Rayleigh-Sommerheld  diffrac¬ 
tion  theory.  The  Rayleigh-Sommerheld  diffraction  formula  was  modified  to  a  numeric  expression 
that  aligned  with  the  discrete  sampling  requirements  adopted  for  each  significant  aperture  in  the 
sensor  optics  system. 

The  algorithms  used  in  the  MatLab®  functions  for  the  software  model  are  explained  in  Chap¬ 
ter  III.  A  step-by-step  process  for  using  the  MatLab®  functions  to  develop  the  system  parameters 
for  a  simulated  software  model  of  the  chromotomographic  sensor  has  been  detailed  in  Chapter  IV. 
Results  for  a  simulated  sensor  optics  system  are  also  presented  in  Chapter  IV,  including  the  re¬ 
quirement  to  improve  the  accuracy  of  the  algorithms  used  to  simulate  the  final  wave  optics  stage 
of  propagation  in  the  detector  optics  system. 

The  software  suite  can  be  adapted  to  simulate  a  sensor  optics  system  suited  for  hyperspectral 
imaging  over  a  different  finite  spectral  range  of  electromagnetic  radiation  from  the  remote  scene. 
This  would  only  require  the  use  of  different  optical  material  characteristics,  say  index  of  refraction 
parameters,  and  other  system  parameters  discussed  in  Chapter  IV  to  achieve  a  different  system 
configuration  with  the  same  software  suite. 
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The  software  suite  developed  during  this  research  is  intended  to  be  a  versatile  baseline  for 
developing  simulated  point  spread  functions  results  for  chromotomographic  sensor  designs.  Use  of 
the  software  suite  presented  in  this  thesis  used  system  parameters  suited  for  the  MWIR  spectral 
range  to  simulate  the  design  configuration  of  the  chromotomographic  sensor.  The  versatility  of 
the  software  suite  enables  different  optical  material  characteristics,  say  index  of  refraction  param¬ 
eters,  and  other  system  parameters  discussed  in  Chapter  IV  to  be  used  to  tailor  the  simulated 
chromotomographic  sensor  design  for  different  spectral  ranges  of  electromagnetic  radiation  from 
the  MWIR. 

The  baseline  reached  during  this  research  for  the  software  suite  also  enables  the  software  model 
design  to  be  adapted  for  further  investigation  into  design  attributes  of  the  chromotomographic 
sensor,  or  possibly  allow  the  option  for  progressively  more  complex  additions  to  be  incorporated  in 
further  revisions  of  the  software  model  with  further  research. 

5.2  Research  Motivation 

Motivation  for  improving  upon  the  MASINT  understanding  of  the  radiant  energy  phenomenol¬ 
ogy  in  the  fireball  resulting  from  an  explosion  or  explosions  in  the  remote  scene  is  the  foundation  of 
interest  behind  this  research  of  the  chromotomographic  sensor.  To  capture  the  spectral  and  spatial 
details  from  transient  events  in  the  remote  scene,  such  as  explosions,  in  the  context  of  hyperspec- 
tral  imagery,  a  sensor  with  the  optical  throughput  characteristics  of  the  chromotomographic  sensor 
is  needed.  Some  questions  of  that  foster  the  motivation  behind  this  research  are  included  in  the 
following  paragraph. 

Is  the  size  of  the  fireball,  as  a  time  resolved  hyperspectral  real  image  data  set,  have  a  spectral 
wavelength  dependence  ?  What  is  the  spectral  and  spatial  profile  of  the  detonation  wavefront  from 
the  explosion  creating  the  fireball  ?  Other  research  efforts  at  AFIT  reveal  from  observations  of 
explosions  in  the  visible  spectrum  that  radiant  energy  decay  rates  do  not  follow  a  Planckian  trend. 
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What  does  the  infrared  spectrum  reveal  about  the  radiant  energy  decay  rates  from  similar  types 
of  explosions  ?  Can  the  radial  temperature  gradient  profile  of  the  fireball  be  obtained  from  the 
spectral  two  dimension  real  images  of  the  remote  scene  constructed  from  the  radiant  energy  data 
collected  at  the  detector  focal  plane  array  of  the  chromotomographic  sensor  ? 

The  chromotomographic  sensor’s  high  optical  throughput  feature  offers  the  ability  to  collect 
the  spectral  electromagnetic  radiation  content  from  the  remote  scene  to  form  a  hyperspectral  real 
image  data  set  representing  the  remote  scene  in  a  significantly  lower  duration  of  time  (seconds) 
compared  with  other  grating  type  conventional  scanning  hyperspectral  imaging  systems  (hours). 
This  lower  time  duration  with  the  staring  chromotomographic  sensor  allows  for  speedier  sequential 
collection  over  time  of  radiant  energy  from  the  remote  scene  to  obtain  a  time  resolved  series  of 
hyperspectral  real  image  data  sets  representing  the  remote  scene.  Achieving  time-resolved  hyper¬ 
spectral  imagery  of  transient  events  elevates  the  dimensionality  of  for  this  type  of  imagery  with  the 
inclusion  of  time  as  a  dimension  to  a  four  dimension  real  image  data  set. 

The  software  model  simulates  the  unit  impulse  response  of  the  chromotomographic  sensor 
as  spectral  point  spread  functions  and  chromotomographic  point  spread  functions.  In  the  context 
of  modelling  the  radiant  energy  from  explosions  in  the  remote  scene,  collected  by  the  chromoto¬ 
mographic  sensor,  the  software  model  may  be  used  to  simulate  an  explosion  (input)  with  unity 
magnitude  across  the  spectral  profile  for  detonation,  the  initial  time  period,  of  the  explosion  in  the 
remote  scene.  This  context  of  modelling  would  also  rely  on  the  radiant  energy  from  the  explosion 
being  the  significantly  dominant  source  of  radiant  energy  in  the  remote  scene. 

With  its  significantly  greater  optical  throughput  feature  the  chromotomographic  sensor  ap¬ 
pears  well  suited  for  collecting  the  spectral  content  of  radiant  energy  from  explosions.  Furthermore 
the  software  suite  developed  during  this  research  is  a  useable  software  tool  to  simulate  and  evaluate 
different  configurations  for  a  chromotomographic  sensor  design. 
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5.3  Analysis  of  system  design  with  software  model 


Analysis  of  chromotomographic  sensor  system  designs  simulated  with  the  software  model  is 
readily  possible.  A  case  in  point  is  the  intensity  profile  of  spectral  point  spread  functions  between 
2.0  /jm  and  5.0  /im  shown  in  Figure  4.15,  also  shown  at  Figure  5.1  (a).  An  expectation  of  spectral 
point  spread  function  spatial  profiles  is  for  the  peak  intensity  to  follow  an  increasing  trend  as  the 
spectral  wavelength  of  the  point  spread  functions  follows  a  decreasing  trend.  The  conservation  of 
energy  collected  at  the  focal  plane  array  is  the  basis  of  this  relationship,  since  the  total  energy  (or 
total  intensity)  available  for  every  spectral  point  spread  function  at  the  focal  plane  array  should  be 
the  same  (or  very  close). 

The  peak  intensity  relationship  evident  from  the  spatial  profile  slice  of  the  simulated  chromo¬ 
tomographic  point  spread  function  is  the  trend  of  increasing  peak  intensity  does  occur  for  spectral 
point  spread  functions  above  3.6  /mi.  This  spectral  wavelength  also  coincided  with  the  undeviated 
wavelength  for  the  simulated  chromotomographic  sensor  design.  For  spectral  wavelengths  below 
3.6  /mi  the  trend  revealed  from  the  software  model  shows  the  trend  of  peak  intensity  rolling  off  and 
tending  to  decrease  with  decreasing  spectral  point  spread  functions.  This  change  in  trend  for  the 
peak  intensity  of  spectral  point  spread  functions  counters  the  intuitive  relationship  based  on  the 
conservation  of  energy.  A  possible  conclusion  assessed  from  this  information  may  link  this  change 
in  trend  for  the  peak  intensity  to  the  undeviated  wavelength  of  the  simulated  chromotomographic 
sensor. 

Further  information  available  from  the  software  model,  specifically  the  total  intensities  for 
every  spectral  wavefront  of  electromagnetic  radiation  as  spectral  point  spread  functions  at  the  field 
stop  and  the  focal  plane  array,  allowed  a  comparison  of  the  total  intensity  throughput  between  the 
field  stop  and  the  focal  plane  array  based  on  the  spectral  content  of  electromagnetic  radiation  to 
be  made.  Figure  5.2  reveals  the  results  of  this  comparison,  where  Figure  5.2  (a)  shows  the  total 
spectral  intensity  at  the  field  stop  and  focal  plane  array,  and  Figure  5.2  (b)  shows  the  percentage 
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Figure  5.1:  Spectral  Point  Spread  Function  spatial  profiles  at  the  focal  plane  array  between  spec¬ 
tral  wavelengths  2.0  pm  and  5.0  pm  for  different  simulated  direct  vision  prism  configurations,  where 
(a)  represents  3.6  pm  undeviated  spectral  wavelength,  prism  entry  surface  angle  of  30°  00'  00", 
prism  exit  surface  angle  of  23°  56'  00",  (b)  represents  2.0  /zm  undeviated  spectral  wavelength, 
prism  entry  surface  angle  of  30°  00'  00",  prism  exit  surface  angle  of  24°  58'  47",  (c)  represents 
3.6  /jm  undeviated  spectral  wavelength,  prism  entry  surface  angle  of  15°  00'  00",  prism  exit  sur¬ 
face  angle  of  11°  44'  32",  and  (d)  represents  2.0  /im  undeviated  spectral  wavelength,  prism  entry 
surface  angle  of  15°  00'  00",  prism  exit  surface  angle  of  12°  17'  41". 


loss  in  total  spectral  intensity  between  the  field  stop  and  focal  plane  array.  Clearly  a  trend  of 
declining  percentage  in  total  intensity  throughput  between  the  field  stop  and  the  focal  plane  array 
as  the  spectral  point  spread  functions  follow  a  trend  of  decreasing  spectral  wavelength.  Furthermore 
this  decrease  gets  more  dramatic,  losing  from  22%  to  58%  of  total  intensity  in  the  spectral  point 
spread  functions  as  the  spectral  wavelength  decreases  from  5.0  pm  to  2.0  pm  respectively. 

Total  spectral  intensities  are  more  revealing  than  following  the  trend  of  peak  intensity  for 
spectral  point  spread  functions  because  the  total  spectral  intensities  revealed  a  definite  trend  of 
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(a)  (b) 

Figure  5.2:  Total  spectral  intensities  between  spectral  wavelengths  2.0  /jm  and  5.0  /jm  at  the  field 
stop  and  the  focal  plane  array  (a),  and  (b)  loss  in  total  spectral  intensity  throughput  between  the 
field  stop  and  the  focal  plane  array  for  spectral  point  spread  functions  between  spectral  wavelengths 
2.0  /jin  and  5.0  pin 

spectral  intensity  throughput  loss  between  the  field  stop  and  the  focal  plane  array.  Understanding 
whether  design  aspects  of  the  simulated  chromotomographic  sensor  influence  the  degree  of  total 
spectral  intensity  throughput  loss  forms  the  motivation  for  the  following  notes  with  investigating 
to  a  limited  extent  certain  attributes  of  the  sensor  design  between  the  field  stop  and  the  focal  plane 
array. 

Investigating  the  phenomenon  of  total  spectral  intensity  throughput  loss  further,  subtle 
changes  to  the  simulated  chromotomographic  sensor  design  were  made  by  varying  the  direct  vision 
prism  entry  surface  angle  and  exit  surface  angle  systematically,  and  note  the  differences  in  the 
chromotomographic  point  spread  function,  spectral  point  spread  functions  and  the  total  intensity 
available  at  the  focal  plane  array.  Figure  5.1  (b)  through  (d)  are  the  chromotomographic  point 
spread  functions  obtained  for  different  simulated  chromotomographic  sensor  design  configurations. 
The  simulated  sensor  design  represented  at  Figure  5.1  (b)  has  the  same  direct  vision  prism  aper¬ 
ture  size  and  same  prism  entry  surface  angle  of  30°  from  Figure  5.1  (a),  however  a  different  prism 
exit  surface  angle  to  accommodate  an  undeviated  spectral  wavelength  of  2.0  /.« n.  The  simulated 
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sensor  design  represented  at  Figure  5.1  (c)  has  the  same  direct  vision  prism  aperture  size  and  same 
undeviated  spectral  wavelength  of  3.6  /iin  from  Figure  5.1  (a),  however  the  prism  entry  surface 
angle  was  changed  to  15°  and  the  prism  exit  surface  angle  changed  accordingly  to  accommodate 
the  undeviatecl  spectral  wavelength.  The  simulated  sensor  design  represented  at  Figure  5.1  (d)  has 
the  same  direct  vision  prism  aperture  size  and  same  prism  entry  surface  angle  from  Figure  5.1  (c), 
however  a  different  prism  exit  surface  angle  to  accommodate  an  undeviated  spectral  wavelength  of 
2.0  /jm. 

Analysis  of  the  total  intensities  for  every  spectral  wavefront  of  electromagnetic  radiation  as 
spectral  point  spread  functions  at  the  field  stop  and  the  focal  plane  array  revealed  no  discernible 
differences  in  total  intensity  values  and  percentage  loss  in  total  intensity  throughput  between  the 
field  stop  and  the  focal  plane  array  trends  shown  at  Figure  5.2.  This  outcome  suggests  that  varying 
the  geometries  of  the  direct  vision  prism  entry  and  exit  surfaces  has  no  discernible  influence  on  the 
total  intensity  throughput  between  the  field  stop  and  the  focal  plane  array,  however  does  influence 
the  spatial  profile  of  the  spectral  point  spread  functions  and  the  chromotomographic  point  spread 
function. 

The  purpose  for  exploring  the  different  configurations  of  the  direct  vision  prism  related  to  in¬ 
vestigating  whether  these  sensor  design  changes  influenced  the  proportion  of  total  spectral  intensity 
throughput  loss  observed  from  analysis  of  the  simulated  chromotomographic  sensor.  The  spectral 
wavefront  fields  simulated  using  wave  optics  propagation  theory  have  some  measure  of  wavefront 
error  and  establishing  system  design  choices  that  influence  the  wavefront  error  may  assist  with 
some  level  of  mitigating  the  loss  in  spectral  total  intensity  throughput  between  the  field  stop  and 
the  focal  plane  array. 

The  aft  lens  aperture  was  the  next  system  parameter  evaluated  as  a  candidate  for  potentially 
influencing  the  total  spectral  intensity  throughput  loss  in  the  simulated  chromotomographic  sensor 
design.  Two  sensor  design  configurations  were  assessed  through  simulation,  the  only  distinguishable 
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(a)  (b) 

Figure  5.3:  Spectral  Point  Spread  Function  spatial  intensity  profiles  for,  (a)  the  22  mm  diame¬ 
ter,  and  (b)  the  30  mm  diameter,  aft  lens  aperture  sensor  design  configurations  between  spectral 
wavelengths  2.0  /r m  and  5.0  /mi 


difference  between  these  configurations  being  the  diameter  for  the  aft  lens  aperture.  Figure  5.3 
shows  the  simulated  spatial  profile  of  the  chromotomographic  point  spread  functions  for  each  sensor 
design  configuration  where  Figure  5.3  (a)  simulates  the  sensor  design  with  an  aft  lens  diameter  of 
22  mm  and  direct  vision  prism  diameter  of  30  mm,  and  Figure  5.3  (b)  simulates  the  sensor  design 
with  an  aft  lens  diameter  of  30  mm  and  direct  vision  prism  diameter  of  38  mm.  The  values  used  for 
the  aft  lens  diameter  were  arbitrarily  chosen,  however  the  direct  vision  prism  diameter  values  used 
were  evaluated  using  the  software  suite  to  ensure  a  sensible  software  model  was  used  to  simulate 
the  sensor  design  configurations. 

Figure  5.3  (b)  clearly  reveals  an  increase  in  the  peak  intensities  for  every  spectral  point  spread 
function  and  an  increase  in  the  peak  intensities  across  the  spatial  profile  of  the  chromotomographic 
point  spread  function  at  the  focal  plane  array.  Viewing  the  total  spectral  intensities  for  the  spectral 
point  spread  functions  at  the  field  stop  and  the  focal  plane  array  respectively  for  the  30  mm  aft 
lens  diameter  sensor  design  configuration  at  Figure  5.4  (a),  and  the  previous  chart  at  Figure  5.2  (a) 
for  the  22  mm  aft  lens  diameter  sensor  design  configuration,  these  charts  show  a  similar  trend  for 
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(a)  (b) 

Figure  5.4:  Total  spectral  intensities  between  spectral  wavelengths  2.0  /rm  and  5.0  /iin  at  the 
field  stop  and  the  focal  plane  array  for,  (a)  the  30  mm  diameter  aft  lens  aperture  sensor  design 
configuration,  and  (b)  comparison  of  loss  in  total  spectral  intensity  throughput  between  the  field 
stop  and  the  focal  plane  array  for  spectral  point  spread  functions  between  spectral  wavelengths 
2.0  /jin  and  5.0  /iin  for  the  22  mm  diameter  and  30  mm  diameter  aft  lens  aperture  sensor  design 
configurations 


total  spectral  intensities.  Comparing  the  percentage  loss  in  total  spectral  intensity  throughput  from 
the  field  stop  to  the  focal  plane  array  for  both  sensor  design  configurations,  shown  in  the  chart 
at  Figure  5.4  (b),  a  noticeable  reduction  in  total  spectral  intensity  throughput  loss  is  apparent 
for  increasing  spectral  wavelengths.  However,  for  spectral  point  spread  functions  tending  toward 
the  minimum  spectral  wavelength  there  is  no  appreciable  difference  in  total  spectral  intensity 
throughput  between  the  simulated  sensor  design  configurations.  This  outcome  is  counter  intuitive 
to  the  expected  results  from  a  physical  sensor  design  with  similar  changes  in  design  configuration. 

Thorough  analysis  of  the  simulated  sensor  design  configurations  represented  with  the  software 
model  was  not  accomplished.  Items  worth  investigating  in  the  software  model  include  the  treatment 
of  wavefront  error  in  the  wave  optics  propagation  for  the  prism  optics  system  and  detector  optics 
system,  plus  the  spatial  sampling  criteria  with  respect  to  the  aperture  sizes  presented  in  the  various 
simulated  sensor  design  configurations.  The  sampling  criteria  used  in  the  software  model  seems 
to  sufficiently  sample  all  spectral  wavefront  fields  nominated  for  simulation  purposes,  however 
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considering  the  spatial  variation  of  the  point  spread  functions  for  lower  spectral  wavelengths  may 
dictate  greater  spatial  sampling  is  required  to  improve  the  sampled  spatial  profile  for  these  spectral 
wavelengths. 

Aside  from  extending  or  investigating  the  software  model  design  more  thoroughly,  there  are 
other  areas  of  further  research  related  to  the  research  effort  reported  in  this  thesis  worthy  of 
consideration  and  introduced  in  the  following  paragraphs. 

5-4  Further  Research 

The  software  model  developed  in  this  research  simulates  the  point  spread  functions  character¬ 
ising  the  sensor  optics  system  for  the  chromotomographic  imaging  system.  The  design  maturity  of 
the  software  model  would  enable  an  optics  technician  or  engineer  to  simulate  a  design  configuration 
for  the  sensor  optics  system,  and  evaluate  the  spectral  real  imaging  and  chromotomographic  real 
imaging  aspects  of  this  system  design.  Further  research  may  be  undertaken  using  this  baseline  of 
the  software  model  to  explore  other  aspects  with  simulating  the  design  configuration  for  a  sensor 
optic  system,  or  evaluating  digital  post-processing  techniques  for  the  chromotomographic  imaging 
system.  Some  areas  of  further  research  are  mentioned  and  may  contain  far  greater  scope  than  the 
few  lines  allocated  in  the  following  paragraphs. 

5.4.1  Chromotomographic  Real  Images.  Using  the  spectral  radiant  energy  characteristics 
for  an  arbitrary  remote  scene,  spectral  real  images  in  the  chromotomographic  real  image  plane  can 
be  obtained  using  the  software  model. 

The  spatial  details  of  the  remote  scene  for  each  spectral  wavelength  of  incident  electromagnetic 
radiation  represents  a  different  system  input  to  the  sensor  optics  system.  The  point  spread  function 
for  each  spectral  wavelength  of  electromagnetic  radiation  represents  a  different  system  impulse 
response  characterising  the  sensor  optics  system.  Using  linear  system  theory,  a  two  dimension 
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convolution  between  each  system  input  and  each  system  impulse  response  may  be  used  to  obtain 
the  system  output  for  the  respective  system  input. 

Using  the  software  model  to  simulate  the  system  impulse  response  for  each  spectral  wavelength 
of  electromagnetic  radiation  at  each  angular  position  of  rotation  for  the  direct  vision  prism,  a 
simulated  system  output  can  be  obtained.  The  system  outputs  being  the  real  image  of  the  remote 
scene  for  each  spectral  wavelength  of  electromagnetic  radiation  at  the  detector  focal  plane  array 
transverse  plane.  The  collection  of  spectral  real  images  for  each  angular  position  of  rotation  for  the 
direct  vision  prism  represents  the  chromotomographic  real  image  of  the  remote  scene. 

The  remote  scene  of  interest  is  approximated  as  a  two  dimension  spatial  image  even  though 
the  scene  is  actually  three-dimensional  event  over  time.  Point  spread  functions  characterising  the 
simulated  sensor  optics  system  for  different  ranges  from  the  forward  lens  can  be  developed  with 
the  software  model  (system  impulse  response)  and  convolved  with  the  associated  features  from  the 
remote  scene  (system  inputs)  to  obtain  a  simulated  representation  of  the  chromotomographic  real 
image  of  the  remote  scene  (system  output) . 

5-4-2  Physical  System  Parameters.  The  software  model  developed  with  this  research  does 
not  have  many  of  the  practical  influences  evident  with  a  physical  sensor  optics  system  included  in 
the  simulated  system.  Transmission  characteristics  of  optical  components,  attenuation  parameters 
characterising  the  extinction  of  radiant  energy  in  the  propagating  spectral  wavefronts  due  to  at¬ 
mospheric  absorption  and  scattering,  the  presence  of  background  radiant  energy  surrounding  the 
remote  scene  and  evident  from  the  optical  components  used  for  the  sensor  optics  system,  collec¬ 
tively  contribute  as  some  of  the  added  influences  evident  in  a  physical  system.  Other  features  of 
physical  components  in  the  sensor  optics  system  that  may  warrant  further  research  include  the 
optical  properties  of  the  bulk  media  used  in  the  direct  vision  prism  varying  due  to  temperatures, 
pressures,  and  operating  environment.  Similarly,  influences  with  focal  plane  array  non-uniformities 
in  detecting  the  radiant  intensity  may  be  worth  further  study. 
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Developing  versions  of  the  software  model  that  contain  these  influences  would  offer  a  closer 
simulation  of  the  physical  sensor  optics  system. 

5-4-3  Maximum  Acceptance  Angle.  The  size  of  the  field  stop  aperture  will  influence  the 
maximum  acceptance  angle  subtended  at  the  forward  lens,  between  the  optic  axis  and  the  location 
of  a  source  of  radiant  electromagnetic  energy  in  the  remote  scene.  The  size  of  the  field  stop 
aperture,  the  range  distance  of  the  remote  scene  from  the  forward  lens,  are  aspects  with  the  sensor 
optics  system  that  can  be  evaluated  with  the  software  model  to  simulate  the  collection  of  radiant 
electromagnetic  energy  for  different  angular  locations  of  remote  scene  sources.  An  evaluation  of  the 
maximum  acceptance  angle  for  the  sensor  optics  system,  based  on  the  scenario  of  range  and  spatial 
extent  used  to  represent  the  remote  scene.  Mean  square  error  of  the  intensity  profiles  between 
point  spread  functions  obtained  for  the  simulated  system  for  point  sources  located  on  the  optic  axis 
compared  with  point  sources  located  some  angular  location  from  the  optic  axis,  may  be  used  to 
evaluate  the  maximum  acceptance  angle. 

5-4-4  Radiometry.  The  software  model  represents  the  wavefront  propagation  for  spectral 
fields  of  electromagnetic  radiation  as  point  spread  functions.  The  radiant  energy  used  for  these 
spectral  fields  is  a  1  watt  source  emitting  spherical  wavefronts  of  electromagnetic  radiation  spatially 
across  a  complete  spherical  surface.  Whatever  fraction  of  this  radiant  energy  incident  on  the  forward 
lens  is  represented  as  an  electromagnetic  field  in  the  software  model.  Photon  energy  emission  is 
the  source  of  this  radiant  energy,  so  characterising  the  radiometric  aspects  of  this  radiant  energy 
for  an  arbitrary  remote  scene  would  provide  a  connection  with  the  photon  phenomena  occurring  in 
the  remote  scene. 

The  spectral  fields  of  electromagnetic  radiation  collected  by  the  sensor  optics  system  originate 
from  the  radiant  energy  emissions  from  features  in  the  remote  scene.  This  radiant  energy  occurs 
from  photon  emissions  in  the  remote  scene,  where  each  photon  emission  is  a  quantised  amount  of 
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energy  according  to  spectral  wavelength  v  characterising  the  photon  energy  emission,  as  described 


by  equation  5.1. 


E  =  hi/  =  h(-^)  (5.1) 

From  the  spectral  real  images  of  the  remote  scene,  the  phenomena  of  photon  emissions  may 
be  derived  from  the  fields  of  spectral  electromagnetic  radiation  represented  in  the  software  model. 
Photons  with  low  oscillation  have  low  energy  and  high  spectral  wavelength,  whereas  photons  with 
high  oscillation  have  high  energy  and  low  spectral  wavelength. 

5-4-5  Digital  Post-Processing  Techniques.  A  further  research  area  extending  from  de¬ 
veloping  chromotomographic  real  images  mentioned  in  section  5.4.1  would  be  exploring  different 
digital  post-processing  techniques  like  the  pseudo  inverse  method  [5,6]  for  assembling  a  hyperspec- 
tral  data  set  of  spectral  real  images  of  the  remote  scene  from  the  collection  of  chromotomographic 
real  images  available  at  the  detector  focal  plane  array. 
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Appendix  A.  Alternate  Interpretations  of  the  Sensor  Optics  System 


A.l  Linear  System  Interpretation 

Aperture  (forward  lens) 


Figure  A.l:  Linear  Systems  Interpretation  of  Sensor  Optics  System 

The  following  stages  of  propagation.  The  symbols  used,  O  represents  the  spectral  wavelength 
electromagnetic  held  as  spatial  information  for  objects  or  subjects  of  interest  observed  in  the  remote 
scene,  B  represents  the  background  spectral  wavelength  electromagnetic  held  propagating  from  the 
surrounding  environment  for  the  stage  of  propagation,  r  represents  the  extinction  function  for  the 
stage  of  propagation,  r  represents  the  position  reference  as  a  radial  distance  from  the  optic  axis  for 
locations  in  the  spectral  wavelength  electromagnetic  held  for  stages  of  propagation,  A  represents 
the  spectral  wavelengths  under  examination  with  the  sensor  optics  system,  g  represents  the  output 
electromagnetic  held  for  a  stage  of  propagation  and  the  input  electromagnetic  held  for  the  following 
stage  of  propagation. 

Implementing  a  simple  version  of  the  model  with  these  expressions  can  be  done  with  consider¬ 
ing  the  background  radiation  parameter  B  set  to  zero,  thereby  removing  any  background  radiation 
effects  in  the  model,  and  the  extinction  function  parameter  r  set  to  unity,  thereby  operating  the 
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system  in  an  atmospheric  environment  with  no  magnitude  attenuation  of  electromagnetic  radiation 
from  propagating  through  the  atmosphere.  Results  of  this  simple  model  are  presented  for  different 
transverse  planes  along  the  optic  axis  that  represent  different  stages  of  the  radiation  (broadband 
signal)  propagation  through  the  forward  optics  system.  The  general  notation  for  each  distinct 
transfer  function  between  stages  of  propagation,  H12,  subscript  (1)  represents  the  commencing  sys¬ 
tem  transverse  plane  and  subscript  (2)  represents  the  concluding  system  transverse  plane  through 
the  distinct  stage  of  propagation  characterised  by  the  transfer  function  H12. 

A.  1.1  Remote  Scene. 

B0(r0,  A)  +  0(ro,  A) 

A.  1.2  Forward  Optics  (after  Forward  Lens). 

ga(ra,  A)  =  H0a(r0  — »  ra,  A)  {  r0(r0,  A)  [  B0(r0,  A)  +  0(ro,  A)  ]  } 

A. 1.3  Forward  Optics  (Field  Stop). 

gb(rb,  A)  =  Hab(ra  — ►  rb,  A)  {  ra(ra,  A)  [  Ba(ra,  A)  +  g(ra,  A)  ]  } 

A.  1.4  Forward  Optics  (after  Aft  Lens). 

gc(rc,  A)  =  Hbc(rb  — ►  rc,  A)  {  rb(rb,  A)  [  Bb(rb,  A)  +  g(rb,  A)  ]  } 

A.  1.5  Prism  Optics  (Entry  Surface). 

gd(rd,  A)  =  Hcd(rd  — >  rc,  A)  {  rc(rc,  A)  [  Bc(rc,  A)  +  g(rc,  A)  ]  } 

A.  1.6  Prism  Optics  (Centre  Surface). 

ge(re,  A)  =  Hde(rd  — »  re,  A)  {  rd(rd,  A)  [  Bd(rd,  A)  +  g(rd,  A)  ]  } 
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A.  1.7  Prism  Optics  (Exit  Surface). 


gf Of,  A)  =  Hef(re  — »  rf,  A)  {  re(re,  A)  [  Be(re,  A)  +  g(re,  A)  ]  } 
A.  1.8  Detector  Optics  (after  Focussing  Lens). 

gkOk,  A)  =  Hfk(rf  — ►  rk,  A)  {  rf  (rf,  A)  [  Bf  (rf ,  A)  +  g(rf,  A)  ]  } 
A. 1.9  Detector  Optics  (Focal  Plane  Array). 

gs(rs,  A)  =  Hks(rk  — »  rs,  A)  {  rk(rk,  A)  [  Bk(rk,  A)  +  g(rk,  A)  ]  } 
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A. 2  Geometric  Optics  Interpretation 


Aperture  (forward  lens) 


Figure  A. 2:  Geometric  Optics  Interpretation  of  Sensor  Optics  System 


Table  A.l:  Matrices  for  Forward  Section  of  the  Sensor  Optics  System 


Matrix 

Expression 

Description  of  Propagation 

Mi 

1  Zi 
0  1 

1st  stage  of  propagation 

m2 

- 

1  O' 

~k  \ 

2nd  stage  of  propagation 

m3 

1  z2 
0  1 

3rd  stage  of  propagation 

m4 

1  z3 

0  1 

4th  stage  of  propagation 

m5 

- 

1  O' 

-k 

5th  stage  of  propagation 
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Table  A. 2:  Matrices  for  Prism  Section  of  the  Sensor  Optics  System 


Table  A. 3:  Matrices  for  Aft  Section  of  the  Sensor  Optics  System 


Matrix 

Expression 

Description  of  Propagation 

M13 

'  1  O' 

rr3  ij 

13th  stage  of  propagation 

M14 

1  z8 

0  1 

14th  stage  of  propagation 
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Appendix  B.  Indexes  of  Refraction  for  Optical  Materials  in  the  Infrared 


Indexes  of  refraction  for  optical  materials  that  have  relatively  high  transmission  properties  for  the 
infrared  spectrum  of  electromagnetic  radiation  have  been  collected.  The  properties  of  optical 
materials  according  to  their  respective  indexes  of  refraction  can  be  modelled  using  Matlab®. 

The  tables  below  contain  the  spectral  indexes  of  refraction  for  Lithium  Fluoride  (LiF),  Barium 
Fluoride  (BaF2),  Sodium  Chloride  (NaCl),  Calcium  Fluoride  (CaF2),  and  Sapphire  (AI2O3),  as 
optical  materials  of  interest. 


Table  B.l:  Spectral  Indexes  of  Refraction  for  Lithium  Fluoride  (LiF) 


index  of  refraction 

wavelength  (A)  in  /nn 

1.38711 

1 

1.37875 

2 

1.36660 

3 

1.34942 

4 

1.32661 

5 

Table  B.2:  Spectral  Indexes  of  Refraction  for  Barium  Fluoride  (BaF2) 


index  of  refraction 

wavelength  (A)  in  /jm 

1.46856 

1 

1.46460 

2 

1.46115 

3 

1.45102 

5 

1.43566 

7 

Table  B.3:  Spectral  Indexes  of  Refraction  for  Sodium  Chloride  (NaCl) 


index  of  refraction 

wavelength  (A)  in  /un 

1.5303 

1.1786 

1.5231 

3.5359 

1.51356 

6.4825 
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Table  B.4:  Spectral  Indexes  of  Refraction  for  Calcium  Fluoride  (CaF2) 


index  of  refraction 

wavelength  (A)  in  /jm 

1.42888 

1 

1.42385 

2 

1.41785 

3 

1.40963 

4 

1.39895 

5 

1.38559 

6 

Table  B.5:  Spectral  Indexes  of  Refraction  for  Sapphire  (AI2O3) 


index  of  refraction 

wavelength  (A)  in  /jm 

1.74618 

1.550 

1.73748 

2.010 

1.73232 

2.249 

1.719 

2.703 

1.712 

2.941 

1.701 

3.333 

1.687 

3.704 

1.674 

4.000 

1.658 

4.348 

1.636 

4.762 

1.623 

5.000 

1.607 

5.263 
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Appendix  C.  Software  Suite 


C.l  Simulated  Chromotomographic  Sensor 

C.1.1  Sensor  Optics  System  Mat-Lab®  Function. 


7,  - 

7.  CTHIS_type_llc_prism_02.m 

7.  - 

"/,  -  computes  the  propagation  of  the  electric  field  (denoted  U) , 

'/,  starting  at  the  transverse  plane  where  the  electric  field 

'/,  is  denoted  U_start ,  and  ending  at  the  transverse  plane  where 

7.  the  electric  field  is  denoted  U_end. 

7.  - 

"/,  EENG799  (Independent  Study)  Hyperspectral  Chromotomography 

7.  Date:  12  February  2004 

7.  FLTLT  Anthony  J.  Dearinger 

'/,  Royal  Australian  Air  Force 

•/,  GE-04M 

7.  - 


’/.  - 

7.  START 
7.  - 


•/.  - 

clear  all 

7.  - 

prm  =  parameters; 

•/.  - 

•/.  - 

sprintf (’ Computation  for  propagation  to  the  FORWARD  APERTURE  PLANE  has  commenced’) 

7.  - 

[U_2a  prm]  =  propagate_01_fwd_aperture (prm) ; 

I  - 

sprintf (’ Computation  for  propagation  to  the  FORWARD  APERTURE  PLANE  has  completed’) 

7.  - 


I  - 

sprintf (’ Computation  for  the  FORWARD  APERTURE  has  commenced’) 

7.  - 

[U_2b  prm]  =  fwd_aperture(prm,  U_2a) ; 

7.  - 

sprintf (’ Computation  for  the  FORWARD  APERTURE  has  completed’) 

7.  - 

7.  - 

sprintf (’ Computation  for  propagation  to  the  FIELD  STOP  PLANE  has  commenced’) 

7.  - 

[U_3a  prm]  =  propagate_02_f ield_stop (prm,  U_2b) ; 

U_3b  =  U_3a; 

7.  - 

sprintf (’ Computation  for  propagation  to  the  FIELD  STOP  PLANE  has  completed’) 

7.  - 


7.  7. 


7.  sprintf (’ Computation  for  the  FIELD  STOP  APERTURE  has  commenced’) 

7.  7.  - 

7.  [U_3b  prm]  =  f ield_stop_aperture(prm,  U_3a)  ; 

7.  7.  - 


7.  sprintf (’ Computation  for  the  FIELD  STOP  APERTURE  has  completed’) 

7.  7.  - 
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•/.  - 

sprintf ( ’Computation  for  propagation  to  the  AFT  APERTURE  PLANE  has  commenced’) 

•/.  - 

[U_4a  prm]  =  propagate_03_aft_aperture (prm,  U_3b) ; 

'/.  - 

sprintf ( ’Computation  for  propagation  to  the  AFT  APERTURE  PLANE  has  completed’) 

•/.  - 

'/.  - 

sprintf ( ’Computation  for  the  AFT  APERTURE  has  commenced’) 

•/.  - 

[U_4b  prm]  =  aft_aperture(prm,  U_4a) ; 

'/.  - 

sprintf ( ’Computation  for  the  AFT  APERTURE  has  completed’) 

•/.  - 


•/.  - 

save  data_test_10a_prm  prm 
clear  ans  prm 

save  data_test_10a_f ield_optics 

'/.  - 

clear  all 

load  data_test_10a_prm 
load  data_test_10a_f ield_optics 
•/.  - 


'/.  - 

sprintf ( ’Computation  for  propagation  through  the  DIRECT  VISION  PRISM  has  commenced’) 

•/.  - 

[U_5  prm  prism]  =  propagate_04_prism(prm,  U_4b) ; 

•/.  - 

sprintf ( ’Computation  for  propagation  through  the  DIRECT  VISION  PRISM  has  completed’) 

7.  - 


7.  - 

save  data_test_10a_prm  prm 

save  data_test_10a_f ield_prism  U_5 

save  data_test_10a_f ield_prism_prm  prism 

7.  - 

clear  all 

load  data_test_10a_prm 

load  data_test_10a_f ield_prism 

7.  - 


7.  - 

sprintf ( ’Computation  for  propagation  to  the  INFRA-RED  DETECTOR  has  commenced’) 

7.  - 

U_6_search_l  =  []  ; 

U_6_search_2  =  [] ; 

U_6  =  []  ; 

7.  - 

[U_6_search_l ,  prm]  =  propagate_05_IR_array_01_search(prm,  U_5) ; 

7.  - 

7.  [U_6_search_2,  prm]  =  propagate_05_IR_array_02_search(prm,  U_5)  ; 

7.  - 

[U_6,  prm]  =  propagate_05_IR_array_03_image (prm,  U_5) ; 

7.  - 

sprintf ( ’Computation  for  propagation  to  the  INFRA-RED  DETECTOR  has  completed’) 

7.  - 


7.  - 

save  data_test_10a_prm  prm 

save  data_test_10a_f ield_detector  U_6 

7.  - 
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clear  all 

load  data_test_10a_prm 

load  data_test_10a_f ield_detector 

•/.  - 


•/.  - 

sprintf ( [’Computation  for  assembling  the  Point  Spread  Functions  ’ 
’on  the  INFRA-RED  DETECTOR  ARRAY  has  commenced’]) 

7,  - 

prm  =  assemble_PSF_on_IR_array (prm,  U_6) ; 

•/.  - 

sprintf ( [’Computation  for  assembling  the  Point  Spread  Functions  ’ 
’on  the  INFRA-RED  DETECTOR  ARRAY  has  completed’]) 

I  - 


•/.  - 

save  data_test_10a_prm  prm 

•/.  - 


V.  — 

•/.  END 

•/. 
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C.2  Wave  Optics  Propagation  MatLab®  Functions 

C.2.1  Propagate  from  Remote  Scene  to  Forward  Lens  Aperture. 


function [U_end,  prm]  =  propagate_01_f wd_aperture (prm) ; 

'/.  - 

7,  propagate_01_fwd_aperture 

'/.  - 

%  -  computes  the  propagation  of  the  electric  field  from  the  OBJECT  PLANE  (ie.  INPUT  SCENE) 

7,  to  the  FORWARD  APERTURE  PLANE. 

7.  -  computes  the  propagation  of  the  electric  field  (denoted  U) ,  starting  at  the  transverse  plane 
7,  where  the  electric  field  is  denoted  U_start,  and  ending  at  the  transverse  plane  where  the 

7.  electric  field  is  denoted  U_end. 

7,  -  this  function  employs  the  Rayleigh-Sommerf ield  diffraction  formula  to  describe  the 

7,  propagation  of  a  diverging  EM  radiation  wavefront 

7.  - 

7,  EENG799  (Independent  Study)  Hyperspectral  Chromo tomography 
7.  Date:  17  November  2003 
7,  FLTLT  Anthony  J.  Dearinger 
7,  Royal  Australian  Air  Force 
7.  GE-04M 

7.  - 


7.  - 

7.  START 

7.  - 

sprintf  (  ’ - ’ ) 

sprintf ( ’Propagation  from  the  POINT  SOURCE  to  the  FORWARD  APERTURE  PLANE’) 

sprintf  (  ’ - ’ ) 

sprintf ( ’Computation  has  commenced  for  all  wavelengths’) 

7.  - 

prop_dist  =  []  ; 

tempX  =  prm .pt_source_x_loc  -  prm. system. x2_index; 
tempY  =  prm.pt_source_y_loc  -  prm. system. y2_index; 
t  empZ  =  prm . z 1 _r ange ; 

prop_dist  =  sqrt (tempX. ~2  +  tempY. ~2  +  tempZ.~2); 

7.  - 


7.  - 

radiant_intensity  =  prm. amplitude  /  (4  *  pi);  7«  [watts  per  steradian] 

7.  - 

7,  solid  angle  -  use  small  angle  approximation  for  computing  the  solid  angle 
7,  solid  angle  -  use  rectangular  aperture  to  help  with  computing  the  distribution 
7,  of  incident  radiant  flux  across  the  aperture 

7.  - 

7.  solid  angle  -  circular  aperture 

7.  solid_angle  =  (pi  /  4)  *  (prm.dia_aperture_fwd~2  /  prm. zl_range~2)  ;  7«  [steradians] 

7.  - 

7,  solid  angle  -  rectangular  aperture 

solid_angle  =  prm. dia_aperture_fwd~2  /  prm.zl_range~2;  l  [steradians] 

7.  - 

incident_f lux  =  radi ant _ intensity  *  solid_angle;  7#  [watts] 

7.  - 

7.  sampled_incident_f lux  -  use  small  angle  approximation,  therefore  assume  flux  is 
7,  uniformly  distributed  across  the  forward  aperture 

7.  - 

sampled_incident_f lux  =  incident_f lux  /  (prm.pts  *  prm.pts)  *  ... 

prm. system. dxl  *  prm. system . dyl ; 

7.  - 

c  =  299792458;  7*  [metre  /  sec]  ,  speed  of  light 

permit ivity_free_space  =  8.85e-12;  7«  [sec~2  Coulombs'^  /  (metres  -  kilogram)] 

7.  - 

sampled_f ield_amplitude  =  sqrt (sampled_incident_f lux  /  c  /  permitivity_f ree_space) ; 

•/#  - 
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•/.  - 

for  si  =  1  :  prm. lambda_bands 

U_end(  :  ,  :  ,  si)  =  sampled_f ield_amplitude ( : , si)  .*  exp(j  *  prm.k(sl)  .*  prop_dist) 

end 

•/.  - 

sprintf (’ Computation  has  completed  for  all  wavelengths') 

sprintf  ( ’ - ' ) 

•/.  - 


•/. 

return 
•/.  — - 
'/,  END 
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C.2.2  Propagate  from  Forward  Lens  Aperture  to  Field  Stop  Aperture. 

function [U_end,  prm]  =  propagate_02_f ield_stop(prm,  U_start) ; 

7.  - 

7.  propagate_02_f  ield_stop  .m 

7.  - 

7*  -  computes  the  propagation  of  the  electric  field  from  the  FORWARD  APERTURE  PLANE 

7.  to  the  FIELD  STOP  PLANE. 

7.  -  computes  the  propagation  of  the  electric  field  (denoted  U) ,  starting  at  the  transverse  plane 
7.  where  the  electric  field  is  denoted  U_start,  and  ending  at  the  transverse  plane  where  the 

7,  electric  field  is  denoted  U_end. 

7.  - 

7,  EENG799  (Independent  Study)  Hyperspectral  Chromo tomography 
7,  Date:  11  February  2004 
7.  FLTLT  Anthony  J.  Dearinger 
7«  Royal  Australian  Air  Force 
7.  GE-04M 

7.  - 


7.  - 

7.  START 
7.  - 


7.  - 

sprintf  (  5 - ’  ) 

sprintf ( ’Propagation  from  the  FORWARD  APERTURE  PLANE  to  the  FIELD  STOP  PLANE’) 

sprintf  (  ’ - ’  ) 

sprintf ( ’Computation  has  commenced  for  all  wavelengths’) 

7.  - 


7.  - 

clear  temp* 

7.  - 


7.  - 

lens_thick_fwd  =  [] ; 

7.  - 

tempX  =  prm. system. x2_index; 
tempY  =  prm. system. y2_index; 

7.  - 

lens_thick_fwd  =  prm.delta_0_fwd  -  (tempX. ~2  +  tempY. ~2  )  *  ... 

(  (prm.Rl_fwd)  /  (2*(prm.Rl_fwd) ~2)  -  (prm.R2_fwd)  /  (2* (prm.R2_fwd) "2)  ); 
7.  - 

prm. lens_thick_fwd  =  lens_thick_f wd; 

7.  - 


7.  - 

7.  expression  A2  -  lens  thickness  function 

7.  - 

7,  apply  the  lens  thickness  function  to  the  phase  of  the  electric  field  of  the  optical 
7,  electromagnetic  radiation 

7.  - 

for  si  =  1  :  length (prm. lambda) 

U_start ( : , : , si)  =  U.start ( : , : , si)  .*  ... 

exp(j  *  prm.k(sl)  .*  (lens_thick_fwd) * (prm.n_fwd  -  1)); 

end 

7.  - 


7.  - 

U_start  =  reshape (U_start,  prm.pts,  (prm.pts)*(prm.lambda_bands)) ; 

7.  - 

clear  temp* 

7.  - 

temp_lambda  =  kr on (prm. lambda,  ones (prm.pts) ) ; 
temp_K  =  kron(prm.k,  ones (prm.pts) ) ; 


120 


•/.  - 

7,  - 

for  sX  =  1  :  prm.pts 

•/.  - 

for  sY  =  1  :  prm.pts 

7.  - 

tempX  =  prm. system. x2_index  -  prm. system. x3_index(sY, sX) ; 
tempY  =  prm. system. y2_index  -  prm. system. y3_index(sY, sX) ; 
tempZ  =  (prm. z2_fwd_f ocal) * (ones (prm. pts) ) ; 

7.  - 

prop_dist  =  []  ; 

prop_dist  =  sqrt(  tempX. ~2  +  tempY. ~2  +  tempZ. ~2  ); 

7.  - 

7.  - 

7,  for  si  =  1  :  length  (prm.  lambda) 

7,  U_end(sY,sX,sl)  =  sum(  sum(  squeeze (U_start (:,:, si) )  .*  ... 

7,  (exp(j  *  prm.k(sl)  *  dist_phase))  )  ); 

7,  end 

7.  - 

7.  comments 

7,  -  U_start  must  not  undergo  the  reshape  at  line  056  for  this  loop  to  work 

7.  - 

7.  - 

7,  expression  A1  -  approximated  Rayleigh-Sommerf ield  diffraction  propagation 
7,  expression  A1  -  was  originally  used  with  lens  thickness  function 
7,  (expression  A2) 

7,  - 

7,  temp_U  =  U_start  ./  repmat (prop_dist ,  1,  prm. lambda_bands)  .*  ... 

7,  exp(  j  *  kron(prm.k,  ones  (prm. pts) )  .*  repmat (prop_dist ,  1,  prm.  lambda_bands)  ); 

7.  - 

7.  - 

7,  expression  B1  -  approximated  Rayleigh-Sommerf ield  diffraction  propagation 
7,  expression  B1  -  more  accurate  numerical  approximation  compared  to 
7,  propagation  using  expression  A1 

7.  - 

temp_U_amplitude  =  sinc(  repmat ((prm. system. x3_index(sY,sX)  *  pi  *  prm. system. dx2) , 

prm.pts,  prm.pts  *  prm. lambda_bands)  ./  repmat (tempZ ,  1,  prm. lambda_bands)  ./  ... 

temp_lambda  )  . *  ... 

sinc(  repmat ( (prm. system. y3_index(sY, sX)  *  pi  *  prm. system. dy2) ,  ... 

prm.pts,  prm.pts  *  prm. lambda_bands)  ./  repmat (tempZ ,  1,  prm. lambda_bands)  ./  ... 

temp_lambda  )  . *  ... 

prm. system. dy2  .*  prm. system. dy2  .*  repmat (tempZ,  1,  prm. lambda_bands)  ./  ... 

j  ./  temp_lambda; 

7.  - 

temp_U  =  U_start  .*  temp_U_amplitude  ./  ... 

repmat  (prop_dist .  "'2,  1,  prm.  lambda_bands)  .*  ... 

exp(  j  *  temp_K  .*  repmat (prop_dist ,  1,  prm. lambda_bands)  ); 

z  - 

temp_U  =  sum(temp_U) ; 

temp_U  =  sum(  reshape (temp_U,  prm.pts,  prm . lambda_bands)  ); 

7o  - 

U_end(sY, sX, : )  =  temp_U; 

7.  - 

clear  temp_U* 

7o  - 

end 

7.  - 

sprintf(JsX  =  7»dJ  ,  sX) 

7.  - 

end 

y#  - 
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sprintf (’ Computation  has  completed  for  all  wavelengths') 

sprintf  ( ’ - ’ ) 

i  - 

l  - 

return 

•/.  - 

•/.  END 

7,  - 
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C.2.3  Propagate  from  Field  Stop  Aperture  to  Aft  Lens  Aperture. 

function [U_end,  prm]  =  propagate_03_af t_aperture (prm,  U_start) ; 

7.  - 

7.  propagate_03_aft_aperture 

7.  - 

7*  -  computes  the  propagation  of  the  electric  field  from  the  FIELD  STOP  PLANE 

7.  to  the  AFT  APERTURE  PLANE. 

7.  -  computes  the  propagation  of  the  electric  field  (denoted  U) ,  starting  at  the  transverse  plane 
7.  where  the  electric  field  is  denoted  U_start,  and  ending  at  the  transverse  plane  where  the 

7,  electric  field  is  denoted  U_end. 

7.  - 

7,  EENG799  (Independent  Study)  Hyperspectral  Chromo tomography 
7,  Date:  11  February  2004 
7.  FLTLT  Anthony  J.  Dearinger 
7,  Royal  Australian  Air  Force 
7.  GE-04M 

7.  - 


7.  - 

7.  START 
7.  - 


7.  - 

sprintf  (  ’ - ’ ) 

sprintf ( ’Propagation  from  the  FIELD  STOP  PLANE  to  the  AFT  APERTURE  PLANE’) 

sprintf  (  ’ - ’ ) 

sprintf ( ’Computation  has  commenced  for  all  wavelengths’) 

7.  - 


7.  - 

clear  temp* 

7.  - 


7.  - 

U_start  =  reshape (U_start,  prm.pts,  (prm. pts) * (prm. lambda_bands) ) ; 

7.  - 

temp_lambda  =  kron (prm. lambda,  ones (prm.pts) ) ; 
temp_K  =  kron(prm.k,  ones (prm. pt s) ) ; 

7.  - 


7.  - 

for  sX  =  1  :  prm.pts 

7.  - 

for  sY  =  1  :  prm.pts 

7.  - 

tempX  =  prm. system. x3_index  -  prm. system. x4_index(sY, sX) ; 
tempY  =  prm. system. y3_index  -  prm. system. y4_index(sY, sX) ; 
tempZ  =  (prm. z3_af t_f ocal) * (ones (prm. pts) ) ; 

7.  - 

prop_dist  =  []  ; 

prop_dist  =  sqrt(  tempX. ~2  +  tempY. "2  +  tempZ. ~2  ); 

7*  - 

/  - 

7,  for  si  =  1  :  length  (prm.  lambda) 

7,  U_end(sY, sX, si)  =  sum(  sum(  squeeze (U_start (:,:, si) )  .*  ... 

7,  (exp(j  *  prm.k(sl)  *  dist_phase))  )  ); 

7,  end 

7.  - 

7,  comments 

7,  -  U_start  must  not  undergo  the  reshape  at  line  034  for  this  loop  to  work 

7o  - 

7.  - 
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7,  expression  A1  -  approximated  Rayleigh-Sommerf ield  diffraction  propagation 
/  expression  A1  -  was  originally  used  with  lens  thickness  function 
7,  (expression  A2) 

7.  - 

7.  temp_U  =  U_start  .  /  repmat  (prop_dist ,  1,  prm.  lambda_bands)  .*  ... 

7.  exp(  j  *  kron(prm.k,  ones  (prm.  pts))  .*  repmat  (prop_dist ,  1,  prm.  lambda_bands)  ); 

7.  - 


7.  - 

7,  expression  B1  -  approximated  Rayleigh-Sommerf ield  diffraction  propagation 
7,  expression  B1  -  more  accurate  numerical  approximation  compared  to 
7.  propagation  using  expression  A1 

7o  - 

temp_U_amplitude  =  sinc(  repmat ( (prm. system. x4_index(sY, sX)  *  pi  *  prm. system. dx3) , 
prm. pts,  prm. pts  *  prm. lambda_bands)  ./  repmat (tempZ ,  1,  prm. lambda_bands)  ./ 
temp_lambda  )  . *  ... 

sinc(  repmat ( (prm. system. y4_index(sY,sX)  *  pi  *  prm. system. dy3) ,  ... 

prm. pts,  prm. pts  *  prm. lambda_bands)  ./  repmat (tempZ ,  1,  prm. lambda_bands)  ./ 

temp_lambda  )  . *  ... 

prm. system. dx3  .*  prm. system. dy3  .*  repmat (tempZ ,  1,  prm. lambda_bands)  ./  ... 

j  ./  temp_lambda; 


temp_U  =  U_start  . *  temp_U_amplitude  ./  ... 

repmat (prop_dist . ~2,  1,  prm. lambda_bands)  .  *  ... 

exp(  j  *  temp_K  .*  repmat (prop_dist ,  1,  prm. lambda_bands)  ); 

7.  - 


7.  - 

temp_U  =  sum(temp_U) ; 

temp_U  =  sum(  reshape (temp_U,  prm. pts,  prm . lambda_bands)  ); 

7.  - 

U_end(sY, sX, : )  =  temp_U; 

7.  - 

clear  temp_U* 

7.  - 

end 

7.  - 

sprintf(JsX  =  7odJ  ,  sX) 

7.  - 

end 

7.  - 


7.  - 

lens_thick_aft  =  [] ; 

7.  - 

tempX  =  prm. system. x4_index; 
tempY  =  prm. system. y4_index; 

7.  - 

lens_thick_aft  =  prm.delta_0_aft  -  (tempX. "2  +  tempY. ~2  )  *  ... 

(  (prm.Rl_aft)/(2*(prm.Rl_aft) "2)  -  (prm.R2_aft)/(2*(prm.R2_aft) ~2)  ); 

7.  - 

prm. lens_thick_aft  =  lens_thick_af t ; 

1  - 


7.  - 

7,  expression  A2  -  lens  thickness  function 

7.  - 

7,  apply  the  lens  thickness  function  to  the  phase  of  the  electric  field  of  the  optical 
7.  electromagnetic  radiation 

7.  - 

for  si  =  1  :  length (prm. lambda) 

U_end( : , : , si)  =  U_end( : , : , si)  .*  ... 

exp(j  *  prm.k(sl)  .*  (lens_thick_af t) * (prm.n_aft  -  1)); 

end 

•/#  - 
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clear  temp*; 

•/.  - 


•/.  - 

sprintf (’ Computation  has  completed  for  all  wavelengths') 

sprintf  ( ’ - ' ) 

•/.  - 


•/. 

return 

I.  — - 

•/,  END 

•/. 
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C.2.4  Propagate  through  the  Direct  Vision  Prism  to  Detector  Focussing  Lens. 


function [U_end,  prm,  prism]  =  propagate_04_prism(prm,  U_start) ; 

7.  - 

7,  propagate_04_prism.m 

7.  - 

%  EENG799  (Independent  Study)  Hyperspectral  Chromo tomography 
7*  Date:  05  January  2004 
7.  FLTLT  Anthony  J.  Dearinger 
7,  Royal  Australian  Air  Force 
7,  GE-04M 

7.  - 


7.  - 

7.  START 
7.  - 


7.  - 

7,  initialise  the  parameters  for  the  direct  vision  prism 

7.  - 

prm  =  prism_01_initialise (prm) ; 

7.  - 


7.  - 

7,  define  the  parameters  for  the  direct  vision  prism 
7.  - 

prm  =  prism_02_parameters(prm) ; 

7.  - 


7.  - 

7,  define  the  index  of  refraction  for  the  centre  wavelength  of  each 
7,  wavelength  band  at  each  stage  of  propagation  through  the 
7,  direct  vision  prism 

7.  - 

prm  =  prism_03_ref ractive_index(prm) ; 

y  - 


y#  - 

7,  compute  the  angles  of  dispersion  for  the  centre  wavelength  of  each 
7,  wavelength  band  at  each  stage  of  propagation  through  the 
7.  direct  vision  prism 

7.  - 

prm  =  prism_04_angles (prm) ; 

7.  - 


y#  - 

7,  rotation  of  axes 

7.  -  for  CW  prism  rotation,  axes  must  rotate  CCW 

7,  -  for  CCW  prism  rotation,  axes  must  rotate  CW 

y#  - 

7,  -  outputs  are  horizontal  distances,  and  vertical  distances  according  to 

7,  the  rotation  of  each  sample  point  of  the  field  with  respect  to  the  prism 

7,  orientation 

y  - 

prm  =  prism_05_rotate_x4y4(prm) ; 

7.  - 


y#  - 

7,  theoretical  angles  of  dispersion 

y#  - 

7,  -  compute  the  theoretical  angles  of  dispersion  for  the  centre  wavelength  of 

7,  each  wavelength  band  due  to  propagating  the  electromagnetic  field  through 

7.  the  direct  vision  prism 

7,  -  the  theoretical  angles  computed  are  the  x-axis  and  y-axis  components  of  the 

7,  angles  of  dispersion  with  respect  to  the  transverse  plane  for  the  sensor 
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7,  coordinate  system 

7.  -  the  theoretical  angles  are  also  computed  for  each  position  of  rotation  for 

7,  the  direct  vision  prism,  for  CCW  and  CW  prism  rotation 

7,  -  the  theoretical  angles  will  serve  to  aid  with  determining  each  location  in 

7,  the  detector  focal  plane  array  of  the  imaging  centre  for  the  point  spread 

7o  function  associated  with  the  centre  wavelength  of  each  wavelength  band 

7.  - 

prm  =  prism_06_tilt_angles (prm) ; 

7.  - 


7.  - 

7,  propagate  all  points  in  transverse  plane  through  the  1ST  stage  of  propagation 

7.  - 

7,  -  propagate  from  aft  aperture  to  prism  1ST  surface  interface 

7.  - 

7,  -  outputs  are  propagation  distance,  and  vertical  displacement 

7.  - 


7.  - 

yl_total  =  reshape (prm. prism. y4_CCW,  ... 

(prm.pts  *  prm.pts  *  prm. prism. number_angles) ,  1); 

7.  - 

zl_offset  =  yl_total  .*  tan(prm. prism. alpha_l) ; 

y, - 

zl_total  =  prm. prism. zl_ref  -  zl_offset;  7*  [metres] 
temp_r_l  =  zl_total; 

7.  - 

7.  —  propagation  distance  through  the  1ST  stage  of  propagation,  r_l 

7.  - 

r_l  =  reshape (temp_r_l ,  prm.pts,  prm.pts,  prm.prism.number_angles) ; 

7.  - 

7,  comments 

7,  -  r_l  ==>  columns  read  from  TOP  =  r_l(l,:,:),  to  BOTTOM  =  r_l  (prm  .pts ,  :  ,  : ) 

7,  for  the  transverse  plane  to  the  optic  axis 

7.  - 

7.  - 

7,  refraction  of  all  points  in  transverse  plane  at  the  prism  1ST  surface  interface 

7.  - 

7,  propagate  all  points  in  transverse  plane  through  the  2ND  stage  of  propagation 

7.  - 

7.  -  outputs  are  propagation  distance,  and  vertical  displacement 

7.  - 


7.  - 

7.  —  refraction  at  the  1ST  surface  interface 

7.  - 

temp_r_2  =  prm. prism. thick_surf_l  ./  ... 

repmat (cos (prm. prism. beta_lb) ,  (prm.pts  *  prm.pts  *  prm. prism. number_angles) ,  1) ; 

7.  - 

7,  —  propagation  distance  through  the  2ND  stage  of  propagation,  r_2 

7.  - 

r_2  =  reshape (temp_r_2,  prm.pts,  prm.pts,  prm. prism. number_angles ,  ... 

prm. lambda_bands)  ;  7«  [metres] 

7.  - 

7,  —  vertical  offset  displacement  at  the  2ND  interface,  y2_total 

7,  due  to  refraction  and  dispersion  through  the  media  of  the 

7,  2ND  stage  of  propagation 

7.  - 

y2_total  =  temp_r_2  . *  ... 

repmat (sin (prm. prism. alpha_l  +  prm. prism. beta_lb) ,  ... 

(prm.pts  *  prm.pts  *  prm.prism.number_angles) ,  1)  +  ... 

repmat  (yl_total ,  1,  prm. lambda_bands)  ;  7«  [metres] 

7.  - 
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•/.  - 

7,  refraction  of  all  points  in  transverse  plane  at  the  prism  2ND  surface  interface 

'/.  - 

%  propagate  all  points  in  transverse  plane  through  the  3RD  stage  of  propagation 

7.  - 

7.  -  outputs  are  propagation  distance,  and  vertical  displacement 

7.  - 


7.  - 

7,  —  refraction  at  the  2ND  surface  interface 

7.  - 

z2_ref  =  (prm. prism. height/2)  *  t an (prm. prism. alpha_l)  +  prm. prism. z2_min;  7*  [metres] 
z2_offset  =  y2_total  .*  tan(prm. prism. alpha_l) ;  %  [metres] 

z2_total  =  z2_ref  +  z2_offset;  7*  [metres] 

'/.  - 

temp_r_3  =  z2_total  ./  repmat (cos (prm. prism. alpha_l  +  prm. prism. beta_2b) ,  ... 

(prm.pts  *  prm.pts  *  prm. prism. number_angles) ,  1);  l  [metres] 

'/.  - 

7,  —  propagation  distance  through  the  3RD  stage  of  propagation,  r_3 

'/.  - 

r_3  =  reshape (temp_r_3 ,  prm.pts,  prm.pts,  prm. prism. number_angles ,  ... 

prm. lambda_bands)  ;  7«  [metres] 

'/.  - 

7.  —  vertical  offset  displacement  at  the  3RD  interface,  y3_total 

7,  due  to  refraction  and  dispersion  through  the  media  of  the 

7,  3RD  stage  of  propagation 

7.  - 

y3_total  =  temp_r_3  .*  repmat (sin (prm. prism. alpha_l  +  prm. prism. beta_2b) ,  ... 

(prm.pts  *  prm.pts  *  prm. prism. number_angles) ,  1)  +  y2_total;  l  [metres] 

7.  - 

7.  comments 

7,  -  r_3  ==>  columns  read  from  TOP  =  r_3(l,:,:),  to  BOTTOM  =  r_3(prm.pts,  :  ,  : ) 

7,  for  the  transverse  plane  to  the  optic  axis 

7.  -  y3_total  ==>  columns  read  from  TOP  =  y3_total  (1 ,  :  )  ,  to  BOTTOM  =  y3_total  (prm  .pts ,  :  ) 

7,  for  the  transverse  plane  to  the  optic  axis 

7,  - 


7.  - 

7,  refraction  of  all  points  in  transverse  plane  at  the  prism  3RD  surface  interface 

•/.  - 

7.  propagate  all  points  in  transverse  plane  through  the  4TH  stage  of  propagation 

7.  - 

7,  -  outputs  are  propagation  distance,  and  vertical  displacement 

'/.  - 


'/.  - 

7.  —  refraction  at  the  3RD  surface  interface 

•/.  - 

z3_ref  =  (prm. prism. height  /  2)  *  tan (prm. prism. alpha_3)  +  prm. prism. z3_min;  %  [metres] 
r4_ref  =  cos (prm. prism. epsilon_3b)  . \  z3_ref ;  %  [metres] 

'/.  - 

y4_ref  =  repmat(  r4_ref  . *  sin (prm. prism. beta_3b) ,  ... 

(prm.pts  *  prm.pts  *  prm. prism. number_angles) ,  1)  +  y3_total;  l  [metres] 

•/.  - 

r4_offset  =  sin (prm. prism. alpha_3)  . *  (  y4_ref  ./  ... 

repmat (sin(pi  /  2  -  prm. prism. alpha_3  +  prm. prism. beta_3b) ,  ... 

(prm.pts  *  prm.pts  *  prm. prism. number_angles) ,  1)  ); 

'/.  - 

temp_r_4  =  repmat (r4_ref ,  (prm.pts  *  prm.pts  *  prm. prism. number_angles) ,  1)  -  r4_offset; 

'/.  - 

7,  —  propagation  distance  through  the  4TH  stage  of  propagation,  r_4 

'/.  - 

r_4  =  reshape (temp_r_4,  prm.pts,  prm.pts,  prm. prism. number_angles ,  ... 

prm. lambda_bands)  ;  7«  [metres] 

'/.  - 
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7,  —  vertical  offset  displacement  at  the  4TH  interface,  y4_total 

7.  due  to  refraction  and  dispersion  through  the  media  of  the 

7,  4TH  stage  of  propagation 

7.  - 

y4_total  =  temp_r_4  .*  repmat (sin (prm. prism. bet a_3b) ,  ... 

(prm.pts  *  prm.pts  *  prm. prism. number_angles ) ,  1)  +  y3_total;  %  [metres] 

7.  - 

7,  comments 

7.  -  y4_ref  ==>  columns  read  from  TOP  =  y4_ref(l,:),  to  BOTTOM  =  y4_ref  (prm.pts ,: ) 

7.  for  the  transverse  plane  to  the  optic  axis 

7,  ~  r4_offset  ==>  columns  read  from  TOP  =  r4_of  f  set  (1 ,  :  )  ,  to  BOTTOM  =  r4_off  set  (prm.pts ,:  ) 

7,  for  the  transverse  plane  to  the  optic  axis 

7,  -  r_4  ==>  columns  read  from  TOP  =  r_4(l,:),  to  BOTTOM  =  r_4(prm.pts ,  : ) 

7,  for  the  transverse  plane  to  the  optic  axis 

7,  -  y4_total  ==>  columns  read  from  TOP  =  y4_total  (1 ,  :  )  ,  to  BOTTOM  =  y4_total  (prm  .pts ,  :  ) 

7,  for  the  transverse  plane  to  the  optic  axis 

7.  - 

7,  comments 

7,  -  propagation  through  the  prism  is  treated  as  a  collimated  beam  propagation. 

7.  -  this  means  all  distances  computed,  with  respect  to  the  refractive  index  of  each  media  of 

7,  propagation,  result  is  a  phase  translation  for  each  point  in  the  cross-section  of  the 

7,  beam  of  propagation. 

7.  -  distance  also  affects  the  magnitude  of  the  radiation  for  each  point  in  the  cross-section 

7,  of  the  beam  of  propagation,  since  distance  has  an  inverse  relation  with  magnitude. 

7.  - 


7.  - 

7,  refraction  of  all  points  in  transverse  plane  at  the  prism  4TH  surface  interface 

7.  - 

7,  propagate  all  points  in  transverse  plane  through  the  5TH  stage  of  propagation 

7.  - 

7,  -  outputs  are  propagation  distance,  and  vertical  displacement 

7  - 


7.  - 

7,  —  refraction  at  the  4TH  surface  interface 

7.  - 

z5_offset  =  y4_total  .  *  tan  (prm.  prism.  alpha_3)  ;  °i  [metres] 

7.  - 

z5_total  =  prm.  prism.  z5_ref  +  z5_offset;  7*  [metres] 

7.  - 

temp_r_5  =  z5_total  ./  repmat (cos (prm. prism. alpha_3  +  prm. prism. beta_4b) ,  ... 

(prm.pts  *  prm.pts  *  prm. prism. number_angles) ,  1);  l  [metres] 

7.  - 

7,  —  propagation  distance  through  the  5TH  stage  of  propagation,  r_5 

7,  to  the  focusing  lens  of  the  IR  Detector 

7.  - 

r_5  =  reshape (temp_r_5 ,  prm.pts,  prm.pts,  prm. prism. number_angles ,  ... 

prm. lambda_bands)  ;  7»  [metres] 

7.  - 

7,  —  vertical  offset  displacement  at  the  focussing  lens  of  the  IR  detector, 

7,  y4_total  due  to  refraction  and  dispersion  through  the  media 

7,  (exiting  the  prism)  of  the  5TH  stage  of  propagation 

7.  - 

y5_total  =  temp_r_5  .*  repmat (sin (prm. prism. alpha_3  +  prm. prism. beta_4b) ,  ... 
(prm.pts  *  prm.pts  *  prm. prism. number_angles) ,  1)  +  y4_total;  l  [metres] 

7.  - 


7, 


prism. r. 

_1  =  permute (r_l, 

a 

2 

4 

3]) 

prism. r. 

_2  =  permute (r_2, 

a 

2 

4 

3]) 

prism. r. 

_3  =  permute (r_3. 

[i 

2 

4 

3]) 

prism. r. 

_4  =  permute (r_4, 

a 

2 

4 

3]) 

prism. r. 

_5  =  permute (r_5, 

a 

2 

4 

3]) 
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/o 

temp  =  []  ; 

temp  =  repmat (yl_total ,  1,  prm. lambda_bands) ; 
prism. yl_total  =  permute (temp ,  [1  2  4  3]); 

7.  [metres] 

/  0 

temp  =  []  ; 

temp  =  reshape (y2_total 

,  prm.pts, 

prm.pts,  prm. prism. number_angles, 

prm. lambda_bands) ; 

7»  [metres] 

prism. y2_total  =  permute (temp ,  [1 

2 

4 

3]); 

7.  [metres] 

/  0 

temp  =  []  ; 

temp  =  reshape (y3_total 

,  prm.pts, 

prm.pts,  prm. prism. number_angles, 

prm. lambda_bands) ; 

7»  [metres] 

prism. y3_total  =  permute (temp ,  [1 

2 

4 

3]); 

7,  [metres] 

/  0 

temp  =  []  ; 

temp  =  reshape (y4_total 

,  prm.pts, 

prm.pts,  prm. prism. number_angles, 

prm. lambda_bands) ; 

7,  [metres] 

prism. y4_total  =  permute (temp ,  [1 

2 

4 

3]); 

7,  [metres] 

/  0 

temp  =  []  ; 

temp  =  reshape (y5_total 

,  prm.pts, 

prm.pts,  prm. prism. number_angles, 

prm. lambda_bands) ; 

7»  [metres] 

prism. y5_total  =  permute (temp ,  [1 

2 

4 

3]); 

7,  [metres] 

/o 

temp  =  []  ; 

temp  =  reshape (prm. prism. x4_CCW, 

(prm.pts  *  prm.pts 

*  prm. prism. number_angles) ,  1);  /  [metres] 

temp  =  repmat (temp,  1, 

prm . lambda_bands) ; 

prism. x_total  =  permute (temp ,  [1 

2 

4 

3]); 

•/.  - 

7,  PART  FIVE 

7.  - 

7,  PHASE  effects  upon  the  IR  electromagnetic  radiation 

7.  - 

7,  -  apply  the  phase  effects  due  to  propagation  through  the  direct  vision  prism 

•/.  - 


•/.  - 

clear  tempU*; 

'/.  - 

tempU  =  U_ start; 

'/.  - 

tempU  =  reshape (tempU,  (prm.pts  *  prm.pts),  prm . lambda_bands) ; 
tempU  =  repmat (tempU,  prm. prism. number _ angles ,  1); 

•/.  - 

tempUl  =  tempU  .*  exp(  j  *  kron((prm.k  .*  prm.prism.nl),  ... 

ones( (prm.pts  *  prm.pts  *  prm. prism. number_angles) ,  1))  . *  ... 

repmat (temp_r_l ,  1,  prm. lambda_bands)  ); 

•/.  - 

tempU2  =  tempUl  . *  exp(  j  *  kron((prm.k  .*  prm. prism. n2) ,  ... 

ones ( (prm.pts  *  prm.pts  *  prm. prism. number _ angles) ,  1))  .*  temp_r_2  ); 

7.  - 

tempU3  =  tempU2  .*  exp(  j  *  kron((prm.k  .*  prm. prism. n3) ,  ... 

ones ( (prm.pts  *  prm.pts  *  prm. prism. number_angles) ,  1))  .*  temp_r_3  ); 

7.  - 

tempU4  =  tempU3  .*  exp(  j  *  kron((prm.k  .*  prm. prism. n4) ,  ... 

ones ( (prm.pts  *  prm.pts  *  prm. prism. number_angles) ,  1))  .*  temp_r_4  ); 

7.  - 

tempU5  =  tempU4  .*  exp(  j  *  kron((prm.k  .*  prm. prism. n5) ,  ... 

ones ( (prm.pts  *  prm.pts  *  prm. prism. number_angles) ,  1))  .*  temp_r_5  ); 

7.  - 
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I  - 

prism. U_5a  =  reshape (tempUl ,  prm.pts,  prm.pts,  prm. prism. number _angles , 
prm. lambda_bands) ; 

prism. U_5a  =  permute (pr ism. U_5a,  [1  2  4  3]); 

'/.  - 

prism. U_5b  =  reshape (tempU2 ,  prm.pts,  prm.pts,  prm. prism. number_angles, 
prm. lambda_bands) ; 

prism. U_5b  =  permute (prism. U_5b,  [1  2  4  3]); 

'/.  - 

prism. U_5c  =  reshape (tempU3 ,  prm.pts,  prm.pts,  prm. prism. number _ angles , 
prm . lambda_bands) ; 

prism. U_5c  =  permute (prism. U_5c,  [1  2  4  3]); 

7.  - 

prism. U_5d  =  reshape (tempU4,  prm.pts,  prm.pts,  prm. prism. number_angles, 
prm . lambda_bands) ; 

prism. U_5d  =  permute (prism. U_5d,  [1  2  4  3]); 

7.  - 

prism. U_5e  =  reshape (tempU5 ,  prm.pts,  prm.pts,  prm. prism. number _angles , 
prm . lambda_bands) ; 

prism. U_5e  =  permute (prism. U_5e,  [1  2  4  3]); 

•/.  - 

U_end  =  prism.U_5e; 

'/.  - 

clear  temp*; 

'/.  - 


'/.  - 

7.  PART  SIX 

7.  - 

7.  TRANSLATION  effects  upon  the  IR  electromagnetic  radiation 
7.  - 

7,  -  apply  the  translation  effects  due  to  propagation  through  the  direct  vision  prism 
7.  -  apply  coordinate  transfer  from  PRISM  coordinates  to  SYSTEM  coordinates 
7.  -  located  at  the  focussing  lens  of  the  IR  detector 

7.  - 


'/.  - 

7,  clockwise  (CW)  rotation  ==> 

7,  pol2cart(  kron(  -(prm.  prism.  rotate_angles)  .  5  +  pi/2,  ... 

7.  ones(  (prm.pts  *  prm.pts),  prm. lambda_bands  )  ),  prism. y_total_of f set  ); 
7.  counterclockwise  (CCW)  rotation  ==> 

7,  pol2cart(  kron(  (prm.  prism.  rotate_angles)  .  1  +  pi/2,  ... 

7,  ones(  (prm.pts  *  prm.pts),  prm. lambda_bands  )  ),  prism. y_total_of f set  ); 

'/.  - 


7.  - 

7.  compute  the  offset  displacements  for  each  discrete  ordinate  of  sampled 
7,  electromagnetic  radiation 

7.  - 

prism. y_total_off set  =  y4_total  -  repmat (yl_total ,  1,  prm. lambda_bands) ; 
prism. y_tilt  =  y5_total  -  y4_total; 

7.  - 

temp_angles  =  kron(  (prm. prism. rotate_angles) . ’  +  pi/2,  ... 

ones(  (prm.pts  *  prm.pts),  prm . lambda_bands  )  ); 

[  temp_x5_of f set  temp_y5_of f set  ]  =  ... 

pol2cart (temp_angles ,  prism. y_total_of f set  +  prism. y_tilt) ; 

7.  - 


7.  - 

7,  compute  the  ordinate  transfer  from  prism  coordinates  to 
7,  system  coordinates  for  each  discrete  ordinate  of  sampled 
7.  electromagnetic  radiation 

7.  - 

temp_system_x5_index  =  repmat(  reshape (prm. system. x4_index,  (prm.pts  *  prm.pts),  1) , 
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prm. prism. number_angles ,  prm. lambda_bands  ); 
temp_system_x5_ index  =  temp_system_x5_index  +  temp_x5_of f set ; 

'/.  - 

temp_system_y5_ index  =  repmat(  reshape (prm. system. y4_index,  (prm.pts  *  prm.pts) ,  1) , 
prm. prism. number_angles ,  prm. lambda_bands  ); 
temp_system_y5_index  =  temp_system_y5_index  +  temp_y5_of f set ; 

'/.  - 

temp_system_x5_index  =  reshape (temp_system_x5_index,  prm.pts,  prm.pts,  ... 

prm . prism . number_angles ,  prm . lambda_bands) ; 
prm. system. x5_index  =  permute (temp_system_x5_index,  [1  2  4  3]); 

7.  - 

temp_system_y5_index  =  reshape (temp_system_y5_index,  prm.pts,  prm.pts,  ... 

prm . prism . number_angles ,  prm . lambda_bands) ; 
prm. system. y5_index  =  permute (temp_system_y5_index,  [1  2  4  3]); 

'/.  - 


7.  - 

7,  compute  the  spatial  intervals  between  samples  ordinates  in 
%  system  coordinates  for  each  discrete  ordinate  of  sampled 
7,  electromagnetic  radiation 

7.  - 

7,  s2  ==>  angle  of  rotation  for  the  direct  vision  prism 
7.  - 

for  s2  =  1  :  prm. prism. number_angles; 

7.  - 

7.  si  ==>  spectral  wavelength  bands  of  electromagnetic  radiation 

7.  - 

for  si  =  1  :  prm . lambda_bands 

7.  - 

temp  =  []  ; 

temp  =  diff (prm. system. y5_index(: ,: ,sl,s2) ,[] ,1) ; 
prm. system. dy5 (si ,s2)  =  abs (temp(l , 1) ) ; 

7.  - 

temp  =  []  ; 

temp  =  diff (prm. system. x5_index( : , : ,sl,s2) ,  [] ,2) ; 
prm. system. dx5 (si ,s2)  =  abs (temp (1 , 1) ) ; 

7,  - 

end 

•/.  - 

end 

•/.  - 


7. 

return 
•/.  — • 
•/,  END 
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C.2.5  Propagate  from  Detector  Focussing  Lens  to  Focal  Plane  Array. 


location. 


C.2.5. 1  Search  for  the  actual  peak  intensity  in  vicinity  of  anticipated  peak  intensity 


function [U_end,  prm]  =  propagate_05_IR_array_01_search(prm,  U_start) ; 

7.  - 

7,  propagate_05_IR_array_01_search.m 

•/.  - 

7,  -  computes  the  propagation  of  the  electric  field  from  the 

7,  IR  DETECTOR  LENS  PLANE  to  the  IR  DETECTOR  InSb  ARRAY  PLANE. 

7.  -  computes  the  propagation  of  the  electric  field  (denoted  U)  , 

7.  starting  at  the  transverse  plane  where  the  electric  field 

7.  is  denoted  U_start,  and  ending  at  the  transverse  plane  where 

7,  the  electric  field  is  denoted  U_end. 

7.  - 

7,  EENG799  (Independent  Study)  Hyperspectral  Chromo tomography 
7,  Date:  16  January  2004 
7.  FLTLT  Anthony  J .  Dearinger 
7,  Royal  Australian  Air  Force 
7.  GE-04M 

7.  - 


7.  - 

7,  START 
7.  - 


7.  - 

sprintf  ( ’ - 5  ) 

sprintf ( ’Propagation  from  the  IR  DETECTOR  FOCUSSING  LENS  PLANE  to’) 
sprintf (’IR  DETECTOR  InSb  ARRAY  PLANE’) 

sprintf  (  ’ - ’  ) 


sprintf ( ’Computation  has  commenced  for  all  wavelengths  and  prism  rotations’) 

7.  - 


7.  - 

7,  establish  some  variables  of  interest 

7.  - 

clear  temp* 

7.  - 

tempX_start  =  reshape (prm. system. x5_index,  ... 

prm.pts,  (prm.pts  *  prm. lambda_bands  *  prm. prism. number _ angles) ) ; 
tempY_start  =  reshape (prm. system. y5_index,  ... 

prm.pts,  (prm.pts  *  prm. lambda_bands  *  prm. prism. number_angles) ) ; 

7.  - 


7.  - 

7,  compute  the  lens  thickness  function 

7,  -  thickness  of  the  lens  across  the  transverse  plane  to  the 

7,  optic  axis  of  the  system 

7.  - 

lens_thick_detector  =  [] ; 

7  - 

lens_thick_detector  =  prm.delta_0_detector  -  ... 

(tempX_start . ~2  +  tempY_start . "2)  *  ... 

(  (prm.Rl_detector)  /  (2*(prm.Rl_detector) "2)  -  ... 

(prm.R2_detector)  /  (2* (prm.R2_detector) "2)  ); 

7.  - 

prm. lens_thick_detector  =  lens_thick_detector ; 

7.  - 


7.  - 

7,  re-establish  some  variables  of  interest 
7.  - 
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clear  temp* 

•/.  - 

tempX_start  =  prm. system. x5_index; 
tempY_start  =  prm. system. y5_index; 
tempZ  =  prm. f ocal_detector ; 

'/.  - 


7.  - 

7,  establish  the  IR  radiation  field,  U_start,  and  the  wavenumber  parameter, 

7*  tempK,  representing  the  wavenumber  for  each  wavelength  band  at  the 
7,  centre  wavelength  for  each  wavelength  band 

•/.  - 

U_start  =  reshape (U_start,  prm.pts,  ... 

(prm.pts  *  prm. lambda_bands  *  prm. prism. number_angles) ) ; 

'/.  - 

tempK  =  ones ( (prm.pts  *  prm.pts  *  prm. prism. number_angles) ,  1)  *  prm.k; 

tempK  =  reshape (tempK,  prm.pts,  prm.pts,  prm. prism. number_angles ,  prm. lambda_bands) ; 

tempK  =  permute (tempK,  [1  2  4  3]); 

•/.  - 


'/.  - 

7,  apply  the  PHASE  effects  introduced  by  the  lens  thickness  function 
7,  to  the  optical  IR  radiation  field  propagating  through  the  lens, 

7,  for  each  wavelength  band  of  the  optical  IR  electromagnetic  radiation 

7.  - 

U_start  =  U_start  .*  exp(  -j  *  (prm.n_detector  -  1)  *  lens_thick_detector  .*  ... 

reshape (tempK,  prm.pts,  (prm.pts  *  prm . lambda_bands  *  prm. prism. number_angles) )  ); 

'/.  - 

U_start  =  reshape (U_start,  prm.pts,  prm.pts,  prm. lambda_bands ,  ... 

prm.prism.number_angles) ; 

'/.  - 


7.  - 

clear  lens_thick_detector  Xcentre  Ycentre 

7.  - 


'/.  - 

7,  compute  the  centre  for  each  location  where  point  spread  function  (PSF) 
7,  is  expected  to  be  imaged  on  the  IR  detector  focal  plane  array 

•/.  - 

7,  -  the  centre  for  each  PSF  is  based  on: 

•/.  - 

7,  —  each  rotation  angle  of  the  direct  vision  prism 

7,  —  each  wavelength  band  represented  by  each  centre  wavelength 

'/.  - 

Xcentre  =  -  tempZ  *  sin (prm. prism. epsilon_4bX_CW) ; 

Ycentre  =  -  tempZ  *  sin(prm. prism. epsilon_4bY) ; 

7.  - 


'/.  - 

7,  establish  some  variables  of  interest 

'/.  - 

dx  =  prm. system. dx6; 
dy  =  prm. system. dy6; 

7.  - 

Xsamples  =  20; 

Ysamples  =  20; 

7.  - 

tempXsamples_total  =  5  *  (-Xsamples)  :  5  *  (Xsamples) ; 
tempYsamples_total  =  5  *  (-Ysamples)  :  5  *  (Ysamples) ; 

7.  - 

tempXsamples_centre  =  5  *  (Xsamples)  +  1; 
tempYsamples_centre  =  5  *  (Ysamples)  +  1; 

7.  - 

prm. detector . Xcentre_temp  =  Xcentre; 
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prm. detector . Ycentre_temp  =  Ycentre; 

prm. detector . Xsamples_total  =  tempXsamples_total ; 

prm. detector . Ysamples_total  =  tempYsamples_total ; 

prm. detector . Xsample_distances  =  tempXsamples_total  *  dx; 

prm. detector . Ysample_distances  =  tempYsamples_total  *  dy; 

prm. detector .X_samples  =  length(prm. detector . Xsamples_total) ; 

prm. detector .Y_samples  =  length(prm. detector . Ysamples_total) ; 

prm. detector . Xsamples_centre  =  tempXsamples_centre ; 

prm. detector . Ysamples_centre  =  tempYsamples_centre ; 

'/.  - 


'/.  - 

7.  s2  ==>  angle  of  rotation  for  the  direct  vision  prism 

7.  - 

for  s2  =  1  :  prm. prism. number_angles ; 

•/.  - 

l  comments 

l  -  only  image  points  along  the  upright  prism  axis  of  dispersion 

•/.  - 

7#  si  ==>  spectral  bands  of  radiation 

•/.  - 

for  si  =  1  :  prm . lambda_bands 

•/.  - 

sX  =  0; 

7.  - 

7,  comments 

7,  -  entire  width  of  imaging  area  is  sX  =  [  -Xsamples  :  Xsamples  ] 

7.  - 

for  sY  =  tempYsamples_total 

•/.  - 

7,  comments 

7,  -  entire  height  of  imaging  area  is  sY  =  [  -Ysamples  :  Ysamples 
7,  -  search  for  Ycentre  is  being  conducted  along  a  length  that  is 
7.  5  times  the  size  of  the  imaging  area  height 

•/.  - 

7.  - 

7,  CCW  (counter-clockwise  rotation) 

•/.  - 

7,  temp_Xrotation  =  cos  (prm.  prism.  rotate_angles  (s2)  )  *  ... 

7,  prm. detector . Xsample_distances (sX  +  tempXsamples_centre)  - 

7,  sin(prm. prism. rotate_angles (s2) )  *  ... 

7,  prm.  detector  .Ysample_di stances  (sY  +  tempYsamples_centre)  ; 

7,  temp_Yrotation  =  sin(prm. prism. rotate_angles (s2) )  *  ... 

7,  prm.  detector  .Xsample_distances(sX  +  tempXsamples_centre)  + 

7,  cos  (prm.  prism.  rotate_angles  (s2) )  *  ... 

7,  prm.  detector  .Ysample_di stances  (sY  +  tempYsamples_centre)  ; 

7.  - 

7.  - 

7,  CW  (clockwise  rotation) 

7.  - 

temp_Xrotation  =  cos (prm. prism. rotate_angles (s2) )  *  ... 

prm. detector . Xsample_distances (sX  +  tempXsamples_centre)  + 
sin(prm. prism. rotate_angles (s2) )  *  ... 

prm. detector . Ysample_distances(sY  +  tempYsamples_centre) ; 
temp_Yrotation  =  -  sin(prm. prism. rotate_angles (s2) )  *  ... 

prm. detector . Xsample_distances (sX  +  tempXsamples_centre)  + 
cos (prm. prism. rotate_angles (s2) )  *  ... 

prm. detector . Ysample_distances(sY  +  tempYsamples_centre) ; 

7.  - 


7.  - 

tempX_end( (sY  +  tempYsamples_centre) ,  si,  s2)  = 
temp_Xrotation  +  Xcentre (si , s2) ; 
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7.  - 

tempY_end( (sY  +  tempYsamples_centre) ,  si,  s2)  = 
temp_Yrotation  +  Ycentre(sl,s2) ; 

•/.  - 


•/.  - 

tempX  =  []  ; 

tempX  =  tempX_end( (sY  +  tempYsamples_centre) ,  si,  s2)  - 
tempX_start ( : , : , si , s2) ; 

7.  - 

tempY  =  []  ; 

tempY  =  tempY_end( (sY  +  tempYsamples_centre) ,  si,  s2)  - 
tempY_start ( : , : , si , s2) ; 

7.  - 

prop_dist  =  []  ; 

prop_dist  =  sqrt(  tempX.  ~2  +  tempY.  "2  +  tempZ.~2  ); 

•/.  - 


•/.  - 

clear  tempU* 

7.  - 

tempU  =  squeeze (U_st art(: ,: ,sl,s2))  ./  prop_dist  . * 
exp(-j  *  prop_dist  . *  tempK( : , : , si , s2) ) ; 

7.  - 

tempU  =  reshape (tempU,  (prm.pts  *  prm.pts),  1); 
tempU  =  sum (tempU) ; 

%  - 

U_end((sY  +  tempYsamples_centre) ,  si,  s2)  =  tempU; 

7.  - 

end 

t  - 

end 

•/.  - 

sprintf  ( ’prism  angular  position  (s2)  =  °/0d  of  Id*  ,  ... 

s2,  prm. prism. number_angles) 

•/.  - 

end 

'/.  - 


'/.  - 

prm. detector .PSF_Xsamples_centre  =  tempX_end; 
prm. detector .PSF_Ysamples_centre  =  tempY_end; 

7.  - 


7.  - 

7,  find  the  peak  intensity  of  the  PSF  samples,  and  the  corresponding 
%  index  in  the  array,  along  the  PSF  centres  radial  for  all  spectral 

7,  wavelengths,  for  each  angular  position  of  prism  rotation 

•/.  - 

tempU_end  =  reshape (U_end,  length (tempYsamples_total) ,  ... 

prm. lambda_bands  *  prm. prism. number_angles) ; 

[tempU_end_peak  tempU_end_peak_index]  =  ... 

max( (tempU_end  . *  conj (tempU_end) ) , [] , 1) ; 
tempU_end_peak  =  reshape (tempU_end_peak, 

prm. lambda_bands ,  prm. prism. number_angles) ; 
tempU_end_peak_index  =  reshape (tempU_end_peak_index, 
prm . lambda_bands ,  prm. prism. number_angles) ; 

'/.  - 

prm. detector .U_6_peak  =  tempU_end_peak; 

prm. detector .U_6_peak_index  =  tempU_end_peak_ index ; 

•/.  - 


7.  - 

7,  use  the  corresponding  index  in  the  array  to  compute  the  offsets, 

7,  Xcentre_off set  and  Ycentre_of f set ,  and  thus  aid  with  designating 
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7,  the  centres  of  the  respective  PSFs  for  all  spectral  wavelengths, 

7.  for  each  angular  position  of  prism  rotation 

7.  - 

7,  s2  ==>  angle  of  rotation  for  the  direct  vision  prism 
7.  - 

for  s2  =  1  :  prm. prism. number_ angles ; 

7.  - 

7.  si  ==>  spectral  bands  of  radiation 

7.  - 

for  si  =  1  :  prm . lambda_bands 

7.  - 

prm. detector . PSF_Xcentre_of f set (si, s2)  =  squeeze (prm . detector . 
PSF_Xsamples_centre (tempU_end_peak_index(sl , s2) ,sl,s2))  - 
Xcentre(sl,s2) ; 

prm. detector . PSF_Ycentre_of f set (si, s2)  =  squeeze (prm . detector . 
PSF_Ysamples_centre (tempU_end_peak_index(sl , s2) ,sl,s2))  - 
Ycentre(sl,s2) ; 

7.  - 

end 

7.  - 

end 

7.  - 


7.  - 

7,  assign  the  ordinates,  Xcentre  and  Ycentre,  as  centres  of  the  respective 
7,  PSFs  for  all  spectral  wavelengths,  for  each  angular  position  of 
7,  prism  rotation 

7.  - 

prm. detector .PSF_Xcentre  =  Xcentre  +  prm. detector .PSF_Xcentre_off set ; 
prm. detector .PSF_Ycentre  =  Ycentre  +  prm. detector .PSF_Ycentre_off set ; 

y.  - 

prm. detector .PSF_FPA_Xcentre  =  round (prm . detector . PSF_Xcentre  /  prm. system. dx6)  * 
prm . system . dx6 ; 

prm. detector .PSF_FPA_Ycentre  =  round (prm . detector . PSF_Ycentre  /  prm. system. dy6)  * 
prm . system . dy6 ; 

7.  - 


7.  - 

clear  temp*  prop_dist; 

7.  - 

sprintf ( ’Computation  has  completed  for  all  wavelengths  and  prism  rotations’) 
sprintf  (  ’ - ’  ) 

7.  - 


7. 

return 

7.  - 

7,  END 
7.  — - 
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C.2.5.2  Compute  field  of  propagation  at  Focal  Plane  Array  around  actual  peak  intensity 

location. 

function [U_end,  prm]  =  propagate_05_IR_array_03_image (prm,  U_start) ; 

7.  - 

7,  propagate_05_IR_array_03_image  .m 

7.  - 

7,  -  computes  the  propagation  of  the  electric  field  from  the 

7.  IR  DETECTOR  LENS  PLANE  to  the  IR  DETECTOR  InSb  ARRAY  PLANE. 

7.  -  computes  the  propagation  of  the  electric  field  (denoted  U)  , 

7,  starting  at  the  transverse  plane  where  the  electric  field 

7,  is  denoted  U_start,  and  ending  at  the  transverse  plane  where 

7.  the  electric  field  is  denoted  U_end. 

7.  - 

7,  EENG799  (Independent  Study)  Hyperspectral  Chromo tomography 
7,  Date:  06  January  2004 
7,  FLTLT  Anthony  J.  Dearinger 
7,  Royal  Australian  Air  Force 
7.  GE-04M 

7.  - 


7.  - 

7.  START 
7.  - 


7.  - 

sprintf  (  ’ - ’  ) 

sprintf ( ’Propagation  from  the  IR  DETECTOR  FOCUSSING  LENS  PLANE  to’) 
sprintf (’IR  DETECTOR  InSb  ARRAY  PLANE’) 

sprintf  (  ’ - ’  ) 


sprintf ( ’Computation  has  commenced  for  all  wavelengths  and  prism  rotations’) 

7.  - 


7.  - 

7,  establish  some  variables  of  interest 

7.  - 

clear  temp* 

7.  - 

tempX_start  =  reshape (prm. system. x5_index,  ... 

prm.pts,  (prm.pts  *  prm. lambda_bands  *  prm. prism. number _angles) ) ; 
tempY_start  =  reshape (prm. system. y5_index,  ... 

prm.pts,  (prm.pts  *  prm. lambda_bands  *  prm. prism. number_angles) ) ; 

7.  - 

dx  =  prm. system. dx6; 
dy  =  prm. system. dy6; 

7.  - 

Xsamples  =  20; 

Ysamples  =  20; 

Xsamples_centre  =  Xsamples  +  1; 

Ysamples_centre  =  Ysamples  +  1; 

y#  - 


y#  - 

7,  compute  the  lens  thickness  function 

7,  -  thickness  of  the  lens  across  the  transverse  plane  to  the 

7.  optic  axis  of  the  system 

y#  - 

lens_thick_detector  =  [] ; 

y#  - 

lens_thick_detector  =  prm.delta_0_detector  -  ... 

(tempX_start . ~2  +  tempY_start . "2)  *  ... 

(  (prm.Rl_detector)  /  (2*(prm.Rl_detector) "2)  -  ... 

(prm.R2_detector)  /  (2* (prm.R2_detector) ~2)  ); 

y#  - 

prm. lens_thick_detector  =  lens_thick_detector ; 
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7. 


7.  - 

%  re-establish  some  variables  of  interest 

7.  - 

clear  temp* 

7.  - 

tempX_start  =  prm. system. x5_index; 
tempY_start  =  prm. system. y5_index; 
tempZ  =  prm. f ocal_detector ; 

7.  - 


7.  - 

%  establish  the  IR  radiation  field,  U_start,  and  the  wavenumber  parameter, 

7*  tempK,  representing  the  wavenumber  for  each  wavelength  band  at  the 
7,  centre  wavelength  for  each  wavelength  band 

'/.  - 

U_start  =  reshape (U_start,  prm.pts,  ... 

(prm.pts  *  prm. lambda_bands  *  prm. prism. number_angles) ) ; 

7.  - 

tempK  =  ones ( (prm.pts  *  prm.pts  *  prm. prism. number_angles) ,  1)  *  prm.k; 

tempK  =  reshape (tempK,  prm.pts,  prm.pts,  prm. prism. number_angles ,  prm. lambda_bands) ; 

tempK  =  permute (tempK,  [1  2  4  3]); 

'/.  - 


7.  - 

7.  apply  the  PHASE  effects  introduced  by  the  lens  thickness  function 
7,  to  the  optical  IR  radiation  field  propagating  through  the  lens, 

7,  for  each  wavelength  band  of  the  optical  IR  electromagnetic  radiation 

7.  - 

U_start  =  U_start  .*  exp(  -j  *  (prm.n_detector  -  1)  *  lens_thick_detector  .*  ... 

reshape (tempK,  prm.pts,  (prm.pts  *  prm . lambda_bands  *  prm. prism. number_angles) )  ); 

7.  - 

U_start  =  reshape (U_start,  prm.pts,  prm.pts,  prm. lambda_bands ,  ... 

prm.prism.number_angles) ; 

7.  - 


7.  - 

clear  lens_thick_detector  Xcentre  Ycentre 

7.  - 

7.  7.  - 

7,  7.  compute  the  centre  for  each  location  where  point  spread  function  (PSF) 
7,  7o  is  expected  to  be  imaged  on  the  IR  detector  focal  plane  array 

7.  7.  - 

7,  7.  -  the  centre  for  each  PSF  is  based  on: 

7.  7,  - 

7.  7o  —  each  rotation  angle  of  the  direct  vision  prism 

7,  7o  —  each  wavelength  band  represented  by  each  centre  wavelength 

7.  7.  - 

7,  Xcentre  =  -  tempZ  *  sin(prm. prism. epsilon_4bX_CW)  ; 

7,  Ycentre  =  -  tempZ  *  sin(prm. prism. epsilon_4bY)  ; 

7.  7,  - 


7.  - 

7,  establish  some  variables  of  interest 

7.  - 

prm. detector .Xsamples  =  [] ; 

prm. detector .Xsamples  =  [-Xsamples  :  Xsamples] ; 
prm. detector .Ysamples  =  [] ; 

prm. detector .Ysamples  =  [-Ysamples  :  Ysamples] ; 
prm. detector . Xsample_distances  =  [] ; 

prm. detector . Xsample_distances  =  [-Xsamples  :  Xsamples]  *  dx; 
prm. detector . Ysample_distances  =  [] ; 

prm. detector . Ysample_distances  =  [-Ysamples  :  Ysamples]  *  dy; 
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•/. 


7,  - 

for  s2  =  1  :  prm. prism. number_aiigles 

i  - 

for  si  =  1  :  prm. lambda_bands 

7,  - 

for  sX  =  prm. detector. Xsamples 

7.  - 

for  sY  =  prm. detector . Ysamples 

y.  - 

tempX_end  =  []  ; 

tempX_end  =  (sX  *  dx)  +  prm. detector. PSF_FPA_Xcentre (si, s2) ; 

7.  - 

tempY_end  =  []  ; 

tempY_end  =  (sY  *  dy)  +  prm. detector. PSF_FPA_Ycentre (si, s2) ; 

7.  - 

tempX  =  []  ; 

tempX  =  tempX_end  -  tempX_start ( : , : , si , s2) ; 

7.  - 

tempY  =  []  ; 

tempY  =  tempY_end  -  tempY_start ( : , : , si , s2) ; 

7.  - 

prop_dist  =  []  ; 

prop_dist  =  sqrt(  tempX. "2  +  tempY. ~2  +  tempZ. ~2  ); 

7,  - 


7.  - 

clear  temp_U* 

7.  - 

7.  - 

7*  expression  A  -  approximated  Rayleigh-Sommerf ield  diffraction  propagation 
7,  expression  A  -  was  originally  used 

7*  - 

7,  temp_U  =  squeeze  (U_start  (:,:  ,sl,s2))  ./  prop_dist  .  *  ... 

7,  exp(-j  *  prop_dist  .  *  tempK(  :  ,  :  ,  si ,  s2)  )  ; 

7.  - 

7.  - 

7,  expression  B  -  approximated  Rayleigh-Sommerf ield  diffraction  propagation 
7.  expression  B  -  more  accurate  numerical  approximation  compared  to 
7,  propagation  using  expression  A 

7.  - 

temp_U_amplitude  =  ... 

sinc(  (sX  *  dx)  *  pi  *  prm. system. dx5 (si , s2)  ./  tempZ  ./  ... 

prm.lambda(sl)  )  .*  ... 

sinc(  (sY  *  dy)  *  pi  *  prm. system. dy5 (si , s2)  ./  tempZ  ./  ... 

prm. lambda (si)  )  .*  ... 

prm. system. dx5 (si, s2)  .*  prm. system. dy5 (si, s2)  .*  tempZ  ./  ... 

j  ./  prm. lambda (si) ; 

7.  - 

temp_U  =  U_start ( : , : ,sl,s2)  .*  temp_U_ amplitude  ./  prop_dist . ~2  .*  ... 

exp(-j  *  prm.k(sl)  .*  prop_dist) ; 

7*  - 


7,  - 

temp_U  =  reshape (temp_U,  (prm.pts  *  prm.pts),  1); 
temp_U  =  sum(temp_U) ; 

7.  - 

U_end((sY  +  Ysamples_centre) , (sX  +  Xsamples_centre) , si , s2)  =  temp_U; 

7.  - 

end 

•/#  - 

end 


140 


#/o  - 

end 

•/.  - 

sprintf  ( ’prism  angular  position  (s2)  =  "/,d  of  ’/,d’ , 
s2,  prm. prism. number_angles) 

’/.  - 

end 

•/.  - 


l  - 

clear  temp*  prop_dist; 

•/.  - 

sprintf (’ Computation  has  completed  for  all  wavelengths  and  prism  rotations’) 

sprintf  ( ’ - ’ ) 

•/.  - 


•/. 

return 

7.  — - 

'/.  END 
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